Tabular list for input.tglf and outputs

Control

input.tglf parameter

Interface parameter

Short description

Default

NS

tglf_ns_in

number of species including both electrons and ions

2

USE_TRANSPORT_MODEL

tglf_use_transport_model_in

.true.

GEOMETRY_FLAG

tglf_geometry_flag_in

geometry type (0= \(s-\alpha\) , 1=Miller, 2=Fourier, 3=ELITE)

1

USE_BPER

tglf_use_bper_in

include transverse magnetic fluctuations, \(\delta A_{\lVert }\)

.false.

USE_BPAR

tglf_use_bpar_in

include compressional magnetic fluctuations, \(\delta B_{\lVert }\)

.false.

USE_BISECTION

tglf_use_bisection_in

use bisection search method to find width that maximizes growth rate

.true.

USE_MHD_RULE

tglf_use_mhd_rule_in

ignore pressure gradient contribution to curvature drift

.true.

USE_INBOARD_DETRAPPED

tglf_use_inboard_detrapped_in

Set trapped fraction to zero if eigenmode is inward ballooning

.false.

SAT_RULE

tglf_sat_rule_in

0=default saturation rule

0

KYGRID_MODEL

tglf_kygrid_model_in

1=standard ky spectrum for transport model, 0=user defined with NKY modes up to KY equal spaced

1

XNU_MODEL

tglf_xnu_model_in

Collision model (2=new)

2

VPAR_MODEL

tglf_vpar_model_in

0=low-Mach-number limit

0

VPAR_SHEAR_MODEL

tglf_vpar_shear_model_in

deprecated parameter

0

SIGN_BT

tglf_sign_bt_in

sign of \(B_{T}\) with respect to CCW toroidal direction from top

1.0

SIGN_IT

tglf_sign_it_in

sign of \(I_{T}\) with respect to CCW toroidal direction from top

1.0

KY

tglf_ky_in

\(k_{\theta }\rho _{s,{\rm {unit}}}\,\!\) for single-mode call to TGLF

0.3

NEW_EIKONAL

tglf_new_eikonal_in

.true. = compute the eikonal, .false. = use the eikonal computed on the last call to TGLF made with tglf_new_eikonal_in = .true.

.true.

VEXB

tglf_vexb_in

normalized of ExB velocity Doppler shift common to all species (not in use, see VPAR)

0.0

VEXB_SHEAR

tglf_vexb_shear_in

normalized toroidal ExB velocity Doppler shift gradient common to all species. For large ExB velocity ordering \(V_{tor}=V_{ExB}\) \(-SIGN(I_{tor}){\frac {r}{ABS(q)}}{\frac {\partial }{\partial r}}({\frac {V_{ExB}}{R}})\frac {a}{c_{s}}\)

0.0

BETAE

tglf_betae_in

\(\beta _{e}\,\!\) defined with respect to \(B_{\rm {unit}}\,\!\)

0.0

XNUE

tglf_xnue_in

electron-ion collision frequency \({\frac {v_{ei}}{c_{s}/a}}\)

0.0

ZEFF

tglf_zeff_in

effective ion charge

1.0

DEBYE

tglf_debye_in

Debye length/gyroradius

0.0

IFLUX

tglf_iflux_in

compute quasilinear weights and mode amplitudes

.true.

IBRANCH

tglf_ibranch_in

0 = find two most unstable modes one for each sign of frequency, electron drift direction (1), ion drift direction (2), -1 = sort the unstable modes by growthrate in rank order

-1

NMODES

tglf_nmodes_in

number of modes to store for tglf_ibranch_in = -1

2

NBASIS_MAX

tglf_nbasis_max_in

maximum number of parallel basis functions

4

NBASIS_MIN

tglf_nbasis_min_in

minimum number of parallel basis functions

2

NXGRID

tglf_nxgrid_in

number of nodes in Gauss-Hermite quadrature

16

NKY

tglf_nky_in

number of poloidal modes in the high-k spectrum of TGLF_TM

12

ADIABATIC_ELEC

tglf_adiabatic_elec_in

use adiabatic electrons

.false.

ALPHA_P

tglf_alpha_p_in

multiplies parallel velocity shear for all species

1.0

ALPHA_MACH

tglf_alpha_mach_in

