# Alphabetical list for input.neo

## ANISO_MODEL_*

**Definition**

Parameter which selects whether to treat a species with an anisotropic temperature model.

**Choices**

ANISO_MODEL_*=1: isotropic temperature model

ANISO_MODEL_*=2: anisotropic temperature model

This option is presently not available for experimental profiles (PROFILE_MODEL = 2).

This model requires ROTATION_MODEL = 2 due to the induced poloidal density asymmetry.

The parallel and perpendicular temperature TEMP_PARA_* and TEMP_PERP_* and the parallel and perpendicular temperature gradient scale lengths DLNTDR_PARA_* and DLNTDR_PERP_* must also be set.

The parameters TEMP_* and DLNTDR_* are not used. The effective Maxwellian temperature in the DKE is determined internally based on the parallel and perpendicular temperatures.

**Comments**

DEFAULT: 1

The anisotropic model of each species 1-11 is set as: ANISO_MODEL_1, ANISO_MODEL_2, ANISO_MODEL_3,…

The subroutine interface parameter is specified as a vector: neo_aniso_model_in(1:11)

## BETA_STAR

**Definition**

The normalized effective pressure gradient:

where \(B_{\rm unit}(r)=(q/r)\psi^\prime\) is the effective magnetic field strength and \(p=\sum_a n_a T_a\) is the total plasma pressure.

**Comments**

DEFAULT: 0.0

NOTE: This parameter is

*not*used in the standard DKE equation! It is only used in the case of an anisotropic temperature species (e.g. ANISO_MODEL_* = 2) to compute \(d\Phi_*/dr\).This is only active with EQUILIBRIUM_MODEL = 2 (the Miller equilibrium model).

When experimental profiles are used (PROFILE_MODEL = 2), \(\beta_*\) is computed internally from the profile parameters in input.profiles and the normalizing length scale is the plasma minor radius.

## BTCCW

**Definition**

Parameter which selects the orientation of the toroidal magnetic field \(B_t\) relative to the toroidal angle \(\varphi\).

**Choices**

BTCCW = 1: Counter-clockwise when viewed from above the torus - negative \(\hat{e}_{\varphi}\) for the right-handed coordinate system \((r,\theta,\varphi)\). Thus, \(B_t\) is oriented along the negative \(\hat{e}_{\varphi}\) direction.

BTCCW = -1: Clockwise when viewed from above the torus - positive \(\hat{e}_{\varphi}\) for the right-handed coordinate system \((r,\theta,\varphi)\). Thus, \(B_t\) is oriented along the positive \(\hat{e}_{\varphi}\) direction.

**Comments**

DEFAULT: -1

In DIII-D, typically BTCCW = 1.

When experimental profiles are used (PROFILE_MODEL = 2), the orientiation of BT is inferred from input.profiles.

## COLLISION_MODEL

**Definition**

Parameter which selects the collision operator model.

**Choices**

COLLISION_MODEL = 1: Connor model.

COLLISION_MODEL = 2: Zeroth-order Hirshman-Sigmar model.

COLLISION_MODEL = 3: Full Hirshman-Sigmar model.

COLLISION_MODEL = 4: Full linearized Fokker-Plank operator.

COLLISION_MODEL = 5: FP test particle operator with ad hoc field particle operator.

**Comments**

DEFAULT: 4

## DELTA

**Definition**

Average triangularity, \(\delta\), of the flux surface:

where \(\delta_{+}\) is the upper triangularity and \(\delta_{-}\) is the lower triangularity.

**Comments**

DEFAULT: 0.0

This is only active with EQUILIBRIUM_MODEL = 2 (the Miller equilibrium model).

When experimental profiles are used (PROFILE_MODEL = 2), the triangularity as a function of radius is read from input.profiles.

## DENS_*

**Definition**

The normalized equilibrium-scale density:

**Commments**

DEFAULT: DENS_1=1.0, DENS_2=DENS_3=…=0.0

The density of each species 1-11 is set as: DENS_1, DENS_2, DENS_3,…

When experimental profiles are used (PROFILE_MODEL = 2), the density as a function of radius is read from input.profiles and the normalizing density is the local density of Species 1, \(n_{\rm norm}(r)=n_{0,{\rm species 1}}\).

When rotation effects are included (ROTATION_MODEL = 2), this parameter is the value at the outboard midplane (\(\theta=0\)).

The subroutine interface parameter is specified as a vector: neo_dens_in(1:11)

## DLNNDR_*

**Definition**

The normalized equilibrium-scale density gradient scale length:

**Commments**

