NEO Output Files

NEO output files are produced only if SILENT_FLAG = 0.

All NEO runtime information is written to out.neo.run.

Standard output files

Filename	Short description
out.neo.equil	Equilibrium/geometry input data
out.neo.f	First-order distribution function
out.neo.grid	Numerical grid parameters
out.neo.phi	Poloidal variation of first-order es potential
out.neo.theory	Neoclassical transport coefficients from analytic theory
out.neo.species	Mass and charge of all species
out.neo.theory_nclass	Neoclassical transport coefficients from the NCLASS code
out.neo.transport	Neoclassical transport coefficients from DKE solve
out.neo.transport_flux	Neoclassical fluxes in GB units from DKE solve
out.neo.transport_gv	Neoclassical fluxes from gyroviscosity
out.neo.vel	Poloidal variation of first-order flows
out.neo.vel_fourier	Poloidal variation of first-order flows (Fourier components)

Experimental profiles output files

Produced only if PROFILE_MODEL = 2.

Filename	Short description
out.neo.transport_exp	Neoclassical transport coefficients from DKE solve (in units)
out.neo.exp_norm	Normalizing experimental parameters (in units)

Rotation output files

Produced only if ROTATION_MODEL = 2.

Filename	Short description
out.neo.rotation	Strong rotation poloidal asymmetry parameters

Subroutine output

When neo is run in subroutine mode, the outputs are contained in a monolithic file named neo_interface. The NEO subroutine output parameters are as follows:

