A EQDSK Files
A EQDSK provides various global plasma parameters of interest. The A EQDSK files can be combined into the E type EQDSK files for viewing using the SUMMARD code, which can then be plotted using REVIEW.
- Variables
- Format for A EQDSK may be found in the Fortran source weqdsku.f ("write (neqdsk,...").
- SUMMARD combines timeslice-based A files (Assssss.ttttt) into shot-based E files (Esssss.fit), i.e., E file contains time-dependent A data for shot ssssss.
- EQVIEW - an IDL utility to tabulate and visualize A and E files.
- reada.pro - an IDL routine boundled in 4d that reads A EQDSK into a structure.
Variables
| AAQ1 | minor radius of q=1 surface in cm, 100 if not found |
| AAQ2 | minor radius of q=2 surface in cm, 100 if not found |
| AAQ3 | minor radius of q=3 surface in cm, 100 if not found |
| ALI | li with normalization average poloidal magnetic defined through Ampere's law |
| ALPHA | Shafranov boundary line integral parameter |
| AOUT | plasma minor radius in cm |
| AREAO | cross sectional area in cm2 |
| BCENTR | vacuum toroidal magnetic field in Tesla at RCENCM |
| BETAP | poloidal b with normalization average poloidal magnetic BPOLAV defined through Ampere's law |
| BETAPD | diamagnetic poloidal b |
| BETAT | toroidal b in % |
| BETATD | diamagnetic toroidal b in % |
| BPOLAV | average poloidal magnetic field in Tesla defined through Ampere's law |
| BTAXP | toroidal magnetic field at magnetic axis in Tesla |
| BTAXV | vacuum toroidal magnetic field at magnetic axis in Tesla |
| CCBRSP | computed external coil currents in Ampere |
| CDFLUX | computed diamagnetic flux in Volt-sec |
| CHIGAMT | total chi2 MSE |
| CHIPRE | total chi2 pressure |
| CHIPRE | total chi2 pressure |
| CJ1AVE | normalized average current density in plasma outer 5% normalized poloidal flux region |
| CJOR0 | normalized axial flux surface average current density |
| CJOR95 | normalized flux surface average current density at 95% of normalized poloidal flux |
| CJOR99 | normalized flux surface average current density at 99% of normalized poloidal flux |
| CPASMA | fitted plasma toroidal current in Ampere-turn |
| CPROF | current profile parametrization parameter |
| CSILOP | computed flux loop signals in Weber |
| DCO2R | line average electron density in cm3 from radial CO2 chord |
| DCO2V | line average electron density in cm3 from vertical CO2 chord |
| DOUTL | upper triangularity |
| DOUTU | upper triangularity |
| ECCURT | measured E-coil current in Ampere |
| ELONGM | elongation at magnetic axis |
| EOUT | Plasma boundary elongation |
| FEXPAN | flux expansion at x point |
| FEXPVS | flux expansion at outer lower vessel hit spot |
| FLUXX | measured diamagnetic flux in Volt-sec |
| JFLAG | error flag, 0 for error |
| LFLAG | error flag, > 0 for error |
| LIMLOC | plasma configuration. IN, OUT, TOP, and BOT for limiter configurations limited at inside, outside, top, and bottom. SNT, SNB, and DN for single null top, single null bottom, and double null configurations. MAR for marginally diverted configurations. |
| MCO2V | number of vertical CO2 density chords |
| MCO2R | number of radial CO2 density chords |
| JFLAG | error flag, 0 for error |
| OBOTT | plasma bottom gap in cm |
| OBOTS | bottom gap of external second separatrix in cm |
| OLEFT | plasma inner gap in cm |
| OLEFS | inner gap of external second separatrix in cm |
| ORIGHT | plasma outer gap in cm |
| ORIGHS | outer gap of external second separatrix in cm |
| ORING | not used |
| OTOP | plasma top gap in cm |
| OTOPS | top gap of external second separatrix in cm |
| PASMAT | measured plasma toroidal current in Ampere |
| PBINJ | neutral beam injection power in Watts |
| PP95 | normalized P'(y) at 95% normalized poloidal flux |
| PSURFA | plasma boundary surface area in m2 |
| PSIREF | reference poloidal flux in VS/rad |
| QMERCI | Mercier stability criterion on axial q(0), q(0) > QMERCI for stability |
| QMFLAG | axial q(0) flag, FIX if constrained and CLC for float |
| QOUT | q at plasma boundary |
| QPSIB | q at 95% of poloidal flux |
| QQMAGX | axial safety factor q(0) |
| QQMIN | minimum safety factor qmin |
| QSTA | equivalent safety factor q* |
| RCENCM | major radius in cm for vacuum field BCENTR |
| RCO2R | path length in cm of radial CO2 density chord |
| RCO2V | path length in cm of vertical CO2 density chord |
| RCURRT | major radius in cm of current centroid |
| RMIDIN | inner major radius in m at Z=0.0 |
| RMIDOUT | outer major radius in m at Z=0.0 |
| RMAGX | major radius in cm at magnetic axis |
| RQQMIN | normalized radius of qmin , square root of normalized volume |
| ROUT | majpr radius of geometric center in cm |
| RSEPS | majpr radius of x point in cm |
| RTTT | Shafranov boundary line integral parameter |
| RVSIN | major radius of vessel inner hit spot in cm |
| RVSOUT | major radius of vessel outer hit spot in cm |
| S1, S2, S3 | Shafranov boundary line integrals |
| SEPEXP | separatrix radial expansion in cm |
| SEPLIM | > 0 for minimum gap in cm in divertor configurations, < 0 absolute value for minimum distance to external separatrix in limiter configurations |
| SEPNOS | radial distance in cm between x point and external field line at ZNOSE |
| SHEARB | magnetic shear at 95% enclosed normalized poloidal flux |
| SIBDRY | plasma boundary poloidal flux in Weber/rad |
| SIMAGX | boundary poloidal flux in Weber/rad at magnetic axis |
| SLANTL | gap to lower outboard limiter in cm |
| SLANTU | gap to upper outboard limiter in cm |
| SSEP | outboard radial distance to external second separatrix in cm for single null configurations. positive for double nulls or single top null and negative for single bottom null. |
| SSI01 | magnetic shear at 1% of normalized poloidal flux |
| SSI95 | magnetic shear at 95% of normalized poloidal flux |
| TAUMHD | energy confinement time in ms |
| TAUDIA | diamagnetic energy confinement time in ms |
| TAVEM | average time in ms for magnetic and MSE data |
| TERROR | equilibrium convergence error |
| TIME | time in ms |
| TSAISQ | total chi2 from magnetic probes, flux loops, Rogowskiand external coils |
| VERTN | vacuum field index at current centroid |
| VLOOPT | measured loop voltage in volt |
| VSURFA | plasma surface loop voltage in volt, E EQDSK only |
| VOUT | plasma volume in cm3 |
| WBDOT | time derivative of poloidal magnetic energy in Watt, E EQDSK only |
| WBDOT | time derivative of poloidal magnetic energy in Watt, E EQDSK only |
| WPDOT | time derivative of plasma stored energy in Watt, E EQDSK only |
| WPLASMD | diamagnetic plasma stored energy in Joule |
| XNDNT | plasma boundary indentation |
| XNNC | vertical stability parameter, vacuum field index normalized to critical index value |
| YYY2 | Shafranov Y2 current moment |
| ZCURRT | Z in cm at current centroid |
| ZMAGX | Z in cm at magnetic axis |
| ZOUT | Z of geometric center in cm |
| ZSEPS | Z of x point in cm |
| ZVSIN | Z of vessel inner hit spot in cm |
| ZVSOUT | Z of vessel outer hit spot in cm |
Format
Reading the A eqdsk file in FORTRAN:
parameter (ntime=400)
parameter (magpri67=29,magpri322=31,magprirdp=8)
parameter (magpri=magpri67+magpri322+magprirdp)
character limloc*4,qmflag*3
dimension time(ntime),jflag(ntime)
. ,eout(ntime),rout(ntime),zout(ntime),doutu(ntime)
. ,doutl(ntime),aout(ntime),vout(ntime),betat(ntime),otop(ntime)
. ,betap(ntime),ali(ntime),oleft(ntime),oright(ntime),qsta(ntime)
. ,rcurrt(ntime),zcurrt(ntime),qout(ntime),olefs(ntime)
. ,orighs(ntime),otops(ntime),sibdry(ntime),areao(ntime)
. ,wplasm(ntime),elongm(ntime),qqmagx(ntime),terror(ntime)
. ,rmagx(ntime),zmagx(ntime),obott(ntime),obots(ntime)
. ,alpha(ntime),rttt(ntime),dbpli(ntime),delbp(ntime),oring(ntime)
. ,rseps(2,ntime),zseps(2,ntime),sepexp(ntime),shearb(ntime)
. ,xtch(ntime),ytch(ntime),qpsib(ntime),vertn(ntime),aaq1(ntime)
. ,aaq2(ntime),aaq3(ntime),btaxp(ntime),btaxv(ntime)
. ,simagx(ntime),jerror(ntime),seplim(ntime)
. ,wbpol(ntime),taumhd(ntime),betapd(ntime),betatd(ntime)
. ,alid(ntime),wplasmd(ntime),taudia(ntime),wbpold(ntime)
. ,qmerci(ntime),slantu(ntime),slantl(ntime),zeff(ntime),
. zeffr(ntime),tave(ntime),rvsin(ntime),zvsin(ntime),
. rvsout(ntime),zvsout(ntime),wpdot(ntime),wbdot(ntime),
. vsurfa(ntime),cjor95(ntime),pp95(ntime),ssep(ntime),
. yyy2(ntime),xnnc(ntime),pasman,betatn,psiq1,betat2,
. wtherm(ntime),wfbeam(ntime),taujd3(ntime),tauthn(ntime)
. qsiwant(ntime),cjorsw(ntime),cjor0(ntime),
. ssiwant(ntime),ssi95(ntime),rexpan,fexpan(ntime),qqmin(ntime),
. fexpvs(ntime),shearc,
. sepnose(ntime),ssi01(ntime),znose,rqqmin(ntime)
c
read (neqdsk,1060) time(jj),jflag(jj),lflag,limloc(jj),
. mco2v,mco2r,qmflag
read (neqdsk,1040) tsaisq(jj),rcencm,bcentr(jj),pasmat(jj)
read (neqdsk,1040) cpasma(jj),rout(jj),zout(jj),aout(jj)
read (neqdsk,1040) eout(jj),doutu(jj),doutl(jj),vout(jj)
read (neqdsk,1040) rcurrt(jj),zcurrt(jj),qsta(jj),betat(jj)
read (neqdsk,1040) betap(jj),ali(jj),oleft(jj),oright(jj)
read (neqdsk,1040) otop(jj),obott(jj),qpsib(jj),vertn(jj)
read (neqdsk,1040) (rco2v(k,jj),k=1,mco2v)
read (neqdsk,1040) (dco2v(jj,k),k=1,mco2v)
read (neqdsk,1040) (rco2r(k,jj),k=1,mco2r)
read (neqdsk,1040) (dco2r(jj,k),k=1,mco2r)
read (neqdsk,1040) shearb(jj),bpolav(jj),s1(jj),s2(jj)
read (neqdsk,1040) s3(jj),qout(jj),olefs(jj),orighs(jj)
read (neqdsk,1040) otops(jj),sibdry(jj),areao(jj),wplasm(jj)
read (neqdsk,1040) terror(jj),elongm(jj),qqmagx(jj),cdflux(jj)
read (neqdsk,1040) alpha(jj),rttt(jj),psiref(jj),xndnt(jj)
read (neqdsk,1040) rseps(1,jj),zseps(1,jj),rseps(2,jj)
. ,zseps(2,jj)
read (neqdsk,1040) sepexp(jj),obots(jj),btaxp(jj),btaxv(jj)
read (neqdsk,1040) aaq1(jj),aaq2(jj),aaq3(jj),seplim(jj)
read (neqdsk,1040) rmagx(jj),zmagx(jj),simagx(jj),taumhd(jj)
read (neqdsk,1040,err=380) betapd(jj),betatd(jj),
. wplasmd(jj),diamag(jj)
read (neqdsk,1040,err=380) vloopt(jj),taudia(jj),qmerci(jj),
. tavem(jj)
if (nvernum.ge.970524) then
read (neqdsk,1041,err=380) nsilop0,magpri0,nfcoil0,nesum0
else
nsilop0 = nsilop
nfcoil0 = nfcoil
nesum0 = nesum
if (ishot.lt.91000) then
magpri0 = magpri67+magpri322
else
magpri0 = magpri
endif
endif
read (neqdsk,1040) ((csilop(k,jj),k=1,nsilop0),
. (cmpr2(k,jj),k=1,magpri0))
read (neqdsk,1040) (ccbrsp(k,jj),k=1,nfcoil0)
read (neqdsk,1040) (eccurt(jj,k),k=1,nesum0)
read (neqdsk,1040) pbinj(jj),rvsin(jj),zvsin(jj),rvsout(jj)
read (neqdsk,1040) zvsout(jj),vsurfa(jj),wpdot(jj),wbdot(jj)
read (neqdsk,1040) slantu(jj),slantl(jj),zuperts(jj),chipre
read (neqdsk,1040) cjor95(jj),pp95(jj),ssep(jj),yyy2(jj)
read (neqdsk,1040) xnnc(jj),cprof,oring(jj),cjor0(jj)
read (neqdsk,1040) fexpan(jj),qqmin(jj),chigamt(jj),ssi01(jj)
read (neqdsk,1040) fexpvs(jj),sepnose(jj),ssi95(jj),rqqmin(jj)
if (nvernum.ge.960412) then
read (neqdsk,1040) cjor99(jj),cj1ave(jj),rmidin(jj),rmidout(jj)
endif
if (nvernum.ge.960724) then
read (neqdsk,1040) psurfa(jj),xdum,xdum,xdum
endif
c
1040 format (1x,4e16.9)
1060 format (1h*,f7.2,10x,i5,11x,i5,1x,a3,1x,i3,1x,i3,1x,a3)
.