DEFAULT: 1.0

The density gradient scale length of each species 1-11 is set as: DLNNDR_1, DLNNDR_2, DLNNDR_3,…

When experimental profiles are used (PROFILE_MODEL = 2), the density as a function of radius is read from input.profiles and the density gradient is computed internally. The normalizing length is the plasma minor radius.

When rotation effects are included (ROTATION_MODEL = 2), this parameter is the value at the outboard midplane (\(\theta=0\)).

The subroutine interface parameter is specified as a vector: neo_dlnndr_in(1:11)

## DLNNDRE_ADE

**Definition**

The normalized equilibrium-scale density gradient scale length of the electrons for the case of adiabatic electrons:

**Commments**

DEFAULT: 1.0

This parameter does not enter the DKE and is used only as a diagnostic for strong rotation (ref:neo_rotation_model = 2), for which it is the value at the outboard midplane (\(\theta=0\)).

This paramter is used only if no species with Z < 0 is specified.

When experimental profiles are used (PROFILE_MODEL = 2), the density as a function of radius is read from input.profiles and the density gradient is computed internally. The normalizing length is the plasma minor radius.

## DLNTDR_*

**Definition**

The normalized equilibrium-scale temperature gradient scale length:

**Commments**

DEFAULT: 1.0

The temperature gradient scale length of each species 1-11 is set as: DLNTDR_1, DLNTDR_2, DLNTDR_3,…

When experimental profiles are used (PROFILE_MODEL = 2), the temperature as a function of radius is read from input.profiles and the temperature gradient is computed internally. The normalizing length is the plasma minor radius.

The subroutine interface parameter is specified as a vector: neo_dlntdr_in(1:11)

## DLNTDR_PARA_*

**Definition**

The normalized equilibrium-scale parallel temperature gradient scale length:

**Commments**

DEFAULT: 1.0

The parallel temperature gradient scale length of each species 1-11 is set as: DLNTDR_PARA_1, DLNTDR_PARA_2, DLNTDR_PARA_3,…

This parameter is used only when the species’ anisotropic flag is set (ANISO_MODEL_* = 2).

The subroutine interface parameter is specified as a vector: neo_dlntdr_para_in(1:11)

## DLNTDR_PERP_*

**Definition**

The normalized equilibrium-scale perpendicular temperature gradient scale length:

**Commments**

DEFAULT: 1.0

The perpendicular temperature gradient scale length of each species 1-11 is set as: DLNTDR_PERP_1, DLNTDR_PERP_2, DLNTDR_PERP_3,…

This parameter is used only when the species’ anisotropic flag is set (ANISO_MODEL_* = 2).

The subroutine interface parameter is specified as a vector: neo_dlntdr_perp_in(1:11)

## DLNTDRE_ADE

**Definition**

The normalized equilibrium-scale temperature gradient scale length of the electrons for the case of adiabatic electrons:

**Commments**

DEFAULT: 1.0

This parameter does not enter the DKE and is used only as a diagnostic for strong rotation (ref:neo_rotation_model = 2), for which it is the value at the outboard midplane (\(\theta=0\)).

This paramter is used only if no species with Z < 0 is specified.

When experimental profiles are used (PROFILE_MODEL = 2), the temperature as a function of radius is read from input.profiles and the temperature gradient is computed internally. The normalizing length is the plasma minor radius.

## DPHI0DR

**Definition**

The normalized equilibrium-scale radial electric field:

such that

**Comments**

DEFAULT: 0.0

When experimental profiles are used (PROFILE_MODEL = 2), this is computed internally from the profile parameters is in input.profiles and the normalizing length scale is the plasma minor radius. See also the parameter PROFILE_ERAD0_MODEL, which allows the simulation to be done with DPHI0DR = 0 regardless of the value in input.profiles.

If sonic rotation effects are included (ROTATION_MODEL = 2), then this parameter is ignored and \(E_r^{(0)}\) is assumed to be zero. With experimental profiles, this means that the \(E_r^{(0)}\) in input.profiles is assumed to be the lowest-order field in sonic rotation theory, i.e. \(E_r^{(-1)}\),and is used to compute the lowest-order sonic toroidal rotation parameters, OMEGA_ROT and OMEGA_ROT_DERIV.

## EPAR0

**Definition**

The normalized equilibrium-scale inductive electric field:

**Comments**

DEFAULT: 0.0

In the neo theory module, the input \(\left< E_\| B \right>\) is used directly.

For the DKE, it into the RHS neoclassical source term as