Parameter name	Short description	Normalization
neo_pflux_dke_out(1:11)	DKE solve particle flux	${\mathrm{\Gamma }}_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}})$
neo_efluxtot_dke_out(1:11)	DKE solve energy flux	$Q_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$
neo_efluxncv_dke_out(1:11)	DKE solve non-convective energy flux	$(Q_{\sigma }-{\omega }_0{\mathrm{\Pi }}_{\sigma })/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$
neo_mflux_dke_out(1:11)	DKE solve momentum flux	${\mathrm{\Pi }}_{\sigma }/(n_{\mathrm{norm}}{T}_{\mathrm{norm}}{a}_{\mathrm{norm}})$
neo_vpol_dke_out(1:11)	DKE solve poloidal flow	${\mathrm{v}}_{\theta ,\sigma }(\theta =0)/{\mathrm{v}}_{\mathrm{norm}}$
neo_vtor_dke_out(1:11)	DKE solve toroidal flow	${\mathrm{v}}_{\phi ,\sigma }(\theta =0)/{\mathrm{v}}_{\mathrm{norm}}$
neo_jpar_dke_out	DKE solve bootstrap current (parallel)	$⟨j_{\Vert }B⟩/(e{n}_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$
neo_jtor_dke_out	DKE solve bootstrap current (toroidal)	$⟨j_{\phi }/R⟩/⟨1/R⟩/(e{n}_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}})$
neo_pflux_gv_out(1:11)	Gyroviscosity particle flux	${\mathrm{\Gamma }}_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}})$
neo_efluxtot_gv_out(1:11)	Gyroviscosity energy flux	$Q_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$
neo_efluxncv_gv_out(1:11)	Gyroviscosity non-convective energy flux	$(Q_{\sigma }-{\omega }_0{\mathrm{\Pi }}_{\sigma })/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$
neo_mflux_gv_out(1:11)	Gyroviscosity momentum flux	${\mathrm{\Pi }}_{\sigma }/(n_{\mathrm{norm}}{T}_{\mathrm{norm}}{a}_{\mathrm{norm}})$
neo_pflux_thHH_out	Hinton-Hazeltine ion particle flux	${\mathrm{\Gamma }}_i/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}})$
neo_eflux_thHHi_out	Hinton-Hazeltine ion energy flux	$Q_i/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$
neo_eflux_thHHe_out	Hinton-Hazeltine electron energy flux	$Q_e/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$
neo_eflux_thCHi_out	Chang-Hinton ion energy flux	$Q_i/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$
neo_pflux_thHS_out(1:11)	Hirshman-Sigmar particle flux	${\mathrm{\Gamma }}_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}})$
neo_eflux_thS_out(1:11)	Hirshman-Sigmar energy flux	$Q_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$
neo_jpar_thS_out	Sauter bootstrap current (parallel)	$⟨j_{\Vert }B⟩/(e{n}_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$
neo_jtor_thS_out	Sauter bootstrap current (toroidal)	$⟨j_{\phi }/R⟩/⟨1/R⟩/(e{n}_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}})$
neo_pflux_nclass_out(1:11)	NCLASS solve particle flux	${\mathrm{\Gamma }}_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}})$
neo_efluxtot_nclass_out(1:11)	NCLASS solve energy flux	$Q_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$
neo_vpol_nclass_out(1:11)	NCLASS solve poloidal flow	${\mathrm{v}}_{\theta ,\sigma }(\theta =0)/{\mathrm{v}}_{\mathrm{norm}}$
neo_vtor_nclass_out(1:11)	NCLASS solve toroidal flow	${\mathrm{v}}_{\phi ,\sigma }(\theta =0)/{\mathrm{v}}_{\mathrm{norm}}$
neo_jpar_nclass_out	NCLASS solve bootstrap current (parallel)	$⟨j_{\Vert }B⟩/(e{n}_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$

---

Detailed description of NEO output files

out.neo.equil

Description

Equilibrium/geometry input data

Format

Rectangular array of ASCII data:

• rows: $\mathtt{N}\mathtt{\_}\mathtt{RADIAL}$

• cols: $7+5\times \mathtt{N}\mathtt{\_}\mathtt{SPECIES}$

1. $r/a$: normalized midplane minor radius

2. $(\partial {\mathrm{\Phi }}_0/\partial r)(ae/T_{\mathrm{norm}})$: normalized equilibrium-scale radial electric field

3. $q$: safety factor

4. ${\rho }_{\ast }=(c\sqrt{m_{\mathrm{norm}}{T}_{\mathrm{norm}}})/(e{B}_{unit}a)$: ratio of Larmor radius of normalizing species to the normalizing length

5. $R_0/a$: normalized flux-surface-center major radius

6. ${\omega }_0(a/{\mathrm{v}}_{\mathrm{norm}})$: normalized toroidal angular frequency

7. $(d{\omega }_0/dr)(a^2/{\mathrm{v}}_{\mathrm{norm}})$: normalized toroidal rotation shear

For each species $\sigma$:

8. $n_{\sigma }/n_{\mathrm{norm}}$: normalized equilibrium-scale density

9. $T_{\sigma }/T_{\mathrm{norm}}$: normalized equilibrium-scale temperature

10. $a/L_{n\sigma }=-a(d\mathrm{l}\mathrm{n}{n}_{\sigma }/dr)$: normalized equilibrium-scale density gradient scale length

11. $a/L_{T\sigma }=-a(d\mathrm{l}\mathrm{n}{T}_{\sigma }/dr)$: normalized equilibrium-scale temperature gradient scale length

12. ${\tau }_{\sigma \sigma }^{-1}(a/{\mathrm{v}}_{\mathrm{norm}})$: normalized collision frequency

---

out.neo.exp_norm

Description

Normalizing experimental parameters (in units)

Format

Rectangular array of ASCII data:

• rows: $\mathtt{N}\mathtt{\_}\mathtt{RADIAL}$

• cols: $7$

1. $r/a$: normalized midplane minor radius

2. $a$: normalizing length (m)

3. $m_{\mathrm{norm}}$: normalizing mass (e-27 kg)

4. $n_{\mathrm{norm}}$: normalizing equilibrium-scale density (e19/m^3)

5. $T_{\mathrm{norm}}$: normalizing equilibrium-scale temperature (keV)

6. ${\mathrm{v}}_{\mathrm{norm}}$: normalizing thermal speed (m/s)

7. $B_{unit}$: normalizing magnetic field (T)

---

out.neo.f

Description

First-order distribution function solution (dimensionless), specifically vector of ${\hat{g}}_{a,ie,ix,it}$ (first-order non-adiabatic distribution function for each species $a$), where

$$g_a(r,\theta ,x_a,\xi )=f_{0a}(r,\theta ,x_a)\sum _{ie=0}^{\mathtt{N}\mathtt{\_}\mathtt{ENERGY}}\sum _{ix=0}^{\mathtt{N}\mathtt{\_}\mathtt{XI}}{L}_{ie}^{k(ix)+1/2}(x_a^2)P_{ix}(\xi ){\hat{g}}_{a,ie,ix,it}(\theta )$$

where $f_{0a}$ is the zeroth-order distribution function (Maxwellian), $L_{ie}$ are associated Laguerre polynomials and $P_{ix}$ are Legendre polynomials, $k(ix)=0$ for ix=0 and $k(ix)=1$ for ix>0, $\xi =\mathrm{v}/{\mathrm{v}}_{\Vert }$ is the cosine of the pitch angle, and $x_a=\mathrm{v}/\sqrt{2{\mathrm{v}}_{ta}}$ is the normalized energy.

Format

Vector of ASCII data:

• $(\mathtt{N}\mathtt{\_}\mathtt{RADIAL})\times (\mathtt{N}\mathtt{\_}\mathtt{SPECIES})\times (\mathtt{N}\mathtt{\_}\mathtt{ENERGY}+1)\times (\mathtt{N}\mathtt{\_}\mathtt{XI}+1)\times (\mathtt{N}\mathtt{\_}\mathtt{THETA}$)

---

out.neo.grid

Description

Numerical grid parameters

Format

Vector of ASCII data:

• $5+\mathtt{N}\mathtt{\_}\mathtt{THETA}+\mathtt{N}\mathtt{\_}\mathtt{RADIAL}$

1. $\mathtt{N}\mathtt{\_}\mathtt{SPECIES}$: number of kinetic species

2. $\mathtt{N}\mathtt{\_}\mathtt{ENERGY}$: number of energy polynomials

3. $\mathtt{N}\mathtt{\_}\mathtt{XI}$: number of $\xi =\mathrm{v}/{\mathrm{v}}_{\Vert }$ (cosine of pitch angle) polynomials

4. $\mathtt{N}\mathtt{\_}\mathtt{THETA}$: number of theta gridpoints

5. ${\theta }_j$: theta gridpoints (j=1..N_THETA)

6. $\mathtt{N}\mathtt{\_}\mathtt{RADIAL}$: number of radial gridpoints

7. $r_j/a$: normalized radial gridpoints (j=1..N_RADIAL)

---

out.neo.phi

Description

Neoclassical first-order electrostatic potential (normalized) vs. $\theta$

Format

Rectangular array of ASCII data:

• rows: $\mathtt{N}\mathtt{\_}\mathtt{RADIAL}$

• cols: $\mathtt{N}\mathtt{\_}\mathtt{THETA}$

1. $\frac{e{\mathrm{\Phi }}_1({\theta }_j)}{T_{\mathrm{norm}}}$: first-order electrostatic potential vs. ${\theta }_j$ (j=1…N_THETA)

---

out.neo.rotation

Description

Strong rotation poloidal asymmetry parameters (normalized)

Define:

• ${\mathrm{\Phi }}_{\ast }={\mathrm{\Phi }}_0-{\mathrm{\Phi }}_0(\theta =0)$

• ${\epsilon }_{\sigma }=\frac{z_{\sigma }e}{T_{\sigma }}-\frac{m_{\sigma }{\omega }_0^2}{2T_{\sigma }}[R^2-R^2(\theta =0)]$

• $e_{0\sigma }=⟨e^{-{\epsilon }_{\sigma }}⟩$

• $e_{1\sigma }=⟨e^{-{\epsilon }_{\sigma }}\frac{z_{\sigma }e{\mathrm{\Phi }}_{\ast }}{T_{\sigma }}⟩$

• $e_{2\sigma }=a_{\mathrm{norm}}⟨e^{-{\epsilon }_{\sigma }}\frac{z_{\sigma }e}{T_{\sigma }}\frac{\partial {\mathrm{\Phi }}_{\ast }}{\partial r}⟩$

• $e_{3\sigma }=\frac{1}{a_{\mathrm{norm}}^2}⟨e^{-{\epsilon }_{\sigma }}[R^2-R^2(\theta =0)]⟩$

• $e_{4\sigma }=\frac{1}{a_{\mathrm{norm}}}⟨e^{-{\epsilon }_{\sigma }}\frac{\partial [R^2-R^2(\theta =0)]}{\partial r}⟩$

• $e_{5\sigma }=a_{\mathrm{norm}}⟨e^{-{\epsilon }_{\sigma }}\frac{\partial \ln \sqrt{g}}{\partial r}⟩-a_{\mathrm{norm}}⟨e^{-{\epsilon }_{\sigma }}⟩⟨\frac{\partial \ln \sqrt{g}}{\partial r}⟩$

• For anisotropic species, all temperatures are interpreted as $T_{\Vert }$, the total energy is modified by ${\epsilon }_{\sigma }\to {\epsilon }_{\sigma }+{\lambda }_{\mathrm{a}\mathrm{n}\mathrm{i}\mathrm{s}\mathrm{o},\sigma }(r,\theta )$, and we define the additional term $e_{6\sigma }=-a_{\mathrm{norm}}⟨e^{-{\epsilon }_{\sigma }}\frac{\partial {\lambda }_{\mathrm{a}\mathrm{n}\mathrm{i}\mathrm{s}\mathrm{o},\sigma }}{\partial r}⟩$

• $F_{V\sigma }=\frac{1}{e_{0\sigma }}[-e_{2\sigma }+e_{3\sigma }{a}_{\mathrm{norm}}^3\frac{{\omega }_0}{{\mathrm{v}}_{t\sigma }}\frac{d{\omega }_0}{dr}+e_{4\sigma }{a}_{\mathrm{norm}}^2\frac{{\omega }_0^2}{2{\mathrm{v}}_{t\sigma }^2}+e_{1\sigma }{a}_{\mathrm{norm}}\frac{d\ln T_{\sigma }}{dr}-e_{3\sigma }{a}_{\mathrm{norm}}^3\frac{d\ln T_{\sigma }}{dr}\frac{{\omega }_0^2}{2{\mathrm{v}}_{t\sigma }^2}+e_{5\sigma }+e_{6\sigma }]$

Format

Rectangular array of ASCII data:

• rows: $\mathtt{N}\mathtt{\_}\mathtt{RADIAL}$

• cols: $2+2\times \mathtt{N}\mathtt{\_}\mathtt{SPECIES}+\mathtt{N}\mathtt{\_}\mathtt{THETA}+2\times \mathtt{N}\mathtt{\_}\mathtt{SPECIES}\times \mathtt{N}\mathtt{\_}\mathtt{THETA}$

Fixed entries:

1. $r/a$: normalized midplane minor radius

2. $\frac{e⟨{\mathrm{\Phi }}_{\ast }⟩}{T_{\mathrm{norm}}}$: difference between the flux-surface-averaged equilibrium-scale potential and the value at the outboard midplane (0 in the diamagnetic ordering limit)

For each species $\sigma$:

3. $\frac{1}{e_{0\sigma }}=\frac{n_{\sigma }}{⟨n_{\sigma }⟩}$: ratio of the density at the outboard midplane to the flux-surface-averaged equilibrium-scale density (1 in the diamagnetic ordering limit)

4. $F_{V\sigma }$: Factor related to the transformation of the particle flux convection (presently only valid in $s-\alpha$ geometry)

For each ${\theta }_j$, j=1..N_THETA

5. $\frac{e{\mathrm{\Phi }}_{\ast }({\theta }_j)}{T_{\mathrm{norm}}}$: difference between the equilibrium-scale potential and the value at the outboard midplane (0 in the diamagnetic ordering limit)

6. $\frac{n_{\sigma }({\theta }_j)}{n_{\sigma }(\theta =0)}$: poloidal variation of the equilibrium-scale density normalized to the value at the outboard midplane (1 in the diamagnetic ordering limit)

---

out.neo.species

Description

Mass and charge of all species

Format

Rectangular array of ASCII data:

• cols: $2\times \mathtt{N}\mathtt{\_}\mathtt{SPECIES}$

1. For each species $\sigma$:

   • $m_{\sigma }/m_{\mathrm{norm}}$: species mass (we suggest always taking deuterium as the normalizing mass)

   • $z_{\sigma }$: species charge

---

out.neo.theory

Description

Neoclassical transport coefficients from analytic theory (normalized)

• Only the Hirshman-Sigmar quantities are meaningful for multiple-ion species plasmas.

• None of the theories are valid with strong rotation effects included.

Theory references

• Hinton-Hazltine flows and fluxes: Rev. Mod. Phys., vol. 48, 239 (1976)

• Chang-Hinton ion heat flux: Phys. Plasmas, vol. 25, 1493 (1982)

• Taguchi ion heat flux (modified with Chang-Hinton collisional interpolation factor): PPCF, vol. 30, 1897 (1988)

• Sauter et al. bootstrap current model: Phys. Plasmas, vol. 6, 2834 (1999)

• Hinton-Rosenbluth potential: Phys. Fluids 16, 836 (1973)

• Hirshman-Sigmar fluxes: Phys. Fluids, vol. 20, 418 (1977)

• Koh et al. bootstrap current model: Phys. Plasmas, vol. 19, 072505 (2012)

Format

Rectangular array of ASCII data:

• rows: $\mathtt{N}\mathtt{\_}\mathtt{RADIAL}$

• cols: $17+2\times \mathtt{N}\mathtt{\_}\mathtt{SPECIES}$

1. $r/a$: normalized midplane minor radius

2. HH ${\mathrm{\Gamma }}_i/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}})$: Hinton-Hazeltine second-order radial particle flux (ambipolar)

