FLUX

Short Description

Program to calculate equilibrium data needed by EPEC.

Keywords

equilibrium

Long Description

Program to read gFile(s) and write plasma equilibrium data needed by NEOCLASSICAL and PHASE to fFile(s).

Stage 1:

Read gFile(s) and output plasma equilibrium data for Stage 2

Stage 2:

Construct flux coordinate system. Calculate metric quantities. Locate rational surfaces. Calculate neoclassical data. Calculate tearing stability matrix. Calculate GGJ data. Calculate island saturation data.

Contents

Inputs:
Flux.nml:

Fortran_90 namelist control file

gFile:

Equilibrium EQDSK gFile

gFiles:
Directory containing gFiles for interpolation:
Index:

List of gFile names and experimental times

gFiles:

Actual gFiles

Outputs:
fFile:

File containing data for NEOCLASSICAL and PHASE

fFiles:
Directory containing fFiles associated with interpolated gFiles
Index:

List of fFile names and experimental times

fFiles:

Actual fFiles

Stage1.nc:

NETCDF data file from Stage 1 calculation

Stage2.nc:

NETCDF data file from Stage 2 calculation

Plots:
*.py:

Python scripts to plot Stage 1 and Stage2 data

gFile Format

read (100, ‘(a48, 3i4)’) string, i3, NRBOX, NZBOX
read (100, ‘(5e16.9)’) RBOXLEN, ZBOXLEN, R0, RBOXLFT, ZOFF
read (100, ‘(5e16.9)’) RAXIS, ZAXIS, PSIAXIS, PSIBOUND, B0
read (100, ‘(5e16.9)’) CURRENT, zero, zero, zero, zero
read (100, ‘(5e16.9)’) zero, zero, zero, zero, zero
read (100, ‘(5e16.9)’) (T (i), i = 1, NRBOX)
read (100, ‘(5e16.9)’) (P (i), i = 1, NRBOX)
read (100, ‘(5e16.9)’) (TTp (i), i = 1, NRBOX)
read (100, ‘(5e16.9)’) (Pp (i), i = 1, NRBOX)
read (100, ‘(5e16.9)’) ((PSI (i, j), i = 1, NRBOX), j = 1, NZBOX)
read (100, ‘(5e16.9)’) (Q (i), i = 1, NRBOX)
read (100, ‘(2i5)’) NBOUND, NLIM
read (100, ‘(5e16.9)’) (RBOUND (i), ZBOUND (i), i = 1, NBOUND)
read (100, ‘(5e16.9)’) (RLIM (i), ZLIM (i), i = 1, NLIM)

Anything after this is ignored