\[{\rm v_\|} \left< E_\| B \right> \frac{B}{\left< B^2 \right>}\]\(E_\|\) is not presently in input.profiles. When experimental profiles are used (PROFILE_MODEL = 2), EPAR0 is read from input.neo and is assumed to be radially constant.

For the Spitzer problem (SPITZER_MODEL = 1), use EPAR0_SPITZER instead.

## EPAR0_SPITZER

**Definition**

The normalized equilibrium-scale inductive electric field for use in the Spitzer problem:

**Comments**

DEFAULT: 1.0

For the DKE, we assume that \(E_\varphi\) is independent of \(\theta\), such that \({\rm v}_\| E_\varphi = {\rm v}_\| {\rm EPAR0\_SPITZER}\).

This parameter is used only for the Spitzer problem (SPITZER_MODEL = 1). For the standard neoclassical problem, use EPAR0 instead.

## EQUILIBRIUM_MODEL

**Definition**

Parameter which selects the geometric equilibrium model.

**Choices**

EQUILIBRIUM_MODEL = 0: s-alpha

EQUILIBRIUM_MODEL = 1: large aspect ratio

EQUILIBRIUM_MODEL = 2: Miller

EQUILIBRIUM_MODEL = 3: General Grad-Shafranov

**Comments**

DEFAULT: 0

For experimental profiles (PROFILE_MODEL = 2), this parameter is ignored and the geometric equilibrium model is instead set by the parameter PROFILE_EQUILIBRIUM_MODEL.

EQUILIBRIUM_MODEL=3 is available via interface. For this option, the number of Fourier coefficients, GEO_NY, must be a positive integer, with the corresponding Fourier coefficients set in GEO_YIN. For input.neo, these parameters are set by the file input.geo. Note that in addition to the fourier coefficients, the input equilibrium parameters RMIN_OVER_A, RMAJ_OVER_A, Q, SHEAR, BETA_STAR, BTCCW, and IPCCW must also be specified.

See the geometry notes for more details about the geometric equilibrium models.

## GEO_NY

**Definition**

Number of Fourier coefficients for general Grad-Shafranov equilibrium.

**Comments**

DEFAULT: 0

This parameter is only available via subroutine interface and not by input.neo.

This parameter is used only if EQUILIBRIUM_MODEL = 3. It must be a positive integer. The Fourier coefficient values themselves are specified by GEO_YIN.

See the geometry notes for more details about the general geometry equilibrium model.

## GEO_YIN

**Definition**

Array of dimension (8,0:32) with the normalized Fourier coefficients \(\{a\_R,b\_R,a\_Z,b\_Z\}/a\) and their radial derivatives \(\{a\_Rp,b\_Rp,a\_Zp,b\_Zp\}\) for general Grad-Shafranov equilibrium.

**Comments**

DEFAULT: 0.0

This parameter is only available via subroutine interface and not by input.neo.

This parameter is used only if EQUILIBRIUM_MODEL = 3. The number of Fourier coefficients is specified by GEO_NY and the coefficients are read-in as geo_yin(8,0:geo_ny).

See the geometry notes for more details about the general geometry equilibrium model.

## IPCCW

**Definition**

Parameter which selects the orientation of the plasma current (and thus the poloidal magnetic field \(B_p\)) relative to the toroidal angle \(\varphi\).

**Choices**

IPCCW = 1: Counter-clockwise when viewed from above the torus - negative \(\hat{e}_{\varphi}\) for the right-handed coordinate system \((r,\theta,\varphi)\). Thus, \(B_p\) is oriented along the negative \(\hat{e}_{\varphi}\) direction.

IPCCW = -1: Clockwise when viewed from above the torus - positive \(\hat{e}_{\varphi}\) for the right-handed coordinate system \((r,\theta,\varphi)\). Thus, \(B_p\) is oriented along the positive \(\hat{e}_{\varphi}\) direction.

**Comments**

DEFAULT: -1

In DIII-D, typically IPCCW = 1.

When experimental profiles are used (PROFILE_MODEL = 2), the orientiation of IP is inferred from input.profiles.

## KAPPA

**Definition**

Elongation, \(\kappa\), of the flux surface.

**Comments**

DEFAULT: 1.0

This is only active with EQUILIBRIUM_MODEL = 2 (the Miller equilibrium model).

When experimental profiles are used (PROFILE_MODEL = 2), the elongation as a function of radius is read from input.profiles.

## MASS_*

**Definition**

The normalized mass:

**Commments**

DEFAULT: 1.0

The mass of each species 1-11 is set as: MASS_1, MASS_2, MASS_3,…

