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*