Initialisation Options

GS2 has a lot of possible options for the initial condition. Of these, only a few are usually used: these are the typical and restarting

Typical options

These options are the most useful for ordinary GS2 simulations:

Less-usual options

There are many more initialisation options for GS2, although the majority are unmaintained and are specific to the studies they were originally written for. Unfortunately, many are also undocumented, so it is difficult to say precisely what they do

  • "alf": Ion distribution function proportional to ; do not use.
  • "all_modes_equal": Initialise with every single parallel and perpendicular mode given the same amplitude. Intended for linear calculations
  • "cgyro": CGYRO-like initialisation, constant amplitude for and zero amplitude for (see ginit_cgyro for details).
  • "cont": Restart, using current value of g in memory; unsupported, do not use.
  • "convect": Unsupported; do not use. Controlled by k0.
  • "gs": Randomly phased sines and cosines in density, parallel velocity, and parallel and perpendicular temperatures. Controlled by:

    • Overall amplitude: phiinit
    • Amplitude of real and imaginary parts of : refac, imfac
    • Density perturbation: den1
    • Parallel velocity perturbation: upar1
    • Parallel temperature perturbation: tpar1
    • Perpendicular temperature perturbation: tperp1
  • "harris": Controlled by k0. Unsupported; do not use.

  • "kpar": Initial perturbation in fields of form , except for parallel velocity which uses is 90 degrees out of phase and uses . Density, parallel velocity, and perpendicular and parallel temperatures can all be controlled independently. Controlled by:

  • "kz0": Constant along a field line, that is initialise only with .

  • "nl": Two modes set to constant amplitude. Controlled by:

    • Overall amplitude: phiinit
    • ky indices of modes: ikk
    • kx indices of modes: itt
    • Parity: even about unless chop_side is true, in which case none
  • "nl2": Two modes proportional to . Controlled by:

    • Overall amplitude: phiinit
    • ky indices of modes: ikk
    • kx indices of modes: itt
    • Parity: even about unless chop_side is true, in which case none
  • "nl3": Two modes with a Gaussian envelope, and perturbation in fields of form , except for parallel velocity which uses is 90 degrees out of phase and uses . Density, parallel velocity, and perpendicular and parallel temperatures can all be controlled independently. Controlled by:

  • "nl3r": Similar to "nl3", but includes apar.

  • "nl4": Unsupported and inconsistent over time; do not use.
  • "nl5": Unsupported; do not use.
  • "nl6": Unsupported; do not use.
  • "nl7": Unsupported; do not use.
  • "ot": Orszag-Tang 2D vortex problem (see aparamp, phiamp, ittt, ikkk)
  • "rh": Maxwellian perturbation in ik == 1 mode; do not use?
  • "stationary_mode": Initialises the distribution function such that it is a time-independent solution to the gyrokinetic equation when collisions are neglected and only one mode (with aky = 0) is present. See Section 4 of the Analytic Geometry Specification documentation for more details.
  • "recon": Unsupported; do not use.
  • "recon3": Unsupported; do not use.
  • "test3": Unsupported; do not use.
  • "xi": Perturbation proportional to pitch angle, . Unlikely to be useful.
  • "xi2": Perturbation proportional to . Unlikely to be useful.
  • "zero": Sets distribution function to zero

Restarting options

These options are used for restarting GS2. See Restarting a simulation for more details.

If GS2 has been built with parallel netCDF, then read_many controls whether a single restart file or one file per processor is used.

You can scale the fields from the restart files with scale. Note that if scale is negative, the fields are scaled by -scale / max(abs(phi))!

Less-usual restarting options

These options are also used for restarting, but have additional behaviours. They are not well documented or maintained, so use at your own risk.

  • "no_zonal": Restart but remove the zonal flow component.
  • "small": Add restarted fields to a background as set by the "noise" option.
  • "smallflat": Similar to "small", but with a different random background.
  • "zonal_only": Restart but set all non-zonal components of the potential and the distribution function to 0. Input parameters:

    • Noise can be added to these other components by setting phiinit > 0.
    • The size of the zonal flows can be adjusted using zf_init.