When experimental profiles are used (PROFILE_MODEL = 2), the normalizing mass is deuterium, \(m_{\rm norm}=m_{D}\) = 3.3452e-27 kg

The subroutine interface parameter is specified as a vector: neo_mass_in(1:11)

## N_ENERGY

**Definition**

The number of energy polynomials - 1 in the computational domain (\(n_{\varepsilon,\rm total}\) = N_ENERGY+1).

**Comments**

DEFAULT: 6

The velocity-space coordinate \(x_a\) is the normalized velocity: \(x_a = \sqrt{\varepsilon} = {\rm v}/(\sqrt{2}{\rm v}_{ta})\).

NEO uses an expansion of associated Laguerre polynomials in \(x_a\), which is coupled with the Legendre expansion in \(\xi\): \(P_l(\xi) L_m^{k(l)+1/2}(x_a^2)x_a^{k(l)}\), where \(k(l)=0\) for Legendre index \(l=0\) and \(k(l)=1\) for Legendre index \(l>0\).

The collocation integrals are formed from the monomial basis elements, \(x_a^{2m+k(l)}\), which can be written in terms of Gamma and Beta functions.

## N_RADIAL

**Definition**

The number of radial gridpoints, \(n_r\) in the computational domain.

**Comments**

DEFAULT: 1

The radial grid is defined on the range RMIN_OVER_A \(\le r/a \le\) RMIN_OVER_A_2. For a local simulation (PROFILE_MODEL = 1), the normalizing length scale \(a\) is arbitrary. For a global simulation (PROFILE_MODEL = 2), \(a\) is the plasma minor radius at the center of the radial simulation domain.

N_RADIAL > 1 requires a global profile model (PROFILE_MODEL = 2). Otherwise, N_RADIAL = 1 and the profile model is local (PROFILE_MODEL = 1).

For solution of only the first-order DKE, which is a radially-local problem, the radial grid is equally-spaced.

## N_SPECIES

**Definition**

The number of kinetic species.

**Comments**

DEFAULT: 1

The maximum allowed N_SPECIES is 11.

Only one species with charge Z < 0 is allowed. If no species with Z < 0 is specified, then an adiabatic electron model is assumed.

For local simulations (PROFILE_MODEL = 1), the order of the species and the normalizing density and temperature are arbitrary.

For experimental profiles (PROFILE_MODEL = 2), the normalizing mass is the mass of deuterium (\(m_D\) = 3.3452e-27 kg), so the input masses should be given relative to this mass. The output quantities are normalized with respect to the density and temperature of the first species in input.neo and \(m_D\), with \({\rm v}_{\rm norm} = \sqrt{T_{0,{\rm species 1}}/m_{D}}\).

The electron species, if kinetic, must be species number N_SPECIES in input.neo.

Of the species-dependent parameters in input.neo, only Z_* and MASS_* are used, while DENS_*, TEMP_*, DLNNDR_*, DLNTDR_*, and NU_1 are determined from the parameters read from input.profiles.

Quasi-neutrality is checked.

See PROFILE_MODEL for more details.

## N_THETA

**Definition**

The number of theta gridpoints, \(n_\theta\) in the computational domain.

**Comments**

DEFAULT: 17

N_THETA must be an odd number

The theta grid range is equally-spaced and defined on the range \(-\pi \le \theta < \pi\).

The theta derivatives in the kinetic equation are treated with a 4th-order centered finite difference scheme. Periodic boundary conditions are assumed.

## N_XI

**Definition**

The number of xi polynomials - 1 in the computational domain (\(n_{\xi,\rm total}\) = N_XI+1).

**Comments**

DEFAULT: 17

The velocity-space coordinate \(\xi\) is the cosine of the pitch angle: \(\xi ={\rm v}_\|/{\rm v}\).

NEO uses an expansion of Legendre polynomials in \(\xi\).

The collocation integrals are done exactly analytically.

## NE_ADE

**Definition**

The normalized equilibrium-scale density of the electrons for the case of adiabatic electrons:

**Commments**

DEFAULT: 1.0

This paramter is used only if no species with Z < 0 is specified.

When experimental profiles are used (PROFILE_MODEL = 2), the density as a function of radius is read from input.profiles and the normalizing density is the local density of Species 1, \(n_{\rm norm}(r)=n_{{\rm species 1}}\).

## NU_1

**Definition**

The normalized collision frequency of the first kinetic species:

where

**Comments**

DEFAULT: 0.1

