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 HAM equilibrium parameters.

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)\). This should span the region \({\hat\rho} \in [0,1]\) using n_grid_exp values. The grid-spacing can be nonuniform.

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}\). For more information, see plasma rotation.

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(:)

The dimensionless plasma elongation, \(\kappa\).

delta(:)

The dimensionless plasma triangularity, \(\delta\).

zeta(:)

Plasma squareness, \(\zeta\).

shape_cos0(:)

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}\).

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.

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}\).