input.gacode

To see what a sample input.gacode looks like, at the command line type

$ cgyro -g reg14 ; cat reg14/input.gacode

Overview

The file input.gacode contains the entire dataset required for specification of experimental profiles. All such profiles are specified on an nexp-point grid. The information included is sufficient to carry out simulations of strongly-shaped up-down asymmetric (i.e., arbitrary) equilibria using the new MXH equilibrium parameters [ACB20].

Profile ordering in input.gacode is arbitrary and comment lines (starting with # in the first column) can be added for convenience. These comment lines are ignored by the parser. To generate input.gacode, please use the profiles_gen command-line tool.

File structure

parameter	description
nion	Total number of ions (thermal and fast).
nexp	Number of experimental data gridpoints.
rho(:)	The dimensionless ONETWO flux-surface label, $\hat{\rho }=\rho (r)/\rho (a)\in [0,1]$ (nonuniform allowed).
rmin(:)	The generalized minor radius, $r$, in units of $\mathrm{m}$. See here for definition.
polflux(:)	Poloidal flux over $2\pi$, in units of Webers/radian.
q(:)	Safety factor, $q$.
omega0(:)	Rotation frequency, ${\omega }_0={\displaystyle \frac{c{E}_r}{R{B}_p}=-c\frac{d\mathrm{\Phi }}{d\psi }}$ in units of $\mathrm{r}\mathrm{a}\mathrm{d}/\mathrm{s}$ (see plasma rotation).
shape	–
rmaj(:)	The generalized major radius, $R_0(r)$, in units of $\mathrm{m}$.
zmag(:)	Flux-surface elevation, $Z_0$, in units of $\mathrm{m}$.
kappa(:)	Flux-surface elongation, $\kappa$.
delta(:)	Flux-surface triangularity, $\delta$.
zeta(:)	Flux-surface squareness, $\zeta$.
shape_cos0(:)	Flux-surface tilt.
shape_cos1(:)
shape_cos2(:)
shape_cos3(:)
shape_sin3(:)
	–
ne(:)	The electron density, $n_e$, in units of ${10}^{19}/{\mathrm{m}}^3$.
te(:)	The electron temperature, $T_e$, in units of $\mathrm{k}\mathrm{e}\mathrm{V}$.
ptot(:)	Total plasma pressure, in units of Pascals.
z_eff(:)	The (dimensionless) effective ion charge, $Z_{\mathrm{e}\mathrm{f}\mathrm{f}}$.
ni(:,:)	Ion density in units of ${10}^{19}/{\mathrm{m}}^3$. There is a column for every ion species.
ti(:,:)	Ion temperature in units of $\mathrm{k}\mathrm{e}\mathrm{V}$. There is a column for every ion species.
jbs(:)	Bootstrap current (parallel) in units of $\mathrm{M}\mathrm{A}/{\mathrm{m}}^2$.
jrf(:)	RF-driven current in units of $\mathrm{M}\mathrm{A}/{\mathrm{m}}^2$.
jnb(:)	Beam-driven current in units of $\mathrm{M}\mathrm{A}/{\mathrm{m}}^2$.
jbstor(:)	Bootstrap current (toroidal) in units of $\mathrm{M}\mathrm{A}/{\mathrm{m}}^2$.
vtor(:,:)	Ion toroidal velocity in units of $\mathrm{m}/\mathrm{s}$. There is a column for every ion species.
vpol(:,:)	Ion poloidal velocity in units of $\mathrm{m}/\mathrm{s}$. There is a column for every ion species.
powers	–
qohme(:)	Ohmic power to electrons in units of $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
qbeame(:)	Beam power to electrons in units of $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
qbeami(:)	Beam power to ions in units of $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
qrfe(:)	RF power to electrons in units of $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
qrfi(:)	RF power to ions in units of $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
qfuse(:)	Fusion power to electrons in units of $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
qfusi(:)	Fusion power to ions in units of $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
qsync(:)	Electron synchrotron radiation in units of $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
qbrem(:)	Bremsstrahlung radiation in units of $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
qline(:)	Electron line radiation in units of $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
qei(:)	Electron-ion exchange $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
qione(:)	Recombination power to electrons in units of $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
qioni(:)	Recombination power to ions in units of $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
qcxi(:)	Charge-exchange power to ions in units of $\mathrm{M}\mathrm{W}/{\mathrm{m}}^3$.
particle sources	–
qpar_beam(:)	Beam-particle source density in units of $1/{\mathrm{m}}^3/\mathrm{s}$.
qpar_wall(:)	Wall-particle source density in units of $1/{\mathrm{m}}^3/\mathrm{s}$.
momentum source source	–
qmom(:)	The total (convected and conducted) torque density in units of $\mathrm{N}/{\mathrm{m}}^2$.

Important

For more information about how powers are summed in TGYRO, please see Scenarios and connection to energy sources.