Only the collision frequency for Species 1 is specified. The collision frequencies for the other species are computed internally in the code using NU_1, Z_*, MASS_*, DENS_*, and TEMP_*.

When rotation effects are included (ROTATION_MODEL = 2), this parameter is the value at the outboard midplane (\(\theta = 0\)).

When experimental profiles are used (PROFILE_MODEL = 2), this is computed internally from the profile parameters read from input.profiles. Also, the normalizing length scale is the plasma minor radius and the normalizing velocity is \({\rm v}_{\rm norm} = \sqrt{T_{\rm species1}/m_{D}}\).

## OMEGA_ROT

**Definition**

The normalized toroidal angular frequency:

where \(\omega_0=-c\frac{d \Phi_{-1}}{d\psi}\)

**Comments**

DEFAULT: 0.0

Used only if sonic rotation effects are included (ROTATION_MODEL = 2).

When experimental profiles are used (PROFILE_MODEL = 2), the toroidal angular frequency as a function of radius is read from input.profiles. The associated \(E_r\) is assumed to be the lowest-order field, \(E_r^{(-1)}\), and \(E_r^{(0)}\) is assumed to be 0.

## OMEGA_ROT_DERIV

**Definition**

The normalized toroidal rotation shear:

where \(\omega_0=-c\frac{d \Phi_{-1}}{d\psi}\) is the torodial angular frequency.

**Comments**

DEFAULT: 0.0

Used only if sonic rotation effects are included (ROTATION_MODEL = 2).

When experimental profiles are used (PROFILE_MODEL = 2), the toroidal angular frequency as a function of radius is read from input.profiles and its gradient is computed internally. The associated \(E_r\) is assumed to be the lowest-order field, \(E_r^{(-1)}\), and \(E_r^{(0)}\) is assumed to be 0.

## PROFILE_DLNNDR_*_SCALE

**Definition**

Scaling factor for the normalized equilibrium-scale density gradient scale length in profile mode:

**Commments**

DEFAULT: 1.0

The scaling factor of each species 1-11 is set as: PROFILE_DLNNDR_1_SCALE, PROFILE_DLNNDR_2_SCALE, PROFILE_DLNNDR_3_SCALE,…

This parameter is only used when experimental profiles are used (PROFILE_MODEL = 2). The density as a function of radius is read from input.profiles and the density gradient is computed internally. The normalizing length is the plasma minor radius. This gradient scale length is then re-scaled.

The subroutine interface parameter is specified as a vector: neo_profile_dlnndr_scale_in(1:11)

## PROFILE_DLNTDR_*_SCALE

**Definition**

Scaling factor for the normalized equilibrium-scale temperature gradient scale length in profile mode:

**Commments**

DEFAULT: 1.0

The scaling factor of each species 1-11 is set as: PROFILE_DLNTDR_1_SCALE, PROFILE_DLNTDR_2_SCALE, PROFILE_DLNTDR_3_SCALE,…

This parameter is only used when experimental profiles are used (PROFILE_MODEL = 2). The temperature as a function of radius is read from input.profiles and the temperature gradient is computed internally. The normalizing length is the plasma minor radius. This gradient scale length is then re-scaled.

The subroutine interface parameter is specified as a vector: neo_profile_dlntdr_scale_in(1:11)

## PROFILE_EQUILIBRIUM_MODEL

**Definition**

Parameter which selects the geometric equilibrium model for experimental profiles.

**Choices**

PROFILE_EQUILIBRIUM_MODEL = 1: Use Miller shaped geometry with the profiles of the geometric parameters as given in input.profiles.

PROFILE_EQUILIBRIUM_MODEL = 2: Use the general Grad-Shafranov geometry with the fourier coefficients specified in input.profiles.geo.

**Comments**

DEFAULT: 1

Used only for experimental profiles (PROFILE_MODEL = 2)

See the geometry notes for more details about the geometric equilibrium models.

## PROFILE_ERAD0_MODEL

**Definition**

Parameter which selects whether to include \(E_r^{(0)}\) for experimental profiles.

**Choices**

PROFILE_ERAD0_MODEL = 0: \(E_r^{(0)}\) is set to zero regardless of the value in input.profiles.

PROFILE_ERAD0_MODEL = 1: \(E_r^{(0)}\) as specified in input.profiles is used.

**Comments**

DEFAULT: 1

Used only for experimental profiles (PROFILE_MODEL = 2).

