RESCALE

Short Description

Program to rescale equilibrium and profile data needed by EPEC.

Keywords

rescale, equilibrium

Long Description

Program to read gFile, pFile, and cFile, rescale equilibrium, and write out new gFile, pFile, and cFile

Contents

/Inputs:

Rescale.nml:

Fortran_90 namelist control file

gFile:

Initial gFile

pFile:

Initial pFile

cFile:

Initial cFile

/Outputs:

gFile:

Rescaled gFile

pFile:

Rescaled pFile

cFile

Rescaled cFile

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

pFile Format

n “psinorm ne(10^20/m^3) dnedpsiN”

for (int i = 0; i < n; i++)

PSI, NE, dNEdPSI

n “psinorm te(KeV) dtedpsiN”

for (int i = 0; i < n; i++)

PSI, TE, dTEdPSI

n “psinorm ni(10^20/m^3) dnidpsiN”

for (int i = 0; i < n; i++)

PSI, NI, dNIdPSI

n “psinorm ti(KeV) dtidpsiN”

for (int i = 0; i < n; i++)

PSI, TI, dTIdPSI

n “psinorm nb(10^20/m^3) dpbdpsiN”

for (int i = 0; i < n; i++)

PSI, NB, dNBdPSI

n “psinorm pb(kPa) dnbdpsiN”

for (int i = 0; i < n; i++)

PSI, PB, dPBdPSI

n “psinorm ptot(kPa) dptotdpsiN”

for (int i = 0; i < n; i++)

PSI, PTOT, dPTOTdPSI

n “psinorm omeg(kRad/s) domeg/dpsiN”

for (int i = 0; i < n; i++)

PSI, WPHI, dWPHIdPSI

n “psinorm omegp(kRad/s) domegp/dpsiN”

for (int i = 0; i < n; i++)

PSI, WTHE, dWTHEdPSI

n “psinorm omegvb(kRad/s) domevb/dpsiN”

for (int i = 0; i < n; i++)

PSI, WVB, dWVBdPSI

n “psinorm omegpp(kRad/s) domepp/dpsiN”

for (int i = 0; i < n; i++)

PSI, WPP, dWPPdPSI

n “psinorm omgeb(kRad/s) domgeb/dpsiN”

for (int i = 0; i < n; i++)

PSI, WEB, dWEBdPSI

n “psinorm er(kV/m) der/dpsiN”

for (int i = 0; i < n; i++)

PSI, ER, dERdPSI

n “psinorm ommvb(kRad/s) dommvb/dpsiN”

for (int i = 0; i < n; i++)

PSI, WMVB, dWMVBdPSI

n “psinorm ommpp(kRad/s) dommpp/dpsiN”

for (int i = 0; i < n; i++)

PSI, WMPP, dWMPPdPSI

n “psinorm omevb(kRad/s) domevb/dpsiN”

for (int i = 0; i < n; i++)

PSI, WEVB, dWEVBdPSI

n “psinorm omepp(kRad/s) domepp/dpsiN”

for (int i = 0; i < n; i++)

PSI, WEPP, dWEPPdPSI

n “psinorm kpol(km/s/T) dkpol/dpsiN”

for (int i = 0; i < n; i++)

PSI, KPOL, dKPOLdPSI

n “psinorm omghb() domghb/dpsiN”

for (int i = 0; i < n; i++)

PSI, WMGB, dWMGBdPSI

n “psinorm nz1(10^20/m^3) dnz1dpsiN”

for (int i = 0; i < n; i++)

PSI, NI, dNIdPSI

n “psinorm vtor1(km/s) dvtor1psiN”

for (int i = 0; i < n; i++)

PSI, VTOR1, dVTOR1dPSI

n “psinorm vpol1(km/s) dvpol1psiN”

for (int i = 0; i < n; i++)

PSI, VPOL1, dVPOL1dPSI

n “N Z A of ION SPECIES”

for (int i = 0; i < n; i++)

N, Z, A (i=0 impurity, i=1 majority; i=2 fast)

PSI:

Normalized poloidal flux

NE:

Electron number density (10^20/m^3)

TE:

Electron temperature (keV)

NI:

Thermal ion number density (10^20/m^3)

TI:

Thermal ion temperature (keV)

NB:

Fast ion number density (10^20/m^3)

WPHI:

Impurity ion toroidal angular velocity on outboard midplane (krad/s)

WTHE:

Impurity ion toroidal angular velocity on outboard midplane (krad/s)

WEB:

ExB frequency (krad/s)

NI:

Impurity ion number density (10^20/m^3)

N:

Ion atomic number

Z:

Ion charge (units of e)

A:

Ion mass number

Fields can occur in any order. Additional fields are ignored.

cFile Format

n

for (int i = 0; i < n; i++)

PSI, CHI_PHI, CHI_E, D_PERP, CHI_I

PSI:

Normalized poloidal flux

CHI_PSI:

Perpendicular toroidal momentum diffusivity (m^2/s)

CHI_E:

Perpendicular electron energy diffusivity (m^2/s)

D_PERP:

Perpendicular particle diffusivity (m^2/s)

CHI_I:

Perpendicular ion energy diffusivity (m^2/s)