This module contains the subroutines which set the initial value of the
| Type | Visibility | Attributes | Name | Initial | |||
|---|---|---|---|---|---|---|---|
| integer, | private | :: | ginitopt_switch | ||||
| integer, | private, | parameter | :: | ginitopt_default | = | 1 | |
| integer, | private, | parameter | :: | ginitopt_noise | = | 2 | |
| integer, | private, | parameter | :: | ginitopt_restart_many | = | 3 | |
| integer, | private, | parameter | :: | ginitopt_kpar | = | 4 | |
| integer, | private, | parameter | :: | ginitopt_nltest | = | 5 | |
| integer, | private, | parameter | :: | ginitopt_kxtest | = | 6 | |
| integer, | private, | parameter | :: | ginitopt_rh | = | 7 | |
| integer, | private, | parameter | :: | ginitopt_remap | = | 8 | |
| real, | public | :: | width0 | ||||
| real, | public | :: | phiinit | ||||
| real, | private | :: | imfac | ||||
| real, | private | :: | refac | ||||
| real, | private | :: | zf_init | ||||
| real, | private | :: | den0 | ||||
| real, | private | :: | upar0 | ||||
| real, | private | :: | tpar0 | ||||
| real, | private | :: | tperp0 | ||||
| real, | private | :: | den1 | ||||
| real, | private | :: | upar1 | ||||
| real, | private | :: | tpar1 | ||||
| real, | private | :: | tperp1 | ||||
| real, | private | :: | den2 | ||||
| real, | private | :: | upar2 | ||||
| real, | private | :: | tpar2 | ||||
| real, | private | :: | tperp2 | ||||
| real, | public | :: | tstart | ||||
| real, | private | :: | scale | ||||
| real, | private | :: | kxmax | ||||
| real, | private | :: | kxmin | ||||
| logical, | private | :: | chop_side | ||||
| logical, | private | :: | left | ||||
| logical, | public | :: | scale_to_phiinit | ||||
| logical, | private | :: | oddparity | ||||
| character(len=300), | public | :: | restart_file | ||||
| character(len=150), | private | :: | restart_dir | ||||
| logical, | private | :: | initialized | = | .false. | ||
| logical, | private | :: | exist |
if maxwwellian_normalization = .true., the pdf is normalized by F0 (which is not the case otherwise) unless reading in g from a restart file, normalise g by F0 for a full flux surface simulation
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| logical, | intent(out) | :: | restarted | |||
| integer, | intent(out) | :: | istep0 |
Initialise the distribution function with random noise. This is the default
gvmu is initialised with a Maxwellian weighting for flux tube simulations, with the Maxwellian evaluated at ia = 1 we are undoing that weighting here, so also need to use ia = 1