NEOCLASSICAL

Short Description

Program to generate profile data needed by EPEC.

Keywords

profile, neoclassical

Long Description

Program to read FLUX data from fFile(s), profile data from pFile(s) and cFiles(s), and write neoclassical/profile data needed by PHASE to nFile(s)

Stage 3:

Read input files, calculate neoclassical/profile data at resonant surfaces, and output nFile(s)

Contents

Inputs:

Neoclassical.nml:

Fortran_90 namelist control file

fFile:

Equiulibrium data from PHASE

fFiles:

Directory containing fFiles for interpolation

Index:

List of fFile names and experimental times

*.f:

Actual fFiles

pFile:

Osborne pFile containing profile data

pFiles:

Directory containing pFiles for interpolation

Index:

List of pFile names and experimental times

pFiles:

Actual pFiles

cFile:

Perpendicular diffusivity data from TRANSP

cFiles:

Directory containing cFiles for interpolation

Index:

List of cFile names and experimental times

cFiles:

Actual cFiles

Outputs:

nFile:

File containing neoclassical/profile data for PHASE

nFiles:

Directory containing interpolated nFiles

Index:

List of nFile names and experimental times

*.n:

Actual nFiles

Stage3.nc:

NETCDF file from Stage 3 calculation

Plots:

*.py:

Python scripts to plot Stage 3 data

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 poloidal 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)