3. HH $Q_i/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$: Hinton-Hazeltine second-order radial energy flux (ion)

4. HH $Q_e/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$: Hinton-Hazeltine second-order radial energy flux (electron)

5. HH $⟨j_{\Vert }B⟩/(e{n}_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$: Hinton-Hazeltine first-order bootstrap current

6. HH $k_i$: Hinton-Hazeltine first-order dimensionless flow coefficient (ion)

7. HH $⟨u_{\Vert ,i}B⟩/({\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$: Hinton-Hazeltine first-order parallel flow (ion)

8. HH ${\mathrm{v}}_{\theta ,i}(\theta =0)/{\mathrm{v}}_{\mathrm{norm}}$: Hinton-Hazeltine first-order poloidal flow at the outboard midplane (ion)

9. CH $Q_i/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$: Chang-Hinton second-order radial energy flux (ion)

10. TG $Q_i/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$: Taguchi second-order radial energy flux (ion)

11. S $⟨j_{\Vert }B⟩/(e{n}_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$: Sauter first-order bootstrap current

12. S $k_i$: Sauter first-order dimensionless flow coefficient (ion)

13. S $⟨u_{\Vert ,i}B⟩/({\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$: Sauter first-order parallel flow (ion)

14. S ${\mathrm{v}}_{\theta ,i}(\theta =0)/{\mathrm{v}}_{\mathrm{norm}}$: Sauter first-order poloidal flow at the outboard midplane (ion)

15. HR $⟨(e{\mathrm{\Phi }}_1/T_{\mathrm{norm}}{)}^2⟩$: Hinton-Rosenbluth first-order electrostatic potential

16. For each species $\sigma$:

    • HS ${\mathrm{\Gamma }}_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}})$: Hirshman-Sigmar second-order radial particle flux

    • HS $Q_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$: Hirshman-Sigmar second-order radial energy flux

18. K $⟨j_{\Vert }B⟩/(e{n}_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$: Koh first-order bootstrap current

19. S $⟨j_{\Vert }B⟩/(e{n}_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$: Sauter first-order bootstrap current

---

out.neo.theory_nclass

Description

Neoclassical transport coefficients from the NCLASS code (normalized)

• Only produced if SIM_MODEL = 1 or 3.

• Note that for local mode (PROFILE_MODEL = 1), it is assumed in the NCLASS calculation that the normalizing mass is the mass of deuterium and that the input collision frequencies are self-consistent across all species.

• NCLASS reference: W.A. Houlberg, et al, Phys. Plasmas, vol. 4, 3230 (1997)

Format

Rectangular array of ASCII data:

• rows: $\mathtt{N}\mathtt{\_}\mathtt{RADIAL}$

• cols: $2+5\times \mathtt{N}\mathtt{\_}\mathtt{SPECIES}$

1. $r/a$: normalized midplane minor radius

2. $⟨j_{\Vert }B⟩/(e{n}_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$: first-order bootstrap current

For each species $\sigma$:

3. ${\mathrm{\Gamma }}_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}})$: second-order radial particle flux

4. $Q_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$: second-order radial energy flux

5. $⟨u_{\Vert ,\sigma }B⟩/({\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$: first-order parallel flow

6. ${\mathrm{v}}_{\theta ,\sigma }(\theta =0)/{\mathrm{v}}_{\mathrm{norm}}$: first-order poloidal flow at the outboard midplane

7. ${\mathrm{v}}_{\phi ,\sigma }(\theta =0)/{\mathrm{v}}_{\mathrm{norm}}$: first-order toroidal flow at the outboard midplane

---

out.neo.transport

Description

Neoclassical transport coefficients from DKE solve (normalized)

Format

Rectangular array of ASCII data:

• rows: $\mathtt{N}\mathtt{\_}\mathtt{RADIAL}$

• cols: $5+8\times \mathtt{N}\mathtt{\_}\mathtt{SPECIES}$

1. $r/a$: normalized midplane minor radius

2. $⟨(e{\mathrm{\Phi }}_1/T_{\mathrm{norm}}{)}^2⟩$: first-order electrostatic potential

3. $⟨j_{\Vert }B⟩/(e{n}_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$: first-order bootstrap current

4. $v_{\phi }^{(0)}(\theta =0)/{\mathrm{v}}_{\mathrm{norm}}$: zeroth-order toroidal flow at the outboard midplane ($v_{\phi }^{(0)}={\omega }_{0}R$)

5. $⟨u_{\Vert }^{(0)}B⟩/({\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$: zeroth-order parallel flow ($u_{\Vert }^{(0)}={\omega }_{0}I/B$)

For each species $\sigma$:

6. ${\mathrm{\Gamma }}_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}})$: second-order radial particle flux

7. $Q_{\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}}{T}_{\mathrm{norm}})$: second-order radial energy flux

8. ${\mathrm{\Pi }}_{\sigma }/(n_{\mathrm{norm}}{T}_{\mathrm{norm}}{a}_{\mathrm{norm}})$: second-order radial momentum flux

9. $⟨u_{\Vert ,\sigma }B⟩/({\mathrm{v}}_{\mathrm{norm}}{B}_{unit})$: first-order parallel flow

10. $k_{\sigma }$: first-order dimensionless flow coefficient

11. $K_{\sigma }/(n_{\mathrm{norm}}{rmv}_{\mathrm{norm}}/B_{unit})$: first-order dimensional flow coefficient

12. ${\mathrm{v}}_{\theta ,\sigma }(\theta =0)/{\mathrm{v}}_{\mathrm{norm}}$: first-order poloidal flow at the outboard midplane

13. ${\mathrm{v}}_{\phi ,\sigma }(\theta =0)/{\mathrm{v}}_{\mathrm{norm}}$: first-order toroidal flow at the outboard midplane

---

out.neo.transport_exp

Description

Neoclassical transport coefficients from DKE solve (in units)

Format

Rectangular array of ASCII data:

• rows: $\mathtt{N}\mathtt{\_}\mathtt{RADIAL}$

• cols: $5+8\times \mathtt{N}\mathtt{\_}\mathtt{SPECIES}$

1. $r$: midplane minor radius ($m$)

2. $⟨({\mathrm{\Phi }}_1{)}^2⟩$: first-order electrostatic potential ($V^2$)

3. $⟨j_{\Vert }B⟩/B_{unit}$: first-order bootstrap current ($A/m^2$)

4. $v_{\phi }^{(0)}(\theta =0)$: zeroth-order toroidal flow at the outboard midplane ($v_{\phi }^{(0)}={\omega }_{0}R$) ($m/s$)

5. $⟨u_{\Vert }^{(0)}B⟩/B_{unit}$: zeroth-order parallel flow ($u_{\Vert }^{(0)}={\omega }_{0}I/B$) ($m/s$)

For each species $\sigma$:

6. ${\mathrm{\Gamma }}_{\sigma }$: second-order radial particle flux ($e19m^{-2}{s}^{-1}$)

7. $Q_{\sigma }$: second-order radial energy flux ($W/m^2$)

8. ${\mathrm{\Pi }}_{\sigma }$: second-order radial momentum flux ($N/m$)

9. $⟨u_{\Vert ,\sigma }B⟩/B_{unit}$: first-order parallel flow ($m/s$)

10. $k_{\sigma }$: first-order dimensionless flow coefficient

11. $K_{\sigma }$: first-order dimensional flow coefficient ($e19/(m^{2}sT)$)

12. ${\mathrm{v}}_{\theta ,\sigma }(\theta =0)$: first-order poloidal flow at the outboard midplane ($m/s$)

13. ${\mathrm{v}}_{\phi ,\sigma }(\theta =0)$: first-order toroidal flow at the outboard midplane ($m/s$)

---

out.neo.transport_flux

Description

Neoclassical fluxes in GB units (defined below) from DKE solve

$$\begin{array}{r}\begin{array}{rl}{\mathrm{\Gamma }}_{GB}& =n_e{c}_s{\left({\rho }_{s,\mathrm{u}\mathrm{n}\mathrm{i}\mathrm{t}}/a\right)}^2\\ Q_{GB}& =n_e{c}_s{T}_e{\left({\rho }_{s,\mathrm{u}\mathrm{n}\mathrm{i}\mathrm{t}}/a\right)}^2\\ {\mathrm{\Pi }}_{GB}& =n_e{T}_{e}a{\left({\rho }_{s,\mathrm{u}\mathrm{n}\mathrm{i}\mathrm{t}}/a\right)}^2\end{array}\end{array}$$

where $c_s=\sqrt{T_e/m_D}$ and ${\displaystyle {\rho }_{s,\mathrm{u}\mathrm{n}\mathrm{i}\mathrm{t}}=\frac{c_s}{e{B}_{\mathrm{u}\mathrm{n}\mathrm{i}\mathrm{t}}/(m_{D}c)}}$.

Format

Rectangular array of ASCII data:

• rows: $\mathtt{N}\mathtt{\_}\mathtt{RADIAL}\times 3\times \mathtt{N}\mathtt{\_}\mathtt{SPECIES}$

• cols: $3$

For each species $\sigma$:

1. row of DKE (${\mathrm{\Gamma }}_{\sigma }/{\mathrm{\Gamma }}_{GB}$, $Q_{\sigma }/Q_{GB}$, ${\mathrm{\Pi }}_{\sigma }/{\mathrm{\Pi }}_{GB}$): second-order radial particle, energy, and momentum fluxes from DKE solve

For each species $\sigma$:

2. row of GV (${\mathrm{\Gamma }}_{\sigma }/{\mathrm{\Gamma }}_{GB}$, $Q_{\sigma }/Q_{GB}$, ${\mathrm{\Pi }}_{\sigma }/{\mathrm{\Pi }}_{GB}$): second-order radial particle, energy, and momentum fluxes from gyroviscosity

For each species $\sigma$:

3. row of TGYRO (${\mathrm{\Gamma }}_{\sigma }/{\mathrm{\Gamma }}_{GB}$, $Q_{\sigma }/Q_{GB}$, ${\mathrm{\Pi }}_{\sigma }/{\mathrm{\Pi }}_{GB}$): : second-order radial particle, energy, and momentum fluxes for transport equations

---

out.neo.transport_gv

Description

Neoclassical fluxes from gyroviscosity (normalized)

• These fluxes are nonzero only for the case of combined sonic rotation with up-down asymmetric flux surfaces.

• In the transport equations, these fluxes should be added to the fluxes from the DKE solve.

• Reference: H. Sugama and W. Horton, Phys. Plasmas, vol. 4, 405 (1997).

Format

Rectangular array of ASCII data:

• rows: $\mathtt{N}\mathtt{\_}\mathtt{RADIAL}$

• cols: $1+3\times \mathtt{N}\mathtt{\_}\mathtt{SPECIES}$

1. $r/a$: normalized midplane minor radius

For each species $\sigma$:

2. ${\mathrm{\Gamma }}_{gv,\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}})$: Gyroviscous second-order radial particle flux

3. $Q_{gv,\sigma }/(n_{\mathrm{norm}}{\mathrm{v}}_{\mathrm{norm}})$: Gyroviscous second-order radial energy flux

4. ${\mathrm{\Pi }}_{gv,\sigma }/(n_{\mathrm{norm}}{T}_{\mathrm{norm}}{a}_{\mathrm{norm}})$: Gyroviscous second-order radial momentum flux

---

out.neo.vel

Description

Poloidal variation of first-order flows (normalized)

Format

Rectangular array of ASCII data:

• rows: $\mathtt{N}\mathtt{\_}\mathtt{RADIAL}$

• cols: $\mathtt{N}\mathtt{\_}\mathtt{SPECIES}\times \mathtt{N}\mathtt{\_}\mathtt{THETA}$

For each species $\sigma$:

1. $u_{\Vert ,\sigma }({\theta }_j)/{\mathrm{v}}_{\mathrm{norm}}$: first-order parallel flow vs. ${\theta }_j(j=1\dots \mathtt{N}\mathtt{\_}\mathtt{THETA})$

---

out.neo.vel_fourier

Description

Poloidal variation of first-order flows (normalized) in Fourier series components

$$u(\theta )=\sum _{j=0}^{\mathtt{N}\mathtt{\_}\mathtt{THETA}}{u}_{cj}\cos (j\theta )+u_{sj}\sin (j\theta )$$

Format

Rectangular array of ASCII data:

• rows: $\mathtt{N}\mathtt{\_}\mathtt{RADIAL}$

• cols: $\mathtt{N}\mathtt{\_}\mathtt{SPECIES}\times 6\times (\mathtt{M}\mathtt{\_}\mathtt{THETA}+1)$ where $\mathtt{M}\mathtt{\_}\mathtt{THETA}=\frac{\mathtt{N}\mathtt{\_}\mathtt{THETA}-1}{2}-1$

For each species $\sigma$:

1. For j=0..M_THETA, $u_{\Vert ,\sigma ,cj}$: cosine-component of first-order parallel flow

2. For j=0..M_THETA, $u_{\Vert ,\sigma ,sj}$: sine-component of first-order parallel flow

3. For j=0..M_THETA, $u_{\theta ,\sigma ,cj}$: cosine-component of first-order poloidal flow

4. For j=0..M_THETA, $u_{\theta ,\sigma ,sj}$: sine-component of first-order poloidal flow

5. For j=0..M_THETA, $u_{\phi ,\sigma ,cj}$: cosine-component of first-order toroidal flow

6. For j=0..M_THETA, $u_{\phi ,\sigma ,sj}$: sine-component of first-order toroidal flow