If sonic rotation effects are included (ROTATION_MODEL = 2) with experimental profiles, then this parameter is ignored and \(E_r^{(0)}\) is assumed to be zero. This means that the \(E_r^{(0)}\) in input.profiles is assumed to be the lowest-order field in sonic rotation theory, i.e. \(E_r^{(-1)}\),and is used to compute the lowest-order sonic toroidal rotation parameters, OMEGA_ROT and OMEGA_ROT_DERIV.

## PROFILE_MODEL

**Definition**

Parameter which selects how the radial profile is defined.

**Choices**

PROFILE_MODEL = 1: local (one radius)

PROFILE_MODEL = 2: global, using experimental profiles

**Comments**

DEFAULT: 1

For PROFILE_MODEL = 1, N_RADIAL must be 1.

For PROFILE_MODEL = 2, experimental profiles are defined in input.profiles. The number of radial gridpoints is specified by N_RADIAL.

Additional models used for this case are specified by PROFILE_EQUILIBRIUM_MODEL and PROFILE_ERAD0_MODEL.

Of the species-dependent parameters in input.neo, only Z_* and MASS_* are used for this case. The normalizing mass is the mass of deuterium (\(m_D\) = 3.3452e-27 kg), so the input masses should be given relative to this mass. The output quantities are normalized with respect to the density and temperature of the first species in input.neo and \(m_D\), with \({\rm v}_{\rm norm} = \sqrt{T_{0,{\rm species 1}}/m_{D}}\).

The electron species, if kinetic, must be species number N_SPECIES in input.neo.

If the density profiles in input.profiles are not quasi-neutral, then the density profile of the first ion species is re-set.

## Q

**Definition**

Magnitude of the safety factor, \(|q|\), of the flux surface:

**Comments**

DEFAULT: 2.0

When experimental profiles are used (PROFILE_MODEL = 2), the safety factor as a function of radius is read from input.profiles.

The orientation of the safety factor is determined by IPCCW and BTCCW.

## RHO_STAR

**Definition**

The ratio of the Larmor radius of the normalizing species to the normalizing length scale:

**Comments**

DEFAULT: 0.001

This parameter must be a positive number. The sign of \(B_{\rm unit}\) is determined by IPCCW and BTCCW.

When experimental profiles are used (PROFILE_MODEL = 2), \(\rho_*\) is computed internally from the profile parameters in input.profiles and the normalizing length scale is the plasma minor radius.

## RMAJ_OVER_A

**Definition**

The ratio of the flux-surface-center major radius, \(R_0\), to the normalizing length scale:math:a.

**Comments**

DEFAULT: 3.0

When experimental profiles are used (PROFILE_MODEL = 2), the flux-surface-center major radius as a function of radius, \(R_0(r)\) is read from input.profiles and the normalizing length scale is the plasma minor radius.

## RMIN_OVER_A

**Definition**

The ratio of the midplane minor radius \(r\) to the normalizing length scale:math:a.

**Comments**

DEFAULT: 0.5

For N_RADIAL > 1, this parameter is the lower bound of the radial grid.

## RMIN_OVER_A_2

**Definition**

The ratio of the midplane minor radius \(r\) to the normalizing length scale:math:a.

**Comments**

DEFAULT: 0.6

For N_RADIAL > 1, this parameter is the upper bound of the radial grid.

For N_RADIAL = 1, this parameter is not used.

## ROTATION_MODEL

**Definition**

Parameter which selects whether to solve the DKE in the diamagnetic ordering limit or in the sonic toroidal rotation ordering limit.

**Choices**

ROTATION_MODEL = 1: sonic rotation effects not included (diamagnetic ordering assumed)

ROTATION_MODEL = 2: sonic rotation effects included (solves the Hinton-Wong generalized DKE which allows for flow speeds on the order of the thermal speed).

The toroidal rotation frequency OMEGA_ROT and the toroidal rotation shear OMEGA_ROT_DERIV must be specified.

**COMMENTS**

DEFAULT: 1

## S_DELTA

**Definition**

Measure of the rate of change of the average triangularity of the flux surface:

**Comments**

DEFAULT: 0.0

This is only active with EQUILIBRIUM_MODEL = 2 (the Miller equilibrium model).

When experimental profiles are used (PROFILE_MODEL = 2), the triangularity as a function of radius is read from input.profiles and the triangularity gradient is computed internally.

## S_KAPPA

**Definition**

Measure of the rate of change of the elongation of the flux surface:

**Comments**

DEFAULT: 0.0

This is only active with EQUILIBRIUM_MODEL = 2 (the Miller equilibrium model).

When experimental profiles are used (PROFILE_MODEL = 2), the elongation as a function of radius is read from input.profiles and the elongation gradient is computed internally.

## S_ZETA

