List of GKV output

GKV output files are:

  • The output directory DIR/

    • cnt/*cnt*

    • fxv/*fxv*

    • phi/\*phi\*, \*Al\*, \*mom\*, \*trn\*, (\*tri\* for nonlinear runs)

    • hst/\*bln\*, \*geq\*, \*gem\*, \*qes\*, \*qem\*, \*wes\*, \*wem\*, \*eng\*, \*men\*, \*dtc\*, \*mtr\*, (\*frq\*, \*dsp\* for linear runs)

    • log/*log*

Their explanations are summarized below.

Explanations on GKV output files

Table 2 cnt/gkvp.cnt.(inum in 3 digits).zarr/

Field

Description

File type

Zarr format binary

Output timing

End of the run

MPI ranks

All

Total files

nprocw*nprocz*nprocv*nprocm*nprocs*(Total run numbers)

GKV unit

ocnt

Stored data

Coordinates:

  • t Simulation time

  • is particle species (integer)

  • mu magnetic moment (or vp perpendicular velocity)

  • vl parallel velocity

  • zz field-aligned coordinate

  • ky Field-line-label (poloidal) wavwenumber

  • kx radial wavenumber

Data variables: cnt[t,is,mu,vl,zz,ky,kx] Perturbed distribution function \(\tilde{f}_{\mathrm{s}\mathbf{k}}\) (double complex).

Table 3 fxv/gkvp.fxv.(inum in 3 digits).zarr/

Field

Description

File type

Zarr format binary

Output timing

dtout_fxv

MPI ranks

All

Total files

nprocw*nprocz*nprocv*nprocm*nprocs*(Total run numbers)

GKV unit

ofxv

Stored data

Coordinates:

  • t Simulation time,

  • is particle species (integer)

  • mu magnetic moment (or vp perpendicular velocity)

  • vl parallel velocity, zz field-aligned coordinate

  • ky Field-line-label (poloidal) wavwenumber

  • kx radial wavenumber

Data variables: fxv[t,is,mu,vl,zz,ky,kx] Perturbed distribution function \(\tilde{f}_{\mathrm{s}\mathbf{k}}\) (double complex) at iz = -nz in each rankz .

Table 4 phi/gkvp.phi.(inum in 3 digits).zarr/

Field

Description

File type

Zarr format binary

Output timing

dtout_ptn

MPI ranks

ranks == 0 .and. vel_rank == 0

Total files

nprocw*nprocz*(Total run numbers)

GKV unit

ophi

Stored data

Coordinates:

  • t Simulation time

  • zz field-aligned coordinate

  • ky Field-line-label (poloidal) wavwenumber

  • kx radial wavenumber

Data variables: phi[t,zz,ky,kx] Perturbed electrostatic potential \(\tilde{\phi}_{\mathbf{k}}\) .

Table 5 phi/gkvp.Al.(inum in 3 digits).zarr/

Field

Description

File type

Zarr format binary

Output timing

dtout_ptn

MPI ranks

ranks == 0 .and. vel_rank == 0

Total files

nprocw*nprocz*(Total run numbers)

GKV unit

oAl

Stored data

Coordinates:

  • t Simulation time

  • zz field-aligned coordinate

  • ky Field-line-label (poloidal) wavwenumber

  • kx radial wavenumber

Data variables: Al[t,zz,ky,kx] Perturbed vector potential \(\tilde{A}_{\parallel\mathbf{k}}\) .

Table 6 phi/gkvp.mom.(inum in 3 digits).zarr/

Field

Description

File type

Zarr format binary

Output timing

dtout_ptn

MPI ranks

vel_rank == 0

Total files

nprocw*nprocz*nprocs*(Total run numbers)

GKV unit

omom

Stored data

Coordinates:

  • t Simulation time

  • is particle species (integer)

  • imom index of fluid moments (integer)

  • zz field-aligned coordinate

  • ky Field-line-label (poloidal) wavwenumber

  • kx radial wavenumber

Data variables: mom[t,is,imom,zz,ky,kx] Perturbed fluid moment (double complex). In the present version nmom = 6 :

  • \(\tilde{n}_{\mathrm{s}\mathbf{k}} = \int dv^3\, J_{0\mathrm{s}\mathbf{k}}\, \tilde{f}_{\mathrm{s}\mathbf{k}}\)

  • \(\tilde{u}_{\parallel\mathrm{s}\mathbf{k}} = \int dv^3\, v_\parallel\, J_{0\mathrm{s}\mathbf{k}}\, \tilde{f}_{\mathrm{s}\mathbf{k}}\)

  • \(\tilde{p}_{\parallel\mathrm{s}\mathbf{k}} = \int dv^3\, \frac{m_\mathrm{s} v_\parallel^2}{2}\, J_{0\mathrm{s}\mathbf{k}}\, \tilde{f}_{\mathrm{s}\mathbf{k}}\)

  • \(\tilde{p}_{\perp\mathrm{s}\mathbf{k}} = \int dv^3\, \mu B\, J_{0\mathrm{s}\mathbf{k}}\, \tilde{f}_{\mathrm{s}\mathbf{k}}\)

  • \(\tilde{q}_{\parallel\parallel\mathrm{s}\mathbf{k}} = \int dv^3\, v_\parallel \frac{m_\mathrm{s} v_\parallel^2}{2}\, J_{0\mathrm{s}\mathbf{k}}\, \tilde{f}_{\mathrm{s}\mathbf{k}}\)

  • \(\tilde{q}_{\parallel\perp\mathrm{s}\mathbf{k}} = \int dv^3\, v_\parallel \mu B\, J_{0\mathrm{s}\mathbf{k}}\, \tilde{f}_{\mathrm{s}\mathbf{k}}\)

  • Normalized by \(\delta_\mathrm{ref}n_\mathrm{ref}\) , \(\delta_\mathrm{ref}n_\mathrm{ref}v_\mathrm{ref}\) , \(\delta_\mathrm{ref}n_\mathrm{ref}T_\mathrm{ref}\) , \(\delta_\mathrm{ref}n_\mathrm{ref}T_\mathrm{ref}\) , \(\delta_\mathrm{ref}n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}\) , \(\delta_\mathrm{ref}n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}\) , respectively.

Table 7 phi/gkvp.trn.(inum in 3 digits).zarr/

Field

Description

File type

Zarr format binary

Output timing

dtout_eng

MPI ranks

zsp_rank == 0 .and. vel_rank == 0

Total files

nprocw*nprocs*(Total run numbers)

GKV unit

otrn

Stored data

Coordinates:

  • t Simulation time

  • is particle species (integer)

  • itrn index of entropy balance diagnostics (integer)

  • ky Field-line-label (poloidal) wavwenumber

  • kx radial wavenumber

Data variables: trn[t,is,itrn,ky,kx]

  • \(S_{\mathrm{s}\mathbf{k}}(-nx:nx,0:ny)\) : Perturbed gyrocenter entropy \([\delta_\mathrm{ref}^2 n_\mathrm{ref} T_\mathrm{ref}]\) (real*8)

  • \(W_{\mathrm{E}\mathbf{k}}(-nx:nx,0:ny)\) : Electrostatic field energy including polarization \([\delta_\mathrm{ref}^2 n_\mathrm{ref} T_\mathrm{ref}]\) (real*8)

  • \(W_{\mathrm{M}\mathbf{k}}(-nx:nx,0:ny)\) : Magnetic field energy \([\delta_\mathrm{ref}^2 n_\mathrm{ref} T_\mathrm{ref}]\) (real*8)

  • \(R_{\mathrm{sE}\mathbf{k}}(-nx:nx,0:ny)\) : Wave–particle interaction \(( W_{\mathrm{E}\mathbf{k}} \rightarrow S_{\mathrm{s}\mathbf{k}} )\) \([\delta_\mathrm{ref}^2 n_\mathrm{ref} T_\mathrm{ref} v_\mathrm{ref}/L_\mathrm{ref}]\) (real*8)

  • \(R_{\mathrm{sM}\mathbf{k}}(-nx:nx,0:ny)\) : Wave–particle interaction \(( W_{\mathrm{M}\mathbf{k}} \rightarrow S_{\mathrm{s}\mathbf{k}} )\) \([\delta_\mathrm{ref}^2 n_\mathrm{ref} T_\mathrm{ref} v_\mathrm{ref}/L_\mathrm{ref}]\) (real*8)

  • \(I_{\mathrm{sE}\mathbf{k}}(-nx:nx,0:ny)\) : Nonlinear entropy transfer by \(\mathbf{E}\times\mathbf{B}\) flow \([\delta_\mathrm{ref}^2 n_\mathrm{ref} T_\mathrm{ref} v_\mathrm{ref}/L_\mathrm{ref}]\) (real*8)

  • \(I_{\mathrm{sM}\mathbf{k}}(-nx:nx,0:ny)\) : Nonlinear entropy transfer by magnetic flutter \([\delta_\mathrm{ref}^2 n_\mathrm{ref} T_\mathrm{ref} v_\mathrm{ref}/L_\mathrm{ref}]\) (real*8)

  • \(D_{\mathrm{s}\mathbf{k}}(-nx:nx,0:ny)\) : Collisional dissipation \([\delta_\mathrm{ref}^2 n_\mathrm{ref} T_\mathrm{ref} v_\mathrm{ref}/L_\mathrm{ref}]\) (real*8)

  • \(\Gamma_{\mathrm{sE}\mathbf{k}}(-nx:nx,0:ny)\) : Particle flux by \(\mathbf{E}\times\mathbf{B}\) flow \([\delta_\mathrm{ref}^2 n_\mathrm{ref} v_\mathrm{ref}]\) (real*8)

  • \(\Gamma_{\mathrm{sM}\mathbf{k}}(-nx:nx,0:ny)\) : Particle flux by magnetic flutter \([\delta_\mathrm{ref}^2 n_\mathrm{ref} v_\mathrm{ref}]\) (real*8)

  • \(Q_{\mathrm{sE}\mathbf{k}}(-nx:nx,0:ny)\) : Energy flux by \(\mathbf{E}\times\mathbf{B}\) flow \([\delta_\mathrm{ref}^2 n_\mathrm{ref} T_\mathrm{ref} v_\mathrm{ref}]\) (real*8)

  • \(Q_{\mathrm{sM}\mathbf{k}}(-nx:nx,0:ny)\) : Energy flux by magnetic flutter \([\delta_\mathrm{ref}^2 n_\mathrm{ref} T_\mathrm{ref} v_\mathrm{ref}]\) (real*8)

See also Supplemental Entropy balance equation for each wavenumber and plasma species .

Table 8 phi/gkvp.tri.mx(mxt in 4 digits)my(myt in 4 digits).(inum in 3 digits).zarr/

Field

Description

File type

Zarr format binary

Output timing

dtout_ptn (when calc_type == "nonlinear" and num_triad_diag > 0 )

MPI ranks

rank == 0

Total files

nprocs*num_triad_diag*(Total run numbers)

GKV unit

otri

Stored data

Coordinates:

  • t Simulation time

  • is particle species (integer)

  • itri index of triad transfer diagnostics (integer)

  • ky Field-line-label (poloidal) wavwenumber

  • kx radial wavenumber

Data variables: tri[t,is,itri,ky,kx]

  • \(J_{\mathrm{sE}\mathbf{k}}^{\mathbf{p,q}}(-nx:nx,-global_{y}:global_{y})\) : Triad transfer function from modes \(\mathbf{p}, \mathbf{q}\) to mode \(\mathbf{k}\) via \(\mathbf{E}\times\mathbf{B}\) nonlinearity \([\delta_\mathrm{ref}^2 n_\mathrm{ref} T_\mathrm{ref} v_\mathrm{ref}/L_\mathrm{ref}]\) (real*8)

  • \(J_{\mathrm{sE}\mathbf{p}}^{\mathbf{q,k}}(-nx:nx,-global_{y}:global_{y})\) : Cyclic change \((\mathbf{k,p,q}) \rightarrow (\mathbf{p,q,k})\) (real*8)

  • \(J_{\mathrm{sE}\mathbf{q}}^{\mathbf{k,p}}(-nx:nx,-global_{y}:global_{y})\) : Cyclic change \((\mathbf{p,q,k}) \rightarrow (\mathbf{q,k,p})\) (real*8)

  • \(J_{\mathrm{sM}\mathbf{k}}^{\mathbf{p,q}}(-nx:nx,-global_{y}:global_{y})\) : Triad transfer function from modes \(\mathbf{p}, \mathbf{q}\) to mode \(\mathbf{k}\) via magnetic-flutter nonlinearity \([\delta_\mathrm{ref}^2 n_\mathrm{ref} T_\mathrm{ref} v_\mathrm{ref}/L_\mathrm{ref}]\) (real*8)

  • \(J_{\mathrm{sM}\mathbf{p}}^{\mathbf{q,k}}(-nx:nx,-global_{y}:global_{y})\) : Cyclic change \((\mathbf{k,p,q}) \rightarrow (\mathbf{p,q,k})\) (real*8)

  • \(J_{\mathrm{sM}\mathbf{q}}^{\mathbf{k,p}}(-nx:nx,-global_{y}:global_{y})\) : Cyclic change \((\mathbf{p,q,k}) \rightarrow (\mathbf{q,k,p})\) (real*8)

Diagnosed for a fixed mode \(\mathbf{k} = (\texttt{mxt,myt})\) and plotted as a 2D function of \(\mathbf{p} = (p_x,p_y)\) , where the triad condition determines \(\mathbf{q} = -\mathbf{k}-\mathbf{p}\) . See also Supplemental Triad transfer function .

Table 9 hst/gkvp_f0.48.bln.(ranks in 1 digits).(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

dtout_eng

MPI ranks

rank == 0

Total files

nprocs*(Total run numbers)

GKV unit

obln

Stored data

time, \(S_{\mathrm{s}}\) , \(W_{\mathrm{E}}\) , \(W_{\mathrm{M}}\) , \(R_{\mathrm{sE}}\) , \(R_{\mathrm{sM}}\) , \(I_{\mathrm{sE}}\) , \(I_{\mathrm{sM}}\) , \(D_{\mathrm{s}}\) , \(\frac{T_\mathrm{s}\Gamma_{\mathrm{sE}}}{L_{p\mathrm{s}}}\) , \(\frac{T_\mathrm{s}\Gamma_{\mathrm{sM}}}{L_{p\mathrm{s}}}\) , \(\frac{\Theta_{\mathrm{sE}}}{L_{T\mathrm{s}}}\) , \(\frac{\Theta_{\mathrm{sM}}}{L_{T\mathrm{s}}}\)

where,

  • time : Simulation time \(t\) [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real*8)

  • \(S_{\mathrm{s}}\) (0:1): Perturbed gyrocenter entropy [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}\) ] (real*8).

  • \(W_{\mathrm{E}}\) (0:1): Electrostatic field energy including polarization [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}\) ] (real*8).

  • \(W_{\mathrm{M}}\) (0:1): Magnetic field energy [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}\) ] (real*8).

  • \(R_{\mathrm{sE}}\) (0:1): Wave-particle interaction ( \(W_{\mathrm{E}\bm{k}} \rightarrow S_{\mathrm{s}\bm{k}}\) ) [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}/L_\mathrm{ref}\) ] (real*8).

  • \(R_{\mathrm{sM}}\) (0:1): Wave-particle interaction ( \(W_{\mathrm{M}\bm{k}} \rightarrow S_{\mathrm{s}\bm{k}}\) ) [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}/L_\mathrm{ref}\) ] (real*8).

  • \(I_{\mathrm{sE}}\) (0:1): Nonlinear entropy transfer by \(\bm{E}\times\bm{B}\) flow [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}/L_\mathrm{ref}\) ] (real*8).

  • \(I_{\mathrm{sM}}\) (0:1): Nonlinear entropy transfer by magnetic flutter [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}/L_\mathrm{ref}\) ] (real*8).

  • \(D_{\mathrm{s}}\) (0:1): Collisional dissipation [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}/L_\mathrm{ref}\) ] (real*8).

  • \(\frac{T_\mathrm{s}\Gamma_{\mathrm{sE}}}{L_{p\mathrm{s}}}\) : Particle flux term by \(\bm{E}\times\bm{B}\) flow [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}/L_\mathrm{ref}\) ] (real*8).

  • \(\frac{T_\mathrm{s}\Gamma_{\mathrm{sM}}}{L_{p\mathrm{s}}}\) : Particle flux term by magnetic flutter [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}/L_\mathrm{ref}\) ] (real*8).

  • \(\frac{\Theta_{\mathrm{sE}}}{L_{T\mathrm{s}}}\) : Heat flux term by \(\bm{E}\times\bm{B}\) flow [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}/L_\mathrm{ref}\) ] (real*8).

  • \(\frac{\Theta_{\mathrm{sM}}}{L_{T\mathrm{s}}}\) : Heat flux term by magnetic flutter [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}/L_\mathrm{ref}\) ] (real*8).

The 0th and 1st components of \(S_\mathrm{s}\) \(D_\mathrm{s}\) correspond to non-zonal ( \(k_y \neq 0\) ) and zonal ( \(k_y = 0\) ) fluctuations, respectively.

Table 10 hst/gkvp_f0.48.ges.(ranks in 1 digits).(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

dtout_eng

MPI ranks

rank == 0

Total files

nprocs*(Total run numbers)

GKV unit

oges

Stored data

time, \(\Gamma_{\mathrm{sE}}\) , \(\Gamma_{\mathrm{sE}k_y}\)

where,

  • time : Simulation time \(t\) [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

  • \(\Gamma_{\mathrm{sE}}\) : Total particle flux by \(\bm{E} \times \bm{B}\) flow [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}v_\mathrm{ref}\) ] (real).

  • \(\Gamma_{\mathrm{sE}k_y}\) (0:global_ny): \(k_y\) spectrum of the particle flux by \(\bm{E} \times \bm{B}\) flow [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}v_\mathrm{ref}\) ] (real).

Table 11 hst/gkvp.gem.(ranks in 1 digits).(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

dtout_eng

MPI ranks

rank == 0

Total files

nprocs*(Total run numbers)

GKV unit

ogem

Stored data

time, \(\Gamma_{\mathrm{sM}}\) , \(\Gamma_{\mathrm{sM}k_y}\)

where,

  • time : Simulation time \(t\) [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

  • \(\Gamma_{\mathrm{sM}}\) : Total particle flux by magnetic flutter [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}v_\mathrm{ref}\) ] (real).

  • \(\Gamma_{\mathrm{sM}k_y}\) (0:global_ny): \(k_y\) spectrum of the particle flux by magnetic flutter [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}v_\mathrm{ref}\) ] (real).

Table 12 hst/gkvp.qes.(ranks in 1 digits).(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

dtout_eng

MPI ranks

rank == 0

Total files

nprocs*(Total run numbers)

GKV unit

oqes

Stored data

time, \(Q_{\mathrm{sE}}\) , \(Q_{\mathrm{sE}k_y}\)

where,

  • time : Simulation time \(t\) [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

  • \(Q_{\mathrm{sE}}\) : Total energy flux by \(\bm{E} \times \bm{B}\) flow [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}\) ] (real).

  • \(Q_{\mathrm{sE}k_y}\) (0:global_ny): \(k_y\) spectrum of the energy flux by \(\bm{E} \times \bm{B}\) flow [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}\) ] (real).

Table 13 hst/gkvp.qem.(ranks in 1 digits).(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

dtout_eng

MPI ranks

rank == 0

Total files

nprocs*(Total run numbers)

GKV unit

oqem

Stored data

time, \(Q_{\mathrm{sM}}\) , \(Q_{\mathrm{sM}k_y}\)

where,

  • time : Simulation time \(t\) [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

  • \(Q_{\mathrm{sM}}\) : Total energy flux by magnetic flutter [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}\) ] (real).

  • \(Q_{\mathrm{sM}k_y}\) (0:global_ny): \(k_y\) spectrum of the energy flux by magnetic flutter [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}v_\mathrm{ref}\) ] (real).

Table 14 hst/gkvp.wes.(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

dtout_eng

MPI ranks

rankg == 0

Total files

(Total run numbers)

GKV unit

owes

Stored data

time, \(W_{\mathrm{E}}\) , \(W_{\mathrm{E}k_y}\)

where,

  • time : Simulation time \(t\) [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

  • \(W_{\mathrm{E}}\) : Total electrostatic field energy including polarization [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}\) ] (real).

  • \(W_{\mathrm{E}k_y}\) (0:global_ny): \(k_y\) spectrum of the electrostatic field energy [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}\) ] (real).

Table 15 hst/gkvp.wem.(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

dtout_eng

MPI ranks

rankg == 0

Total files

(Total run numbers)

GKV unit

owem

Stored data

time, \(W_{\mathrm{M}}\) , \(W_{\mathrm{M}k_y}\)

where,

  • time : Simulation time \(t\) [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

  • \(W_{\mathrm{M}}\) : Total magnetic field energy [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}\) ] (real).

  • \(W_{\mathrm{M}k_y}\) (0:global_ny): \(k_y\) spectrum of the magnetic field energy [ \(\delta_\mathrm{ref}^2n_\mathrm{ref}T_\mathrm{ref}\) ] (real).

Table 16 hst/gkvp.eng.(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

dtout_eng

MPI ranks

rankg == 0

Total files

(Total run numbers)

GKV unit

oeng

Stored data

time, \(\sum_{k_x,k_y} \langle |\tilde{\phi}_{\bm{k}}|^2 \rangle\) , \(\sum_{k_x} \langle |\tilde{\phi}_{\bm{k}}|^2 \rangle\)

where,

  • time : Simulation time \(t\) [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

  • \(\sum_{k_x,k_y} \langle |\tilde{\phi}_{\bm{k}}|^2 \rangle\) : Squared amplitude of the perturbed electrostatic potential [ \((\delta_\mathrm{ref}T_\mathrm{ref}/e_\mathrm{ref})^2\) ] (real).

  • \(\sum_{k_x} \langle |\tilde{\phi}_{\bm{k}}|^2 \rangle\) (0:global_ny): \(k_y\) spectrum of the squared amplitude of the perturbed electrostatic potential [ \((\delta_\mathrm{ref}T_\mathrm{ref}/e_\mathrm{ref})^2\) ] (real).

Table 17 hst/gkvp.men.(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

dtout_eng

MPI ranks

rankg == 0

Total files

(Total run numbers)

GKV unit

omen

Stored data

time, \(\sum_{k_x,k_y} \langle |\tilde{A}_{\parallel\bm{k}}|^2 \rangle\) , \(\sum_{k_x} \langle |\tilde{A}_{\parallel\bm{k}}|^2 \rangle\)

where,

  • time : Simulation time \(t\) [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

  • \(\sum_{k_x,k_y} \langle |\tilde{A}_{\parallel\bm{k}}|^2 \rangle\) : Squared amplitude of the perturbed electrostatic potential [ \((\delta_\mathrm{ref}\rho_\mathrm{ref}B_\mathrm{ref})^2\) ] (real).

  • \(\sum_{k_x} \langle |\tilde{A}_{\parallel\bm{k}}|^2 \rangle\) (0:global_ny): \(k_y\) spectrum of the squared amplitude of the perturbed electrostatic potential [ \((\delta_\mathrm{ref}\rho_\mathrm{ref}B_\mathrm{ref})^2\) ] (real).

Table 18 hst/gkvp.dtc.(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

dtout_dtc

MPI ranks

rankg == 0

Total files

(Total run numbers)

GKV unit

odtc

Stored data

time, dt, dt_limit, dt_nl

where,

  • time : Simulation time \(t\) [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

  • dt : Time step size [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

  • dt_limit : Estimation of time step size limit [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

  • dt_nl : Estimation of time step size limit from nonlinear advection [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

Table 19 hst/gkvp.mtr.(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

Beginning of the run

MPI ranks

rankg == 0

Total files

(Total run numbers)

GKV unit

omtr

Stored data

time, \(\theta\) (or \(\varphi\) ), \(B, \frac{\partial B}{\partial x}, \frac{\partial B}{\partial y}, \frac{\partial B}{\partial z}, g^{xx}, g^{xy}, g^{xz}, g^{yy}, g^{yz}, g^{zz}, \sqrt{g}\)

where,

  • time : Simulation time \(t\) [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

  • \(\theta\) : Poloidal angle (or Toroidal angle \(\varphi\) when equib_type = “vmec”) (real)

  • \(B\) : Magnetic field strength [ \(B_\mathrm{ref}\) ] (real)

  • \(\frac{\partial B}{\partial x}, \frac{\partial B}{\partial y}, \frac{\partial B}{\partial z}\) : Derivative of \(B\) [ \(B_\mathrm{ref}/L_\mathrm{ref}\) ] (real)

  • \(g^{xx}, g^{xy}, g^{xz} [L_\mathrm{ref}^{-1}], g^{yy}, g^{yz} [L_\mathrm{ref}^{-1}], g^{zz} [L_\mathrm{ref}^{-2}]\) : Metric tensor (real)

  • \(\sqrt{g}\) : Jacobian [ \(L_\mathrm{ref}\) ] (real)

Table 20 hst/gkvp.frq.(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

dtout_eng (when calc_type == "linear" or "lin_freq" )

MPI ranks

rankg == 0

Total files

(Total run numbers)

GKV unit

ofrq

Stored data

time, omega

where,

  • time : Simulation time \(t\) [ \(L_\mathrm{ref}/v_\mathrm{ref}\) ] (real)

  • omega(1:global_ny) : \(k_y\) spectrum of complex linear frequency \(\omega\) = (real frequency, growthrate) [ \(v_\mathrm{ref}/L_\mathrm{ref}\) ] (real, real)

Complex frequency for \(k_x = 0\) at each time is evaluated as \(\omega = \omega_\mathrm{r} + i \gamma = \frac{\ln \tilde{\phi}_{\bm{k}}(t+\Delta t) - \ln \tilde{\phi}_{\bm{k}}(t)}{-i \Delta t}\) by assuming \(\tilde{\phi}_{\bm{k}}(t) \propto e^{-i\omega t}\) .

Table 21 hst/gkvp.dsp.(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

End of the run (when calc_type == "linear" or "lin_freq" )

MPI ranks

rankg == 0

Total files

(Total run numbers)

GKV unit

odsp

Stored data

ky, omega, diff, 1-ineq

where,

  • ky : Field-line-label (poloidal) wavenumber \(k_y\) [ \(\rho_\mathrm{ref}^{-1}\) ] (real)

  • omega : Complex linear frequency \(\omega\) = (real frequency, growthrate) [ \(v_\mathrm{ref}/L_\mathrm{ref}\) ] (real, real)

  • diff : Relative residual error \(\frac{\omega(t) - \omega(t-\Delta t)}{\omega(t)}\) (real, real)

  • 1-ineq : Convergence check based on Schwartz inequality (real)

At the end of run, estimated complex frequency for \(k_x = 0\) are dumped. If some modes are not yet converged, they are commented out.

Table 22 log/gkvp.(rankg in 6 digits).(ranks in 1 digit).log.(inum in 3 digits)

Field

Description

File type

Ascii

Output timing

As needed

MPI ranks

All

Total files

nprocw*nprocz*nprocv*nprocm*nprocs*(Total run numbers)

GKV unit

olog

Stored data

Simulation log