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 \({\rm 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 \Phi}{d \psi}\) in units of \({\rm rad/s}\) (see plasma rotation).
shape	–
`rmaj(:)`	The generalized major radius, \(R_0(r)\), in units of \({\rm m}\).
`zmag(:)`	Flux-surface elevation, \(Z_0\), in units of \({\rm 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}/{\rm m}^3\).
`te(:)`	The electron temperature, \(T_e\), in units of \({\rm keV}\).
`ptot(:)`	Total plasma pressure, in units of Pascals.
`z_eff(:)`	The (dimensionless) effective ion charge, \(Z_{\rm eff}\).
`ni(:,:)`	Ion density in units of \(10^{19}/{\rm m}^3\). There is a column for every ion species.
`ti(:,:)`	Ion temperature in units of \({\rm keV}\). There is a column for every ion species.
`jbs(:)`	Bootstrap current (parallel) in units of \({\rm MA/m^2}\).
`jrf(:)`	RF-driven current in units of \({\rm MA/m^2}\).
`jnb(:)`	Beam-driven current in units of \({\rm MA/m^2}\).
`jbstor(:)`	Bootstrap current (toroidal) in units of \({\rm MA/m^2}\).
`vtor(:,:)`	Ion toroidal velocity in units of \({\rm m/s}\). There is a column for every ion species.
`vpol(:,:)`	Ion poloidal velocity in units of \({\rm m/s}\). There is a column for every ion species.
powers	–
`qohme(:)`	Ohmic power to electrons in units of \({\rm MW/m^3}\).
`qbeame(:)`	Beam power to electrons in units of \({\rm MW/m^3}\).
`qbeami(:)`	Beam power to ions in units of \({\rm MW/m^3}\).
`qrfe(:)`	RF power to electrons in units of \({\rm MW/m^3}\).
`qrfi(:)`	RF power to ions in units of \({\rm MW/m^3}\).
`qfuse(:)`	Fusion power to electrons in units of \({\rm MW/m^3}\).
`qfusi(:)`	Fusion power to ions in units of \({\rm MW/m^3}\).
`qsync(:)`	Electron synchrotron radiation in units of \({\rm MW/m^3}\).
`qbrem(:)`	Bremsstrahlung radiation in units of \({\rm MW/m^3}\).
`qline(:)`	Electron line radiation in units of \({\rm MW/m^3}\).
`qei(:)`	Electron-ion exchange \({\rm MW/m^3}\).
`qione(:)`	Recombination power to electrons in units of \({\rm MW/m^3}\).
`qioni(:)`	Recombination power to ions in units of \({\rm MW/m^3}\).
`qcxi(:)`	Charge-exchange power to ions in units of \({\rm MW/m^3}\).
particle sources	–
`qpar_beam(:)`	Beam-particle source density in units of \({\rm 1/m^3/s}\).
`qpar_wall(:)`	Wall-particle source density in units of \({\rm 1/m^3/s}\).
momentum source source	–
`qmom(:)`	The total (convected and conducted) torque density in units of \({\rm N/m^2}\).

Important

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