**Definition**

Measure of the rate of change of the squareness of the flux surface:

**Comments**

DEFAULT: 0.0

This is only active with EQUILIBRIUM_MODEL = 2 (the Miller equilibrium model).

When experimental profiles are used (PROFILE_MODEL = 2), the squareness as a function of radius is read from input.profiles and the squareness gradient is computed internally.

## S_ZMAG

**Definition**

Measure of the rate of change of the elevation of the flux surface:

**Comments**

DEFAULT: 0.0

This is only active with EQUILIBRIUM_MODEL = 2 (the Miller equilibrium model).

When experimental profiles are used (PROFILE_MODEL = 2), the flux-surface elevation as a function of radius, \(Z_0(r)\), is read from input.profiles and its derivative is computed internally.

## SHEAR

**Definition**

Magnetic shear, \(s\), of the flux surface:

**Comments**

DEFAULT: 1.0

NOTE: This parameter is

*not*used in the standard DKE equation! It is only used in the case of an anisotropic temperature species (e.g. ANISO_MODEL_* = 2) to compute \(d\Phi_*/dr\).This is only active with EQUILIBRIUM_MODEL = 2 (the Miller equilibrium model).

When experimental profiles are used (PROFILE_MODEL = 2), the safety factor as a function of radius is read from input.profiles and the safety factor gradient is computed internally.

## SHIFT

**Definition**

Shafranov shift, \(\Delta\), of the flux surface:

**Comments**

DEFAULT: 0.0

This is only active with EQUILIBRIUM_MODEL = 2 (the Miller equilibrium model).

When experimental profiles are used (PROFILE_MODEL = 2), the flux-surface-center major radius as a function of radius, \(R_0(r)\), is read from input.profiles and its derivative is computed internally.

## SHAPE_COS0

**Definition**

0th antisymmetric moment.

DEFAULT = 0.0

## SHAPE_S_COS0

**Definition**

0th antisymmetric moment shear.

DEFAULT = 0.0

## SHAPE_COS1

**Definition**

1st antisymmetric moment.

DEFAULT = 0.0

## SHAPE_S_COS1

**Definition**

1th antisymmetric moment shear.

DEFAULT = 0.0

## SHAPE_COS2

**Definition**

2nd antisymmetric moment.

DEFAULT = 0.0

## SHAPE_S_COS2

**Definition**

2th antisymmetric moment shear.

DEFAULT = 0.0

## SHAPE_COS3

**Definition**

3rd antisymmetric moment.

DEFAULT = 0.0

## SHAPE_S_COS3

**Definition**

3rd antisymmetric moment.

DEFAULT = 0.0

## SHAPE_SIN3

**Definition**

3rd symmetric moment.

DEFAULT = 0.0

## SHAPE_S_SIN3

**Definition**

3rd symmetric moment shear.

DEFAULT = 0.0

## SILENT_FLAG

**Definition**

Parameter which selects how much data to print out.

**Choices**

SILENT_FLAG = 0: output files are written.

SILENT_FLAG > 0: no output files are written.

**Comments**

DEFAULT: 0

## SIM_MODEL

**Definition**

Parameter which selects whether to determine the neoclassical transport from analytic theory or from numerical solution of the DKE.

**Choices**

SIM_MODEL = 0: analytic theory only.

SIM_MODEL = 1: numerical solution and analytic theory and NCLASS.

SIM_MODEL = 2: numerical solution and analytic theory only.

SIM_MODEL = 3: analytic theory and NCLASS only.

SIM_MODEL = 4: neural network of NEO DKE solution.

**Comments**

DEFAULT: 2

## SPITZER_MODEL

**Definition**

Parameter which selects whether to solve the standard neoclassical transport problem or the Spitzer problem.

**Choices**

SPITZER_MODEL = 0: solve the standard neoclassical transport problem.

SPITZER_MODEL = 1: solve the Spitzer problem.

Must be run with an electron species and an ion species.

The Spitzer coefficients (L11, L12, L21, L22) are output in the file out.neo.spitzer.

**Comments**

– DEFAULT: 0

## TE_ADE

**Definition**

The normalized equilibrium-scale temperature of the electrons for the case of adiabatic electrons:

**Commments**

DEFAULT: 1.0

This paramter is used only if no species with Z < 0 is specified.

When experimental profiles are used (PROFILE_MODEL = 2), the temperature as a function of radius is read from input.profiles and the normalizing temperature is the local temperature of Species 1, \(T_{\rm norm}(r)=T_{{\rm species 1}}\).

## TEMP_*