multiplies parallel velocity for all species

0.0

ALPHA_E

tglf_alpha_e_in

multiplies ExB velocity shear for spectral shift model

1.0

ALPHA_QUENCH

tglf_alpha_quench_in

1.0 = use quench rule, 0.0 = use new spectral shift model

0.0

XNU_FACTOR

tglf_xnu_factor_in

multiplies the trapped/passing boundary electron-ion collision terms

1.0

DEBYE_FACTOR

tglf_debye_factor_in

multiplies the debye length

1.0

ETG_FACTOR

tglf_etg_factor_in

exponent for ETG saturation rule

1.25

WRITE_WAVEFUNCTION_FLAG

tglf_write_wavefunction_flag_in

Self-explanatory

0

Species Vectors

input.tglf parameter

Interface parameter

Short description

Default

ZS_1

tglf_zs_in(:)

species charge numbers

-1.0,1.0

MASS_1

tglf_mass_in(:)

species masses normalized to \(m_{D}\)

me/md,1.0

RLNS_1

tglf_rlns_in(:)

species array of normalized density gradients \(-{\frac {a}{n_{s}}}{\frac {dn_{s}}{dr}}\)

1.0,1.0

RLTS_1

tglf_rlts_in(:)

species array of normalized temperature gradients \(-{\frac {a}{T_{s}}}{\frac {dT_{s}}{dr}}\)

3.0,3.0

TAUS_1

tglf_taus_in(:)

species array of temperatures \({\frac {T_{s}}{T_{e}}}\)

1.0,1.0

AS_1

tglf_as_in(:)

species array of densities \({\frac {n_{s}}{n_{e}}}\)

1.0,1.0

VPAR_1

tglf_vpar_in(:)

species array of parallel velocities \(SIGN(I_{tor}){\frac {R_{maj}V_{tor}}{Rc_{s}}}\)

0.0,0.0

VPAR_SHEAR_1

tglf_vpar_shear_in(:)

normalized parallel velocity gradient \(-SIGN(I_{tor})R_{maj}{\frac {\partial }{\partial r}}({\frac {V_{tor}}{R}}){\frac {a}{c_{s}}}\)

0.0,0.0

Gaussian width parameters

input.tglf parameter

Interface parameter

Short description

Default

WIDTH

tglf_width_in

maximum width of the Gaussian measure for the parallel Hermite polynomial basis

1.65

WIDTH_MIN

tglf_width_min_in

minimum width used in search for maximum growth rate

0.3

NWIDTH

tglf_nwidth_in

maximum number of widths used in search for maximum growth rate

21

FIND_WIDTH

tglf_find_width_in

.true. = find the width that maximizes the growth rate, .false. = use width

.true.

Miller geometry parameters

input.tglf parameter

Interface parameter

Short description

Default

RMIN_LOC

tglf_rmin_loc_in

flux surface centroid minor radius \(r/a\,\!\)

0.5

RMAJ_LOC

tglf_rmaj_loc_in

flux surface centroid major radius \(R_{maj}/a\,\!\)

3.0

ZMAJ_LOC

tglf_zmaj_loc_in

flux surface centroid elevation \(Z_{maj}/a\,\!\)

0.0

Q_LOC

tglf_q_loc_in

absolute value of the safety factor, \(ABS(q)\,\!\)

2.0

Q_PRIME_LOC

tglf_q_prime_loc_in

\({\frac {q^{2}a^{2}}{r^{2}}}s\)

16.0

P_PRIME_LOC

tglf_p_prime_loc_in

\({\frac {qa^{2}}{rB_{unit}^{2}}}{\frac {\partial p}{\partial r}}\)

0.0

DRMINDX_LOC

tglf_drmindx_loc_in

allows for x different than r \({\frac {\partial r}{\partial x}}\)

1.0

DRMAJDX_LOC

tglf_drmajdx_loc_in

\({\frac {\partial R_{maj}}{\partial x}}\)

0.0

DZMAJDX_LOC

tglf_dzmajdx_loc_in

\({\frac {\partial Z_{maj}}{\partial x}}\)

0.0

KAPPA_LOC

tglf_kappa_loc_in

elongation of flux surface, \(\kappa \,\!\)

1.0

S_KAPPA_LOC

