Jump to: navigation, search

Description of Fortran Global Variables

Various Platforms | Directory Structure | Installing EFIT | Logic Flow Chart | Fortran Global Variables


Categories of Global Variables

1. Response arrays
2. Currents computed by EFIT
3. Weight arrays
4. Chisq variables
5. Controlling flags (may be set by user)
  • iacoil
  • ivesel
  • fitsiref
  • ifitvs
  • iecurr
  • ifcurr
  • symmetrize
5.1 Controlling flags (may be set by user) --- controls output
  • iconsi (set to 1 for poloidal contour plots)
  • iexcal (set to 1 to plot measured v. predicted signals)
  • itek (set to 5 to write output to "plotout.dat")
  • keqdsk (determines format of A and G files)
  • kwripre (set to 1 for the "q" output files)
6. ERROR flags set while reading in data
7. Variables that record which constraints are actually used


brsp computed f-coil currents
pcurrtcomputed grid currents
vcurrtcomputed vessel segements currents
ecurrtcomputed e-coil currents
caccurtcomputed a-coil currents
cecurrcomputed e-coil currents
ccbrsp2-d array of computed f-coil currents [computed in routine CHISQ]
chi2rm 1-d array store the value of the largest contribution to chisq at each time point [computed in routine CHISQ]
cpasma 1-d array of plasma current, stored at each time time point [computed in routine CHISQ]
cmpr2 2-d array of computed magnetic detector signals for each time point.
csilop 2-d array of computed flux loop signals for each time point. [computed in routine CHISQ]
expmp2 1-D array contains array of magnetic signals for reading or writing to namelist. Data values are transferred to/from array "expmpi". This variable is not global.
expmpi2-D array array of magnetic signals for each time point.
fitsirefIf .true. , then ???
gmp2pc(m,n)signals at mag det "m" due to a unit current at grid node "n" Computed in routine GREEN.
gsilpc(m,n)signals at flux loop "m" due to a unit current at grid node "n". Computed in routine GREEN.
iacoilIf greater then 0. , then a-coils are included.
idlopc[ntime]1-d time array is set to a value 1 if the diamagnetic loop weight is .gt.0. Otherwise value is 0. Used for informative output. Set in routine DATA_INPUT.
iec[ntime]1-d time array counts the number of e-coils with weighting .gt.0. Used for informative output. Set in routine DATA_INPUT.
iecurr if .eq. 1 : currents from e-coils are fixed at the value(s) read-in (stored in the 1-d array ECURRT).
  • 2 : currents from e-coils are computed by EFIT and stored in the 1-d array CECURR.
  • otherwise: All e-currents are disabled (i.e. not represented).
  • default: IECURR=1
ierpsi[nsilop] contains of error flags for flux loops. No error for loop "i" if ierpsi(i)=0
ifc[ntime]1-d time array counts the number of f-coils with weighting .gt.0. Used for informative output. Set in routine DATA_INPUT.
ifitvs If .true. , then ???
ifcurr if .gt.0 : then programme exits from routine FCURRT at start of this routine. This is useful for running EFIT in equilibrium mode, and in the absense of a plasma, where we just want to compute signals due to fixed f-coil currents.
  • otherwise: programme goes all the way though routine FCURRT.
  • default: IFCURR=0
imag2[ntime]1-d time array counts the number of magnetic detectors with weighting .gt.0. Used for informative output. Set in routine DATA_INPUT.
iplasm[ntime]1-d time array is set to a value 1 if the plasma current weight is .gt.0. Otherwise value is 0. Used for informative output. Set in routine DATA_INPUT.
ipsi[ntime]1-d time array counts the number of flux loops with weighting .gt.0. Used for informative output. Set in routine DATA_INPUT.
ivesel If greater then 0. , then vessel segements are included.
kdata Mode of operating EFIT, supllied by user.
kstark Variable counts the number of mse measurements at this time slice with weighting .gt.0. Used by many routines, but not output to fitout.dat file like arrays IPSI, IFC, IEC and IMAG.
nqaxis Variable is set to a value 1 if the weight constraining q(0) is .gt.0. Otherwise value is 0. Used for informative output. Set in routine DATA_INPUT.
rmp2ac(m,n)signals at mag det "m" due to a unit current in a-coil "n" Computed in routine GREEN.
rmp2ec(m,n)signals at mag det "m" due to a unit current in e-coil (iecurr=1) or c-coil (iecurr=2) "n". Computed in routine GREEN.
rmp2fc(m,n)signals at mag det "m" due to a unit current in f-coil "n" Computed in routine GREEN.
rmp2pc Response array, computed in routine GREEN.
rmp2vs(m,n) signals at mag det "m" due to a unit current in vessel segement "n" Computed in routine GREEN.
rbdrvsResponse array, computed in routine GREEN.
rgampc Response array, computed in routine GREEN.
rgamvs Response array, computed in routine GREEN.
rsilac(m,n)signals at flux loop "m" due to a unit current in a-coil "n" Computed in routine GREEN.
rsilec(m,n)signals at flux loop "m" due to a unit current in e-coil (iecurr=1) or c-coil (iecurr=2) "n". Computed in routine GREEN.
rsilfc(m,n)signals at flux loop "m" due to a unit current in f-coil "n". Computed in routine GREEN.
rsilpc Response array, computed in routine GREEN.
rsilvs(m,n) signals at flux loop "m" due to a unit current in vessel segment "n". Computed in routine GREEN.
saiec 1-d array of chisq for each e-coil current (coils with zero weight have their chisq set to zero) [computed in routine CHISQ]
saifc 1-d array of chisq for each f-coil current (coils with zero weight have their chisq set to zero) [computed in routine CHISQ]
saiip variable contains the chisq of the plasma current [computed in routine CHISQ]
saimpi1-d array of chisq for each magnetic detector (detectors with zero weight have their chisq set to zero) [computed in routine CHISQ]
saisil1-d array of chisq for each flux loop (flux loops with zero weight have their chisq set to zero) [computed in routine CHISQ]
saisq this variable contains the sum of all chisq contributions (NSILOP flux loops + MAGPRI magnetic detectors + plasma current [computed in routine CHISQ]
swrcur saved weight of plasma current?
swtmp2saved weight of mag detectors?
swtsi saved weight of flux loops?
symmetrizeThis options, set by the user in namelist "INK", averages the computed flux values above and below the mid-line of the grid. However, it doesn't symmetrize the response functions, therefore,it is not a properly symmetrizing facility i. (Default value=.false.)