**Definition**

The normalized equilibrium-scale temperature:

**Commments**

DEFAULT: 1.0

The temperature of each species 1-11 is set as: TEMP_1, TEMP_2, TEMP_3,…

When experimental profiles are used (PROFILE_MODEL = 2), the temperature as a function of radius is read from input.profiles and the normalizing temperature is the local temperature of Species 1, \(T_{\rm norm}(r)=T_{{\rm species 1}}\).

The subroutine interface parameter is specified as a vector: neo_temp_in(1:11)

## TEMP_PARA_*

**Definition**

The normalized equilibrium-scale parallel temperature:

**Commments**

DEFAULT: 1.0

The parallel temperature of each species 1-11 is set as: TEMP_PARA_1, TEMP_PARA_2, TEMP_PARA_3,…

This parameter is used only when the species’ anisotropic flag is set (ANISO_MODEL_* = 2).

The subroutine interface parameter is specified as a vector: neo_temp_para_in(1:11)

## TEMP_PERP_*

**Definition**

The normalized equilibrium-scale perpendicular temperature:

**Commments**

DEFAULT: 1.0

The perpendicular temperature of each species 1-11 is set as: TEMP_PERP_1, TEMP_PERP_2, TEMP_PERP_3,…

This parameter is used only when the species’ anisotropic flag is set (ANISO_MODEL_* = 2).

The subroutine interface parameter is specified as a vector: neo_temp_perp_in(1:11)

## THREED_MODEL

**Definition**

Parameter which selects whether to solve the DKE in toroidally axisymmetric limit (2D) or with nonaxisymmetric effects (3D).

**Choices**

THREED_MODEL = 0: toroidally axisymmetric limit (2D).

THREED_MODEL = 1: toroidally nonaxisymmetric effects are included (3D).

This option is presently not available for experimental profiles (PROFILE_MODEL = 2).

The local 3D equilibrium solver LE3 must be run first. All of the equilibrium parameters, including the spatial dimensions for \((\theta,\varphi)\), are read from the LE3 output file.

Of the plasma equilibrium/geometry NEO input paramters, only RHO_STAR, DPHI0DR, and RMIN_OVER_A are used.

Of the numerical resolution NEO input parameters, only N_XI and N_ENERGY are used.

**COMMENTS**

DEFAULT: 0

## THREED_EXB_MODEL

**Definition**

Parameter which selects whether to include the higher-order \({\bf E} \times {\bf B}\) drift velocity in the DKE with nonaxisymmetric effects (3D).

**Choices**

THREED_EXB_MODEL = 0: higher-order \({\bf E} \times {\bf B}\) drift velocity not included.

THREED_EXB_MODEL = 1: higher-order \({\bf E} \times {\bf B}\) drift velocity included.

Used only if toroidal nonaxisymmetric effects are included (THREED_MODEL = 1).

The value of the equilibrium potential in the higher-order \({\bf E} \times {\bf B}\) drift velocity is specified by THREED_EXB_DPHI0DR. Note that this does not affect the equilibrium potential in the neoclassical source term, which is specified by DPHI0DR.

**COMMENTS**

DEFAULT: 0

## THREED_EXB_DPHI0DR

**Definition**

The normalized equilibrium-scale radial electric field in the higher-order \({\bf E} \times {\bf B}\) drift velocity:

such that

**Comments**

DEFAULT: 0.0

Used only if toroidal nonaxisymmetric effects (3D) are included (THREED_MODEL = 1).

This does not affect the equilibrium potential in the neoclassical source term, which is specified by DPHI0DR.

## Z_*

**Definition**

The species’ charge.

**Commments**

DEFAULT: 1.0

The charge of each species 1-11 is set as: Z_1, Z_2, Z_3,…

The subroutine interface parameter is specified as a vector: neo_z_in(1:11)

## ZETA

**Definition**

Squareness, \(\zeta\), of the flux surface.

**Comments**

DEFAULT: 0.0

This is only active with EQUILIBRIUM_MODEL = 2 (the Miller equilibrium model).

When experimental profiles are used (PROFILE_MODEL = 2), the squareness as a function of radius is read from input.profiles.

## ZMAG_OVER_A

**Definition**

The ratio of the elevation of the flux surface, \(Z_0\), to the normalizing length scale \(a\).

**Comments**

DEFAULT: 0.0

When experimental profiles are used (PROFILE_MODEL = 2), the flux-surface elevation as a function of radius, \(Z_0(r)\) is read from input.profiles and the normalizing length scale is the plasma minor radius.

Return to table of inputs