tglf_s_kappa_loc_in

shear in elongation, \({\frac {r}{\kappa }}{\frac {\partial \kappa }{\partial r}}\)

0.0

DELTA_LOC

tglf_delta_loc_in

0.0

S_DELTA_LOC

tglf_s_delta_loc_in

shear in triangularity, \(r{\frac {\partial \delta }{\partial r}}\)

0.0

ZETA_LOC

tglf_zeta_loc_in

squareness, \(\zeta \,\!\) , of flux surface

0.0

S_ZETA_LOC

tglf_s_zeta_loc_in

shear in squareness, \(r{\frac {\partial \zeta }{\partial r}}\)

0.0

KX0_LOC

tglf_kx0_in

kx0/ky ballooning mode offset

0.0

s-alpha geometry parameters

input.tglf parameter

Interface parameter

Short description

Default

RMIN_SA

tglf_rmin_sa_in

normalized minor radius of flux surface \(r/a\,\!\)

0.5

RMAJ_SA

tglf_rmaj_sa_in

normalized major radius of flux surface \(R_{maj}/a\,\!\)

3.0

Q_SA

tglf_q_sa_in

absolute value of safety factor

2.0

SHAT_SA

tglf_shat_sa_in

magnetic shear \({\frac {r}{q}}{\frac {\partial q}{\partial r}}\)

1.0

ALPHA_SA

tglf_alpha_sa_in

normalized pressure gradient

0.0

XWELL_SA

tglf_xwell_sa_in

magnetic well

0.0

THETA0_SA

tglf_theta0_sa_in

\(\theta _{0}={\frac {k_{x}}{sk_{y}}}\)

0.0

B_MODEL_SA

tglf_b_model_sa_in

0/1 to exclude/include the B(theta) factor in k_per

1

FT_MODEL_SA

tglf_ft_model_sa_in

1 uses trapped fraction at the outboard midplane

1

Change-at-your-own-risk parameters

input.tglf parameter

Interface parameter

Short description

Default

THETA_TRAPPED

tglf_theta_trapped_in

parameter to adjust trapped fraction model

0.7

PARK

tglf_park_in

multiplies the parallel gradient term

1.0

GHAT

tglf_ghat_in

multiplies the curvature drift closure terms

1.0

GCHAT

tglf_gchat_in

multiplies the curvature drift irreducible terms

1.0

WD_ZERO

tglf_wd_zero_in

cutoff for curvature drift eigenvalues to prevent zero

0.1

LINSKER_FACTOR

tglf_linsker_factor_in

multiplies the Linsker terms

0.0

GRADB_FACTOR

tglf_gradB_factor_in

multiplies the gradB terms

0.0

FILTER

tglf_filter_in

sets threshold for frequency/drift frequency to filter out non-driftwave instabilities

2.0

TGLF output parameters

Interface Parameter

Definition

Description

tglf_elec_pflux_out

\(\overline{\Gamma_e}/\Gamma_{GB}\)

Gyrobohm-normalized electron particle flux density

tglf_elec_mflux_out

\(\overline{\Pi_e}/\Pi_{GB}\)

Gyrobohm-normalized electron angular momentum flux density

tglf_elec_eflux_out

\(\overline{Q_e}/Q_{GB}\)

Gyrobohm-normalized electron energy flux density

tglf_ion1_pflux_out

\(\overline{\Gamma_1}/\Gamma_{GB}\)

Gyrobohm-normalized main ion particle flux density

tglf_ion1_mflux_out

\(\overline{\Pi_1}/\Pi_{GB}\)

Gyrobohm-normalized main ion angular momentum flux density

tglf_ion1_eflux_out

\(\overline{Q_1}/Q_{GB}\)

Gyrobohm-normalized main ion energy flux density

tglf_ion2_pflux_out

\(\overline{\Gamma_2}/\Gamma_{GB}\)

Gyrobohm-normalized 2nd ion particle flux density

tglf_ion2_mflux_out

\(\overline{\Pi_2}/\Pi_{GB}\)

Gyrobohm-normalized 2nd ion angular momentum flux density

tglf_ion2_eflux_out

\(\overline{Q_2}/Q_{GB}\)

Gyrobohm-normalized 2nd ion energy flux density