Skip to content

physics_variables

Module containing tokamak plasma physics routines

N/A This module contains all the primary plasma physics routines for a tokamak device.

Module containing global variables relating to the plasma physics

iscz = None module-attribute

err242 = None module-attribute

err243 = None module-attribute

f_p_plasma_separatrix_rad = None module-attribute

Separatrix radiation fraction

e_plasma_beta = None module-attribute

[J]

p_plasma_heating_total_mw = None module-attribute

[W]

t_energy_confinement_beta = None module-attribute

[s]

ptarmw = None module-attribute

lambdaio = None module-attribute

drsep = None module-attribute

fio = None module-attribute

fli = None module-attribute

flo = None module-attribute

fui = None module-attribute

fuo = None module-attribute

plimw = None module-attribute

plomw = None module-attribute

puimw = None module-attribute

puomw = None module-attribute

rho_star = None module-attribute

nu_star = None module-attribute

beta_mcdonald = None module-attribute

itart_r = None module-attribute

first_call = None module-attribute

N_CONFINEMENT_SCALINGS = 51 module-attribute

number of energy confinement time scaling laws

m_beam_amu = None module-attribute

beam ion mass (amu)

m_fuel_amu = None module-attribute

average mass of fuel portion of ions (amu)

m_ions_total_amu = None module-attribute

average mass of all ions (amu)

m_plasma_fuel_ions = None module-attribute

Mass of the plasma fuel ions (kg)

m_plasma_ions_total = None module-attribute

Mass of all ions in plasma (kg)

m_plasma_alpha = None module-attribute

Mass of the alpha particles in the plasma (kg)

m_plasma_electron = None module-attribute

Mass of the electrons in the plasma (kg)

m_plasma = None module-attribute

Total mass of the plasma (kg)

alphaj = None module-attribute

current profile index

alphaj_wesson = None module-attribute

Wesson-like current profile index

alphan = None module-attribute

density profile index

alphap = None module-attribute

pressure profile index

fusden_alpha_total = None module-attribute

Alpha particle production rate per unit volume, from plasma and beams [particles/m3/sec]

fusden_plasma_alpha = None module-attribute

Alpha particle production rate per unit volume, just from plasma [particles/m3/sec]

alphat = None module-attribute

temperature profile index

aspect = None module-attribute

aspect ratio (iteration variable 1)

beamfus0 = None module-attribute

multiplier for beam-background fusion calculation

beta_total_vol_avg = None module-attribute

Volume averaged total plasma beta (iteration variable 5) (calculated if stellarator)

beta_fast_alpha = None module-attribute

fast alpha beta component

beta_vol_avg_max = None module-attribute

Max allowable volume averaged beta

beta_vol_avg_min = None module-attribute

Minimum allowable volume averaged beta

beta_beam = None module-attribute

neutral beam beta component

beta_poloidal_vol_avg = None module-attribute

poloidal beta

beta_poloidal_eps = None module-attribute

Poloidal beta and inverse aspcet ratio product

beta_toroidal_vol_avg = None module-attribute

Plasma volume averaged toroidal beta

beta_thermal_toroidal_profile = None module-attribute

toroidal beta profile

beta_thermal_vol_avg = None module-attribute

Plasma volume averaged thermal beta

beta_thermal_poloidal_vol_avg = None module-attribute

Plasma volume averaged poloidal thermal beta

beta_thermal_toroidal_vol_avg = None module-attribute

Plasma volume averaged toloidal thermal beta

beta_norm_total = None module-attribute

normaised total beta

beta_norm_thermal = None module-attribute

normaised thermal beta

beta_norm_toroidal = None module-attribute

normaised toroidal beta

beta_norm_poloidal = None module-attribute

normaised poloidal beta

e_plasma_beta_thermal = None module-attribute

Plasma thermal energy derived from thermal beta

betbm0 = None module-attribute

leading coefficient for NB beta fraction

b_plasma_surface_poloidal_average = None module-attribute

Plasma surface average poloidal field (T)

b_plasma_toroidal_on_axis = None module-attribute

Plasma toroidal field on axis (T) (iteration variable 2)

b_plasma_inboard_toroidal = None module-attribute

Plasma inboard toroidal field (T)

b_plasma_outboard_toroidal = None module-attribute

Plasma outboard toroidal field (T)

b_plasma_toroidal_profile = None module-attribute

toroidal field profile in plasma (T)

b_plasma_total = None module-attribute

Sum of plasma total toroidal + poloidal field (T)

e_plasma_magnetic_stored = None module-attribute

Plasma stored magnetic energy (J)

burnup = None module-attribute

fractional plasma burnup

burnup_in = None module-attribute

fractional plasma burnup user input

b_plasma_vertical_required = None module-attribute

Vertical field needed for plasma equilibrium (T)

c_beta = None module-attribute

Destabalisation parameter for i_beta_norm_max=4 beta limit

csawth = None module-attribute

coeff. for sawteeth effects on burn V-s requirement

f_vol_plasma = None module-attribute

multiplying factor for the plasma volume (normally=1)

f_r_conducting_wall = None module-attribute

maximum ratio of conducting wall distance to plasma minor radius for vertical stability (constraint equation 23)

nd_plasma_electrons_vol_avg = None module-attribute

Plasma volume averaged electron density (/m3) (iteration variable 6)

nd_plasma_fuel_ions_vol_avg = None module-attribute

Plasma volume averaged fuel ion density (/m3)

dlamee = None module-attribute

electron-electron coulomb logarithm

dlamie = None module-attribute

ion-electron coulomb logarithm

nd_plasma_electron_max_array = None module-attribute

Array of plasma electron density upper limits values (/m3)

nd_plasma_alphas_vol_avg = None module-attribute

Plasma volume averaged thermal alpha density (/m3)

nd_beam_ions = None module-attribute

hot beam ion density, variable (/m3)

nd_beam_ions_out = None module-attribute

hot beam ion density from calculation (/m3)

beta_norm_max = None module-attribute

Troyon-like coefficient for beta scaling

beta_norm_max_wesson = None module-attribute

Wesson-like coefficient for beta scaling

beta_norm_max_menard = None module-attribute

Menard-like coefficient for beta scaling

beta_norm_max_original_scaling = None module-attribute

Original scaling coefficient for beta scaling

beta_norm_max_tholerus = None module-attribute

Tholerus-like coefficient for beta scaling

beta_norm_max_stambaugh = None module-attribute

Stambaugh-like coefficient for beta scaling

nd_plasma_electrons_max = None module-attribute

Plasma electron max density limit (/m3)

nd_plasma_ions_total_vol_avg = None module-attribute

Plasma volume averaged total ion density (/m3)

nd_plasma_electron_line = None module-attribute

Plasma line averaged electron density (/m3)

nd_plasma_protons_vol_avg = None module-attribute

Plasma volume averaged proton ash density (/m3)

ntau = None module-attribute

Fusion double product (s/m3)

nTtau = None module-attribute

Lawson triple product [keV s / m3]

nd_plasma_impurities_vol_avg = None module-attribute

Plasma volume averaged impurity (Z > 2) ion density (/m3)

gradient_length_ne = None module-attribute

Max. normalised gradient length in el. density (i_plasma_pedestal==0 only)

gradient_length_te = None module-attribute

Max. normalised gradient length in el. temperature (i_plasma_pedestal==0 only)

beta_poloidal_eps_max = None module-attribute

maximum (eps*beta_poloidal) (constraint equation 6). Note: revised issue #346 "Operation at the tokamak equilibrium poloidal beta-limit in TFTR", 1992 Nucl. Fusion 32 1468

eps = None module-attribute

inverse aspect ratio

f_c_plasma_auxiliary = None module-attribute

fraction of plasma current produced by auxiliary current drive

f_c_plasma_inductive = None module-attribute

fraction of plasma current produced inductively

f_alpha_electron = None module-attribute

fraction of alpha energy to electrons

f_p_alpha_plasma_deposited = None module-attribute

Fraction of alpha power deposited in plasma. Default of 0.95 taken from https://doi.org/10.1088/0029-5515/39/12/305.

f_alpha_ion = None module-attribute

fraction of alpha power to ions

f_plasma_fuel_deuterium = None module-attribute

Plasma deuterium fuel fraction

f_p_div_lower = None module-attribute

fraction of power to the lower divertor in double null configuration (i_single_null = 0 only) (default assumes SN)

ffwal = None module-attribute

factor to convert plasma surface area to first wall area in neutron wall load calculation (i_pflux_fw_neutron=1)

f_nd_plasma_pedestal_greenwald = None module-attribute

fraction of Greenwald density to set as pedestal-top density. If <0, pedestal-top density set manually using nd_plasma_pedestal_electron (i_plasma_pedestal==1). (iteration variable 145)

f_nd_plasma_separatrix_greenwald = None module-attribute

fraction of Greenwald density to set as separatrix density. If <0, separatrix density set manually using nd_plasma_separatrix_electron (i_plasma_pedestal==1). (iteration variable 152)

f_plasma_fuel_helium3 = None module-attribute

Plasma Helium-3 fuel fraction

figmer = None module-attribute

physics figure of merit (= plasma_current*aspect**sbar, where sbar=1)

fkzohm = None module-attribute

Zohm elongation scaling adjustment factor (i_plasma_geometry=2, 3)

f_plasma_fuel_tritium = None module-attribute

Plasma tritium fuel fraction

fusden_total = None module-attribute

fusion reaction rate density, from beams and plasma (reactions/m3/sec)

fusrat_total = None module-attribute

fusion reaction rate, from beams and plasma (reactions/sec)

fusrat_plasma_dt_profile = None module-attribute

Profile of D-T fusion reaction rate in plasma, (reactions/sec)

fusrat_plasma_dd_triton_profile = None module-attribute

Profile of D-D fusion reaction rate (tritium branch) in plasma, (reactions/sec)

fusrat_plasma_dd_helion_profile = None module-attribute

Profile of D-D fusion reaction rate (helium branch) in plasma, (reactions/sec)

fusrat_plasma_dhe3_profile = None module-attribute

Profile of D-3He fusion reaction rate in plasma, (reactions/sec)

fusden_plasma = None module-attribute

fusion reaction rate, just from plasma (reactions/m3/sec)

f_c_plasma_non_inductive = None module-attribute

fraction of the plasma current produced by non-inductive means (iteration variable 44)

ejima_coeff = None module-attribute

Ejima coefficient for resistive startup V-s formula

f_beta_alpha_beam_thermal = None module-attribute

ratio of (fast alpha + neutral beam beta) to thermal beta

hfac = None module-attribute

H factors for an ignited plasma for each energy confinement time scaling law

hfact = None module-attribute

H factor on energy confinement times, radiation corrected (iteration variable 10).

hstar = None module-attribute

H* non-radiation corrected H factor on energy confinement times

t_plasma_energy_confinement_max = None module-attribute

Maximum allowed energy confinement time (s)

i_bootstrap_current = None module-attribute

switch for bootstrap current scaling - =1 ITER 1989 bootstrap scaling (high R/a only) - =2 for Nevins et al general scaling - =3 for Wilson et al numerical scaling - =4 for Sauter et al scaling - =5 for Sakai et al scaling - =6 for ARIES scaling - =7 for Andrade et al scaling - =8 for Hoang et al scaling - =9 for Wong et al scaling - =10 for Gi-I et al scaling - =11 for Gi-II et al scaling - =12 for Sugiyama (L-mode) et al scaling - =13 for Sugiyama (H-mode) et al scaling

i_beta_component = None module-attribute

switch for beta limit scaling (constraint equation 24) - =0 apply limit to total beta - =1 apply limit to thermal beta - =2 apply limit to thermal + neutral beam beta - =3 apply limit to toroidal beta

i_plasma_current = None module-attribute

switch for plasma current scaling to use - =1 Peng analytic fit - =2 Peng double null divertor scaling (ST) - =3 simple ITER scaling (k = 2.2, d = 0.6) - =4 later ITER scaling, a la Uckan - =5 Todd empirical scaling I - =6 Todd empirical scaling II - =7 Connor-Hastie model - =8 Sauter scaling allowing negative triangularity - =9 FIESTA ST fit

i_diamagnetic_current = None module-attribute

switch for diamagnetic current scaling - =0 Do not calculate - =1 Use original TART scaling - =2 Use SCENE scaling

i_density_limit = None module-attribute

switch for density limit to enforce (constraint equation 5) - =1 old ASDEX - =2 Borrass model for ITER (I) - =3 Borrass model for ITER (II) - =4 JET edge radiation - =5 JET simplified - =6 Hugill-Murakami Mq limit - =7 Greenwald limit - =8 ASDEX New

i_beta_fast_alpha = None module-attribute

switch for fast alpha pressure calculation - =0 ITER physics rules (Uckan) fit - =1 Modified fit (D. Ward) - better at high temperature

i_plasma_ignited = None module-attribute

switch for ignition assumption. Obviously, i_plasma_ignited must be zero if current drive is required. If i_plasma_ignited is 1, any auxiliary power is assumed to be used only during plasma start-up, and is excluded from all steady-state power balance calculations. - =0 do not assume plasma ignition - =1 assume ignited (but include auxiliary power in costs)

i_plasma_pedestal = None module-attribute

switch for pedestal profiles: - =0 use original parabolic profiles - =1 use pedestal profile

i_pfirsch_schluter_current = None module-attribute

switch for Pfirsch-Schlüter current scaling (issue #413): - =0 Do not calculate - =1 Use SCENE scaling

nd_plasma_pedestal_electron = None module-attribute

electron density of pedestal [m-3] (`i_plasma_pedestal==1)

nd_plasma_separatrix_electron = None module-attribute

electron density at separatrix [m-3] (`i_plasma_pedestal==1)

alpha_crit = None module-attribute

critical ballooning parameter value

nd_plasma_separatrix_electron_eich_max = None module-attribute

Eich critical electron density at separatrix [m-3]

plasma_res_factor = None module-attribute

plasma resistivity pre-factor

radius_plasma_pedestal_density_norm = None module-attribute

Normalised radius of density pedestal (i_plasma_pedestal==1)

radius_plasma_pedestal_temp_norm = None module-attribute

Normalised radius of temperature pedestal (i_plasma_pedestal==1)

rho_te_max = None module-attribute

r/a where the temperature gradient is largest (i_plasma_pedestal==0)

rho_ne_max = None module-attribute

r/a where the density gradient is largest (i_plasma_pedestal==0)

tbeta = None module-attribute

temperature profile index beta (`i_plasma_pedestal==1)

temp_plasma_pedestal_kev = None module-attribute

Plasma electron temperature of pedestal (keV) (i_plasma_pedestal==1)

temp_plasma_separatrix_kev = None module-attribute

Plasma electron temperature at separatrix (keV) (i_plasma_pedestal==1) calculated if reinke criterion is used (icc=78)

i_beta_norm_max = None module-attribute

Switch for maximum normalised beta scaling:

i_ind_plasma_internal_norm = None module-attribute

Switch for plasma normalised internal inductance scaling:

i_alphaj = None module-attribute

Switch for current profile index scaling:

i_rad_loss = None module-attribute

switch for radiation loss term usage in power balance (see User Guide): - =0 total power lost is scaling power plus radiation - =1 total power lost is scaling power plus core radiation only - =2 total power lost is scaling power only, with no additional allowance for radiation. This is not recommended for power plant models.

i_confinement_time = None module-attribute

switch for energy confinement time scaling law

i_plasma_wall_gap = None module-attribute

Switch for plasma-first wall clearances at the mid-plane: - =0 use 10% of plasma minor radius - =1 use input (dr_fw_plasma_gap_inboard and dr_fw_plasma_gap_outboard)

i_plasma_geometry = None module-attribute

switch for plasma elongation and triangularity calculations: - =0 use input kappa, triang to calculate 95% values - =1 scale q95_min, kappa, triang with aspect ratio (ST) - =2 set kappa to the natural elongation value (Zohm ITER scaling), triang input - =3 set kappa to the natural elongation value (Zohm ITER scaling), triang95 input - =4 use input kappa95, triang95 to calculate separatrix values - =5 use input kappa95, triang95 to calculate separatrix values based on MAST scaling (ST) - =6 use input kappa, triang to calculate 95% values based on MAST scaling (ST) - =7 use input kappa95, triang95 to calculate separatrix values based on fit to FIESTA (ST) - =8 use input kappa, triang to calculate 95% values based on fit to FIESTA (ST) - =9 set kappa to the natural elongation value, triang input - =10 set kappa to maximum stable value at a given aspect ratio (2.6<A<3.6)), triang input (#1399) - =11 set kappa Menard 2016 aspect-ratio-dependent scaling, triang input (#1439)

i_plasma_shape = None module-attribute

switch for plasma boundary shape: - =0 use original PROCESS 2-arcs model - =1 use the Sauter model

itart = None module-attribute

switch for spherical tokamak (ST) models: - =0 use conventional aspect ratio models - =1 use spherical tokamak models

itartpf = None module-attribute

switch for Spherical Tokamak PF models: - =0 use Peng and Strickler (1986) model - =1 use conventional aspect ratio model

i_pflux_fw_neutron = None module-attribute

switch for neutron wall load calculation: - =1 use scaled plasma surface area - =2 use first wall area directly

plasma_square = None module-attribute

plasma squareness used by Sauter plasma shape

kappa = None module-attribute

plasma separatrix elongation (calculated if i_plasma_geometry = 1-5, 7 or 9-10)

kappa95 = None module-attribute

plasma elongation at 95% surface (calculated if i_plasma_geometry = 0-3, 6, or 8-10)

kappa_ipb = None module-attribute

Separatrix elongation calculated for IPB scalings

nd_plasma_electron_on_axis = None module-attribute

central electron density (/m3)

nd_plasma_ions_on_axis = None module-attribute

central ion density (/m3)

m_s_limit = None module-attribute

margin to vertical stability

pres_plasma_thermal_on_axis = None module-attribute

Plasma central thermal pressure (no fast ions or beam pressure) (Pa)

pres_plasma_thermal_total_profile = None module-attribute

Profile of total pressure in plasma (Pa)

pres_plasma_electron_profile = None module-attribute

Profile of electron pressure in plasma (Pa)

pres_plasma_ion_total_profile = None module-attribute

Profile of ion pressure in plasma (Pa)

pres_plasma_fuel_profile = None module-attribute

Profile of fuel pressure in plasma (Pa)

j_plasma_on_axis = None module-attribute

Central plasma current density (A/m2)

j_plasma_bootstrap_sauter_profile = None module-attribute

Profile of bootstrap current density in plasma using Sauter et al scaling (A/m2)

n_plasma_profile_elements = None module-attribute

Number of elements in plasma profile

pres_plasma_thermal_vol_avg = None module-attribute

Volume averaged thermal plasma pressure (no fast ions or beam pressure) (Pa)

f_dd_branching_trit = None module-attribute

branching ratio for DD -> T

pden_plasma_alpha_mw = None module-attribute

Alpha power per volume just from plasma [MW/m3]

pden_alpha_total_mw = None module-attribute

Alpha power per volume from plasma and beams [MW/m3]

f_pden_alpha_electron_mw = None module-attribute

Alpha power per volume to electrons [MW/m3]

p_fw_alpha_mw = None module-attribute

alpha power escaping plasma and reaching first wall (MW)

f_pden_alpha_ions_mw = None module-attribute

alpha power per volume to ions (MW/m3)

p_plasma_alpha_mw = None module-attribute

Alpha power from only the plasma (MW)

p_alpha_total_mw = None module-attribute

Total alpha power from plasma and beams (MW)

p_beam_alpha_mw = None module-attribute

alpha power from hot neutral beam ions (MW)

p_beam_neutron_mw = None module-attribute

neutron power from hot neutral beam ions (MW)

p_beam_dt_mw = None module-attribute

D-T fusion power from hot neutral beam ions (MW)

p_non_alpha_charged_mw = None module-attribute

non-alpha charged particle fusion power (MW)

p_charged_particle_mw = None module-attribute

Total charged particle fusion power [MW]

pden_non_alpha_charged_mw = None module-attribute

Non-alpha charged particle fusion power per volume [MW/m3]

f_temp_plasma_electron_density_vol_avg = None module-attribute

Ratio of density weighted plasma electron tempertaurature to volume averaged (Profile Factor)

p_plasma_inner_rad_mw = None module-attribute

radiation power from inner zone (MW)

pden_plasma_core_rad_mw = None module-attribute

total core radiation power per volume (MW/m3)

p_dd_total_mw = None module-attribute

deuterium-deuterium fusion power (MW)

p_dhe3_total_mw = None module-attribute

deuterium-helium3 fusion power (MW)

p_plasma_separatrix_mw = None module-attribute

power to conducted to the divertor region (MW)

p_plasma_separatrix_rmajor_mw = None module-attribute

Power to conducted to the divertor region per major radius (MW/m)

p_div_bt_q_aspect_rmajor_mw = None module-attribute

EU DEMO divertor protection parameter (MW/T/m)

p_div_lower_separatrix_mw = None module-attribute

Separatrix power conducted to the lower divertor region (calculated if i_single_null = 0) (MW)

p_div_upper_separatrix_mw = None module-attribute

Separatrix power conducted to the upper divertor region (calculated if i_single_null = 0) (MW)

p_div_separatrix_max_mw = None module-attribute

Separatrix power conducted to the divertor with most load (calculated if i_single_null = 0) (MW)

p_dt_total_mw = None module-attribute

Total deuterium-tritium fusion power, from plasma and beams [MW]

p_plasma_dt_mw = None module-attribute

Deuterium-tritium fusion power, just from plasma [MW]

p_plasma_outer_rad_mw = None module-attribute

radiation power from outer zone (MW)

pden_plasma_outer_rad_mw = None module-attribute

edge radiation power per volume (MW/m3)

vs_plasma_internal = None module-attribute

internal plasma V-s

pflux_fw_rad_mw = None module-attribute

Nominal mean radiation load on inside surface of reactor (MW/m2)

pden_ion_electron_equilibration_mw = None module-attribute

ion/electron equilibration power per volume (MW/m3)

plasma_current = None module-attribute

plasma current (A)

c_plasma_peng_analytic = None module-attribute

Peng analytic plasma current (A)

c_plasma_peng_double_null = None module-attribute

Peng double null divertor plasma current (A)

c_plasma_cyclindrical = None module-attribute

Cylindrical plasma current (A)

c_plasma_ipdg89 = None module-attribute

ITER IPDG89 plasma current (A)

c_plasma_todd_empirical_i = None module-attribute

Todd empirical plasma current I (A)

c_plasma_todd_empirical_ii = None module-attribute

Todd empirical plasma current II (A)

c_plasma_connor_hastie = None module-attribute

Connor-Hastie plasma current (A)

c_plasma_sauter = None module-attribute

Sauter plasma current (A)

c_plasma_fiesta_st = None module-attribute

FIESTA ST plasma current (A)

p_plasma_neutron_mw = None module-attribute

Neutron fusion power from just the plasma [MW]

p_neutron_total_mw = None module-attribute

Total neutron fusion power from plasma and beams [MW]

pden_neutron_total_mw = None module-attribute

neutron fusion power per volume from beams and plasma (MW/m3)

pden_plasma_neutron_mw = None module-attribute

neutron fusion power per volume just from plasma (MW/m3)

p_plasma_ohmic_mw = None module-attribute

ohmic heating power (MW)

pden_plasma_ohmic_mw = None module-attribute

ohmic heating power per volume (MW/m3)

p_plasma_loss_mw = None module-attribute

heating power (= transport loss power) (MW) used in confinement time calculation

p_fusion_total_mw = None module-attribute

fusion power (MW)

len_plasma_poloidal = None module-attribute

plasma poloidal perimeter (m)

p_plasma_rad_mw = None module-attribute

total radiation power from inside LCFS (MW)

pden_plasma_rad_mw = None module-attribute

total radiation power per volume (MW/m3)

pradsolmw = None module-attribute

radiation power from SoL (MW)

proton_rate_density = None module-attribute

Proton production rate [particles/m3/sec]

psolradmw = None module-attribute

SOL radiation power (MW) (stellarator only)

pden_plasma_sync_mw = None module-attribute

synchrotron radiation power per volume (MW/m3)

p_plasma_sync_mw = None module-attribute

Total synchrotron radiation power from plasma (MW)

i_l_h_threshold = None module-attribute

switch for L-H mode power threshold scaling to use (see l_h_threshold_powers for list)

p_l_h_threshold_mw = None module-attribute

L-H mode power threshold (MW) (chosen via i_l_h_threshold, and enforced if constraint equation 15 is on)

l_h_threshold_powers = None module-attribute

L-H power threshold for various scalings (MW) - =1 ITER 1996 scaling: nominal - =2 ITER 1996 scaling: upper bound - =3 ITER 1996 scaling: lower bound - =4 ITER 1997 scaling: excluding elongation - =5 ITER 1997 scaling: including elongation - =6 Martin 2008 scaling: nominal - =7 Martin 2008 scaling: 95% upper bound - =8 Martin 2008 scaling: 95% lower bound - =9 Snipes 2000 scaling: nominal - =10 Snipes 2000 scaling: upper bound - =11 Snipes 2000 scaling: lower bound - =12 Snipes 2000 scaling (closed divertor): nominal - =13 Snipes 2000 scaling (closed divertor): upper bound - =14 Snipes 2000 scaling (closed divertor): lower bound - =15 Hubbard et al. 2012 L-I threshold scaling: nominal - =16 Hubbard et al. 2012 L-I threshold scaling: lower bound - =17 Hubbard et al. 2012 L-I threshold scaling: upper bound - =18 Hubbard et al. 2017 L-I threshold scaling - =19 Martin 2008 aspect ratio corrected scaling: nominal - =20 Martin 2008 aspect ratio corrected scaling: 95% upper bound - =21 Martin 2008 aspect ratio corrected scaling: 95% lower bound

p_electron_transport_loss_mw = None module-attribute

electron transport power (MW)

pden_electron_transport_loss_mw = None module-attribute

electron transport power per volume (MW/m3)

p_ion_transport_loss_mw = None module-attribute

ion transport power (MW)

pscalingmw = None module-attribute

Total transport power from scaling law (MW)

pden_ion_transport_loss_mw = None module-attribute

ion transport power per volume (MW/m3)

q0 = None module-attribute

Safety factor on axis

q95 = None module-attribute

Safety factor at 95% flux surface (iteration variable 18) (unless icurr=2 (ST current scaling), in which case q95 = mean edge safety factor qbar)

molflow_plasma_fuelling_required = None module-attribute

plasma fuelling rate (nucleus-pairs/s)

tauratio = None module-attribute

tauratio /1.0/ : ratio of He and pellet particle confinement times

q95_min = None module-attribute

lower limit for edge safety factor

qstar = None module-attribute

cylindrical safety factor

rad_fraction_sol = None module-attribute

SoL radiation fraction

rad_fraction_total = None module-attribute

Radiation fraction total = SoL + LCFS radiation / total power deposited in plasma

f_nd_alpha_electron = None module-attribute

thermal alpha density/electron density (iteration variable 109)

f_nd_protium_electrons = None module-attribute

Seeded f_nd_protium_electrons density / electron density.

ind_plasma_internal_norm = None module-attribute

Plasma normalised internal inductance

ind_plasma_internal_norm_iter_3 = None module-attribute

Plasma normalised internal inductance (ITER type 3)

ind_plasma_internal_norm_wesson = None module-attribute

Wesson-like plasma normalised internal inductance

ind_plasma_internal_menard = None module-attribute

Menard-like plasma normalised internal inductance

ind_plasma = None module-attribute

plasma inductance (H)

rmajor = None module-attribute

plasma major radius (m) (iteration variable 3)

rminor = None module-attribute

plasma minor radius (m)

f_nd_beam_electron = None module-attribute

hot beam density / n_e (iteration variable 7)

f_nd_plasma_carbon_electron = None module-attribute

n_carbon / n_e

rndfuel = None module-attribute

fuel burnup rate (reactions/second)

f_nd_plasma_iron_argon_electron = None module-attribute

n_highZ / n_e

f_nd_plasma_oxygen_electron = None module-attribute

n_oxygen / n_e

f_res_plasma_neo = None module-attribute

neo-classical correction factor to res_plasma

res_plasma = None module-attribute

plasma resistance (ohm)

t_plasma_res_diffusion = None module-attribute

plasma current resistive diffusion time (s)

a_plasma_surface = None module-attribute

plasma surface area

a_plasma_surface_outboard = None module-attribute

outboard plasma surface area

i_single_null = None module-attribute

switch for single null / double null plasma: - =0 for double null - =1 for single null (diverted side down)

f_sync_reflect = None module-attribute

synchrotron wall reflectivity factor

t_electron_energy_confinement = None module-attribute

electron energy confinement time (sec)

tauee_in = None module-attribute

Input electron energy confinement time (sec) (i_confinement_time=48 only)

t_energy_confinement = None module-attribute

global thermal energy confinement time (sec)

t_ion_energy_confinement = None module-attribute

ion energy confinement time (sec)

t_alpha_confinement = None module-attribute

alpha particle confinement time (sec)

f_alpha_energy_confinement = None module-attribute

alpha particle to energy confinement time ratio

temp_plasma_electron_vol_avg_kev = None module-attribute

volume averaged electron temperature (keV) (iteration variable 4)

temp_plasma_electron_on_axis_kev = None module-attribute

central electron temperature (keV)

temp_plasma_electron_density_weighted_kev = None module-attribute

density weighted average electron temperature (keV)

temp_plasma_ion_vol_avg_kev = None module-attribute

volume averaged ion temperature (keV). N.B. calculated from temp_plasma_electron_vol_avg_kev if f_temp_plasma_ion_electron > 0.0

temp_plasma_ion_on_axis_kev = None module-attribute

central ion temperature (keV)

temp_plasma_ion_density_weighted_kev = None module-attribute

density weighted average ion temperature (keV)

f_temp_plasma_ion_electron = None module-attribute

ion temperature / electron temperature(used to calculate temp_plasma_ion_vol_avg_kev if f_temp_plasma_ion_electron > 0.0

triang = None module-attribute

plasma separatrix triangularity (calculated if i_plasma_geometry = 1, 3-5 or 7)

triang95 = None module-attribute

plasma triangularity at 95% surface (calculated if i_plasma_geometry = 0-2, 6, 8 or 9)

vol_plasma = None module-attribute

plasma volume (m3)

vs_plasma_burn_required = None module-attribute

V-s needed during flat-top (heat + burn times) (Wb)

vs_plasma_ramp_required = None module-attribute

V-s needed during ramp-up (Wb)

v_plasma_loop_burn = None module-attribute

Plasma loop voltage during flat-top (V)

vs_plasma_ind_ramp = None module-attribute

Total plasma inductive flux consumption for plasma current ramp-up (Vs)(Wb)

vs_plasma_res_ramp = None module-attribute

Plasma resistive flux consumption for plasma current ramp-up (Vs)(Wb)

vs_plasma_total_required = None module-attribute

total V-s needed (Wb)

pflux_fw_neutron_mw = None module-attribute

average neutron wall load (MW/m2)

pflux_plasma_surface_neutron_avg_mw = None module-attribute

Average neutron flux at plasma surface (MW/m2)

wtgpd = None module-attribute

mass of fuel used per day (g)

a_plasma_poloidal = None module-attribute

plasma poloidal cross-sectional area [m^2]

n_charge_plasma_effective_vol_avg = None module-attribute

Volume averaged plasma effective charge

n_charge_plasma_effective_profile = None module-attribute

Profile of plasma effective charge

n_charge_plasma_effective_mass_weighted_vol_avg = None module-attribute

Plasma mass-weighted volume averaged plasma effective charge

len_plasma_debye_electron_profile = None module-attribute

Profile of electron Debye length in plasma (m)

radius_plasma_deuteron_toroidal_larmor_isotropic_profile = None module-attribute

Profile of deuteron toroidal Larmor radius in plasma, assuming equal speeds in all directions (m)

radius_plasma_deuteron_toroidal_larmor_isotropic_vol_avg = None module-attribute

Volume averaged deuteron toroidal Larmor radius in plasma, assuming equal speeds in all directions (m)

radius_plasma_triton_toroidal_larmor_isotropic_profile = None module-attribute

Profile of triton toroidal Larmor radius in plasma, assuming equal speeds in all directions (m)

radius_plasma_triton_toroidal_larmor_isotropic_vol_avg = None module-attribute

Volume averaged triton toroidal Larmor radius in plasma, assuming equal speeds in all directions (m)

len_plasma_debye_electron_vol_avg = None module-attribute

Volume averaged electron Debye length in plasma (m)

vel_plasma_electron_profile = None module-attribute

Profile of electron thermal velocity in plasma (m/s)

vel_plasma_deuteron_vol_avg = None module-attribute

Volume averaged deuteron thermal velocity in plasma (m/s)

vel_plasma_electron_vol_avg = None module-attribute

Volume averaged electron thermal velocity in plasma (m/s)

vel_plasma_deuteron_profile = None module-attribute

Profile of deuteron thermal velocity in plasma (m/s)

vel_plasma_triton_profile = None module-attribute

Profile of triton thermal velocity in plasma (m/s)

vel_plasma_triton_vol_avg = None module-attribute

Volume averaged triton thermal velocity in plasma (m/s)

vel_plasma_alpha_thermal_profile = None module-attribute

Profile of thermal alpha particle velocity in plasma (m/s)

vel_plasma_alpha_thermal_vol_avg = None module-attribute

Volume averaged thermal alpha particle velocity in plasma (m/s)

vel_plasma_alpha_birth = None module-attribute

Birth velocity of alpha particles in plasma (m/s)

t_plasma_electron_alpha_spitzer_slow_profile = None module-attribute

Profile of electron-alpha Spitzer slowing down time in plasma (s)

t_plasma_electron_alpha_spitzer_slow_vol_avg = None module-attribute

Volume averaged electron-alpha Spitzer slowing down time in plasma (s)

freq_plasma_electron_profile = None module-attribute

Electron plasma frequency profile (Hz)

freq_plasma_electron_vol_avg = None module-attribute

Volume averaged electron plasma frequency (Hz)

freq_plasma_deuteron_profile = None module-attribute

Deuteron plasma frequency profile (Hz)

freq_plasma_larmor_toroidal_electron_profile = None module-attribute

Profile of electron Larmor frequency in plasma due to toroidal magnetic field (Hz)

freq_plasma_larmor_toroidal_electron_vol_avg = None module-attribute

Volume averaged electron Larmor frequency in plasma due to toroidal magnetic field (Hz)

freq_plasma_larmor_toroidal_deuteron_profile = None module-attribute

Profile of deuteron Larmor frequency in plasma due to toroidal magnetic field (Hz)

freq_plasma_larmor_toroidal_deuteron_vol_avg = None module-attribute

Volume averaged deuteron Larmor frequency in plasma due to toroidal magnetic field (Hz)

freq_plasma_larmor_toroidal_triton_profile = None module-attribute

Profile of triton Larmor frequency in plasma due to toroidal magnetic field (Hz)

freq_plasma_larmor_toroidal_triton_vol_avg = None module-attribute

Volume averaged triton Larmor frequency in plasma due to toroidal magnetic field (Hz)

freq_plasma_upper_hybrid_profile = None module-attribute

Profile of upper hybrid frequency in plasma (Hz)

freq_plasma_upper_hybrid_vol_avg = None module-attribute

Volume averaged upper hybrid frequency in plasma (Hz)

t_plasma_electron_electron_collision_profile = None module-attribute

Profile of electron-electron collision time in plasma (s)

t_plasma_electron_electron_collision_vol_avg = None module-attribute

Volume averaged electron-electron collision time in plasma (s)

t_plasma_electron_deuteron_collision_profile = None module-attribute

Profile of electron-deuteron collision time in plasma (s)

t_plasma_electron_deuteron_collision_vol_avg = None module-attribute

Volume averaged electron-deuteron collision time in plasma (s)

t_plasma_electron_triton_collision_profile = None module-attribute

Profile of electron-triton collision time in plasma (s)

t_plasma_electron_triton_collision_vol_avg = None module-attribute

Volume averaged electron-triton collision time in plasma (s)

t_plasma_electron_alpha_thermal_collision_profile = None module-attribute

Profile of electron-alpha collision time in plasma (s)

t_plasma_electron_alpha_thermal_collision_vol_avg = None module-attribute

Volume averaged electron-alpha collision time in plasma (s)

freq_plasma_electron_electron_collision_profile = None module-attribute

Profile of electron-electron collision frequency in plasma (Hz)

freq_plasma_electron_electron_collision_vol_avg = None module-attribute

Volume averaged electron-electron collision frequency in plasma (Hz)

freq_plasma_electron_deuteron_collision_profile = None module-attribute

Profile of electron-deuteron collision frequency in plasma (Hz)

freq_plasma_electron_deuteron_collision_vol_avg = None module-attribute

Volume averaged electron-deuteron collision frequency in plasma (Hz)

freq_plasma_electron_triton_collision_profile = None module-attribute

Profile of electron-triton collision frequency in plasma (Hz)

freq_plasma_electron_triton_collision_vol_avg = None module-attribute

Volume averaged electron-triton collision frequency in plasma (Hz)

freq_plasma_electron_alpha_thermal_collision_profile = None module-attribute

Profile of electron-alpha collision frequency in plasma (Hz)

freq_plasma_electron_alpha_thermal_collision_vol_avg = None module-attribute

Volume averaged electron-alpha collision frequency in plasma (Hz)

len_plasma_electron_electron_mean_free_path_profile = None module-attribute

Profile of electron-electron mean free path in plasma (m)

len_plasma_electron_electron_mean_free_path_vol_avg = None module-attribute

Volume averaged electron-electron mean free path in plasma (m)

len_plasma_electron_deuteron_mean_free_path_profile = None module-attribute

Profile of electron-deuteron mean free path in plasma (m)

len_plasma_electron_deuteron_mean_free_path_vol_avg = None module-attribute

Volume averaged electron-deuteron mean free path in plasma (m)

len_plasma_electron_triton_mean_free_path_profile = None module-attribute

Profile of electron-triton mean free path in plasma (m)

len_plasma_electron_triton_mean_free_path_vol_avg = None module-attribute

Volume averaged electron-triton mean free path in plasma (m)

len_plasma_electron_alpha_thermal_mean_free_path_profile = None module-attribute

Profile of electron-alpha mean free path in plasma (m)

len_plasma_electron_alpha_thermal_mean_free_path_vol_avg = None module-attribute

Volume averaged electron-alpha mean free path in plasma (m)

res_plasma_fuel_spitzer_profile = None module-attribute

Profile of plasma Spitzer resistivity due to fuel ions (ohm m)

res_plasma_fuel_spitzer_vol_avg = None module-attribute

Volume averaged plasma Spitzer resistivity due to fuel ions (ohm m)

init_physics_module()

Initialise the physics module

Source code in process/data_structure/physics_variables.py 
1480
1481
1482
1483
1484
1485
1486
1487
1488
1489
1490
1491
1492
1493
1494
1495
1496
1497
1498
1499
1500
1501
1502
1503
1504
1505
1506
1507
1508
1509
1510
1511
1512
1513
1514
1515
1516
1517
1518
1519
1520
1521
1522
1523
1524
1525
1526
1527
1528
1529
1530
1531
1532
1533
1534
1535
def init_physics_module():
    """Initialise the physics module"""
    global \
        first_call, \
        iscz, \
        err242, \
        err243, \
        f_p_plasma_separatrix_rad, \
        e_plasma_beta, \
        p_plasma_heating_total_mw, \
        t_energy_confinement_beta, \
        ptarmw, \
        lambdaio, \
        drsep, \
        fio, \
        fli, \
        flo, \
        fui, \
        fuo, \
        plimw, \
        plomw, \
        puimw, \
        puomw, \
        rho_star, \
        rho_ne_max, \
        rho_te_max, \
        nu_star, \
        beta_mcdonald, \
        itart_r

    first_call = 1
    iscz = 0
    err242 = 0
    err243 = 0
    f_p_plasma_separatrix_rad = 0.0
    e_plasma_beta = 0.0
    p_plasma_heating_total_mw = 0.0
    t_energy_confinement_beta = 0.0
    ptarmw = 0.0
    lambdaio = 0.0
    drsep = 0.0
    fio = 0.0
    fli = 0.0
    flo = 0.0
    fui = 0.0
    fuo = 0.0
    plimw = 0.0
    plomw = 0.0
    puimw = 0.0
    puomw = 0.0
    rho_ne_max = 0.0
    rho_te_max = 0.0
    rho_star = 0.0
    nu_star = 0.0
    beta_mcdonald = 0.0
    itart_r = 0.0

init_physics_variables()

Source code in process/data_structure/physics_variables.py 
1538
1539
1540
1541
1542
1543
1544
1545
1546
1547
1548
1549
1550
1551
1552
1553
1554
1555
1556
1557
1558
1559
1560
1561
1562
1563
1564
1565
1566
1567
1568
1569
1570
1571
1572
1573
1574
1575
1576
1577
1578
1579
1580
1581
1582
1583
1584
1585
1586
1587
1588
1589
1590
1591
1592
1593
1594
1595
1596
1597
1598
1599
1600
1601
1602
1603
1604
1605
1606
1607
1608
1609
1610
1611
1612
1613
1614
1615
1616
1617
1618
1619
1620
1621
1622
1623
1624
1625
1626
1627
1628
1629
1630
1631
1632
1633
1634
1635
1636
1637
1638
1639
1640
1641
1642
1643
1644
1645
1646
1647
1648
1649
1650
1651
1652
1653
1654
1655
1656
1657
1658
1659
1660
1661
1662
1663
1664
1665
1666
1667
1668
1669
1670
1671
1672
1673
1674
1675
1676
1677
1678
1679
1680
1681
1682
1683
1684
1685
1686
1687
1688
1689
1690
1691
1692
1693
1694
1695
1696
1697
1698
1699
1700
1701
1702
1703
1704
1705
1706
1707
1708
1709
1710
1711
1712
1713
1714
1715
1716
1717
1718
1719
1720
1721
1722
1723
1724
1725
1726
1727
1728
1729
1730
1731
1732
1733
1734
1735
1736
1737
1738
1739
1740
1741
1742
1743
1744
1745
1746
1747
1748
1749
1750
1751
1752
1753
1754
1755
1756
1757
1758
1759
1760
1761
1762
1763
1764
1765
1766
1767
1768
1769
1770
1771
1772
1773
1774
1775
1776
1777
1778
1779
1780
1781
1782
1783
1784
1785
1786
1787
1788
1789
1790
1791
1792
1793
1794
1795
1796
1797
1798
1799
1800
1801
1802
1803
1804
1805
1806
1807
1808
1809
1810
1811
1812
1813
1814
1815
1816
1817
1818
1819
1820
1821
1822
1823
1824
1825
1826
1827
1828
1829
1830
1831
1832
1833
1834
1835
1836
1837
1838
1839
1840
1841
1842
1843
1844
1845
1846
1847
1848
1849
1850
1851
1852
1853
1854
1855
1856
1857
1858
1859
1860
1861
1862
1863
1864
1865
1866
1867
1868
1869
1870
1871
1872
1873
1874
1875
1876
1877
1878
1879
1880
1881
1882
1883
1884
1885
1886
1887
1888
1889
1890
1891
1892
1893
1894
1895
1896
1897
1898
1899
1900
1901
1902
1903
1904
1905
1906
1907
1908
1909
1910
1911
1912
1913
1914
1915
1916
1917
1918
1919
1920
1921
1922
1923
1924
1925
1926
1927
1928
1929
1930
1931
1932
1933
1934
1935
1936
1937
1938
1939
1940
1941
1942
1943
1944
1945
1946
1947
1948
1949
1950
1951
1952
1953
1954
1955
1956
1957
1958
1959
1960
1961
1962
1963
1964
1965
1966
1967
1968
1969
1970
1971
1972
1973
1974
1975
1976
1977
1978
1979
1980
1981
1982
1983
1984
1985
1986
1987
1988
1989
1990
1991
1992
1993
1994
1995
1996
1997
1998
1999
2000
2001
2002
2003
2004
2005
2006
2007
2008
2009
2010
2011
2012
2013
2014
2015
2016
2017
2018
2019
2020
2021
2022
2023
2024
2025
2026
2027
2028
2029
2030
2031
2032
2033
2034
2035
2036
2037
2038
2039
2040
2041
2042
2043
2044
2045
2046
2047
2048
2049
2050
2051
2052
2053
2054
2055
2056
2057
2058
2059
2060
2061
2062
2063
2064
2065
2066
2067
2068
2069
2070
2071
2072
2073
2074
2075
2076
2077
2078
2079
2080
2081
2082
2083
2084
2085
2086
2087
2088
2089
2090
2091
2092
2093
2094
2095
2096
2097
2098
2099
2100
2101
2102
2103
2104
2105
2106
2107
2108
2109
2110
2111
2112
2113
2114
2115
2116
2117
2118
2119
2120
2121
2122
2123
2124
2125
2126
2127
2128
2129
2130
2131
2132
2133
2134
2135
2136
2137
2138
2139
2140
2141
2142
2143
2144
2145
2146
2147
2148
2149
2150
2151
2152
2153
2154
2155
2156
2157
2158
2159
2160
2161
2162
2163
2164
2165
2166
2167
2168
2169
2170
2171
2172
2173
2174
2175
2176
2177
2178
2179
2180
2181
2182
2183
2184
2185
2186
2187
2188
2189
2190
2191
2192
2193
2194
2195
2196
2197
2198
2199
2200
2201
def init_physics_variables():
    global \
        m_beam_amu, \
        m_fuel_amu, \
        m_ions_total_amu, \
        m_plasma_fuel_ions, \
        m_plasma_ions_total, \
        m_plasma_alpha, \
        m_plasma_electron, \
        m_plasma, \
        alphaj, \
        i_alphaj, \
        alphan, \
        alphap, \
        fusden_alpha_total, \
        fusden_plasma_alpha, \
        alphat, \
        aspect, \
        beamfus0, \
        beta_total_vol_avg, \
        beta_fast_alpha, \
        beta_vol_avg_max, \
        beta_vol_avg_min, \
        beta_beam, \
        beta_poloidal_vol_avg, \
        beta_poloidal_eps, \
        beta_toroidal_vol_avg, \
        beta_thermal_toroidal_profile, \
        beta_thermal_vol_avg, \
        beta_thermal_poloidal_vol_avg, \
        beta_thermal_toroidal_vol_avg, \
        beta_norm_total, \
        beta_norm_thermal, \
        beta_norm_poloidal, \
        e_plasma_beta_thermal, \
        beta_norm_toroidal, \
        betbm0, \
        b_plasma_surface_poloidal_average, \
        b_plasma_toroidal_on_axis, \
        b_plasma_toroidal_inboard, \
        b_plasma_toroidal_outboard, \
        b_plasma_toroidal_profile, \
        b_plasma_total, \
        e_plasma_magnetic_stored, \
        burnup, \
        burnup_in, \
        b_plasma_vertical_required, \
        c_beta, \
        csawth, \
        f_vol_plasma, \
        f_r_conducting_wall, \
        nd_plasma_electrons_vol_avg, \
        nd_plasma_fuel_ions_vol_avg, \
        dlamee, \
        dlamie, \
        nd_plasma_electron_max_array, \
        nd_plasma_alphas_vol_avg, \
        nd_beam_ions, \
        nd_beam_ions_out, \
        beta_norm_max, \
        beta_norm_max_wesson, \
        beta_norm_max_menard, \
        beta_norm_max_original_scaling, \
        beta_norm_max_tholerus, \
        beta_norm_max_stambaugh, \
        nd_plasma_electrons_max, \
        nd_plasma_ions_total_vol_avg, \
        nd_plasma_electron_line, \
        nd_plasma_protons_vol_avg, \
        ntau, \
        nTtau, \
        nd_plasma_impurities_vol_avg, \
        beta_poloidal_eps_max, \
        eps, \
        f_c_plasma_auxiliary, \
        f_c_plasma_inductive, \
        f_alpha_electron, \
        f_p_alpha_plasma_deposited, \
        f_alpha_ion, \
        f_plasma_fuel_deuterium, \
        f_p_div_lower, \
        ffwal, \
        f_nd_plasma_pedestal_greenwald, \
        f_nd_plasma_separatrix_greenwald, \
        f_plasma_fuel_helium3, \
        figmer, \
        fkzohm, \
        f_plasma_fuel_tritium, \
        fusden_total, \
        fusrat_total, \
        fusrat_plasma_dt_profile, \
        fusrat_plasma_dd_triton_profile, \
        fusrat_plasma_dd_helion_profile, \
        fusrat_plasma_dhe3_profile, \
        fusden_plasma, \
        f_c_plasma_non_inductive, \
        ejima_coeff, \
        f_beta_alpha_beam_thermal, \
        hfac, \
        hfact, \
        hstar, \
        t_plasma_energy_confinement_max, \
        i_bootstrap_current, \
        i_beta_component, \
        i_plasma_current, \
        i_diamagnetic_current, \
        i_density_limit, \
        i_beta_fast_alpha, \
        i_plasma_ignited, \
        i_plasma_pedestal, \
        i_pfirsch_schluter_current, \
        nd_plasma_pedestal_electron, \
        nd_plasma_separatrix_electron, \
        alpha_crit, \
        nd_plasma_separatrix_electron_eich_max, \
        plasma_res_factor, \
        radius_plasma_pedestal_density_norm, \
        radius_plasma_pedestal_temp_norm, \
        tbeta, \
        temp_plasma_pedestal_kev, \
        temp_plasma_separatrix_kev, \
        i_beta_norm_max, \
        i_rad_loss, \
        i_confinement_time, \
        i_plasma_wall_gap, \
        i_plasma_geometry, \
        i_plasma_shape, \
        itart, \
        itartpf, \
        i_pflux_fw_neutron, \
        plasma_square, \
        kappa, \
        kappa95, \
        kappa_ipb, \
        nd_plasma_electron_on_axis, \
        nd_plasma_ions_on_axis, \
        m_s_limit, \
        pres_plasma_thermal_on_axis, \
        pres_plasma_thermal_total_profile, \
        pres_plasma_electron_profile, \
        pres_plasma_ion_total_profile, \
        pres_plasma_fuel_profile, \
        j_plasma_on_axis, \
        n_plasma_profile_elements, \
        f_dd_branching_trit, \
        pden_plasma_alpha_mw, \
        pden_alpha_total_mw, \
        f_pden_alpha_electron_mw, \
        p_fw_alpha_mw, \
        f_pden_alpha_ions_mw, \
        p_alpha_total_mw, \
        p_plasma_alpha_mw, \
        p_beam_alpha_mw, \
        p_beam_neutron_mw, \
        p_beam_dt_mw, \
        p_non_alpha_charged_mw, \
        pden_non_alpha_charged_mw, \
        f_temp_plasma_electron_density_vol_avg, \
        p_plasma_inner_rad_mw, \
        pden_plasma_core_rad_mw, \
        p_dd_total_mw, \
        p_dhe3_total_mw, \
        p_plasma_separatrix_mw, \
        p_plasma_separatrix_rmajor_mw, \
        p_div_bt_q_aspect_rmajor_mw, \
        p_div_lower_separatrix_mw, \
        p_div_upper_separatrix_mw, \
        p_div_separatrix_max_mw, \
        p_dt_total_mw, \
        p_plasma_dt_mw, \
        p_plasma_outer_rad_mw, \
        pden_plasma_outer_rad_mw, \
        p_charged_particle_mw, \
        vs_plasma_internal, \
        pflux_fw_rad_mw, \
        pden_ion_electron_equilibration_mw, \
        plasma_current, \
        c_plasma_peng_analytic, \
        c_plasma_peng_double_null, \
        c_plasma_cyclindrical, \
        c_plasma_ipdg89, \
        c_plasma_todd_empirical_i, \
        c_plasma_todd_empirical_ii, \
        c_plasma_connor_hastie, \
        c_plasma_sauter, \
        c_plasma_fiesta_st, \
        p_plasma_neutron_mw, \
        p_neutron_total_mw, \
        pden_neutron_total_mw, \
        pden_plasma_neutron_mw, \
        p_plasma_ohmic_mw, \
        pden_plasma_ohmic_mw, \
        p_plasma_loss_mw, \
        p_fusion_total_mw, \
        len_plasma_poloidal, \
        p_plasma_rad_mw, \
        pden_plasma_rad_mw, \
        pradsolmw, \
        proton_rate_density, \
        psolradmw, \
        pden_plasma_sync_mw, \
        p_plasma_sync_mw, \
        i_l_h_threshold, \
        p_l_h_threshold_mw, \
        l_h_threshold_powers, \
        p_electron_transport_loss_mw, \
        pden_electron_transport_loss_mw, \
        p_ion_transport_loss_mw, \
        pscalingmw, \
        pden_ion_transport_loss_mw, \
        q0, \
        q95, \
        molflow_plasma_fuelling_required, \
        tauratio, \
        q95_min, \
        qstar, \
        rad_fraction_sol, \
        rad_fraction_total, \
        f_nd_alpha_electron, \
        f_nd_protium_electrons, \
        ind_plasma_internal_norm, \
        ind_plasma_internal_norm_wesson, \
        ind_plasma_internal_norm_menard, \
        ind_plasma_internal_norm_iter_3, \
        i_ind_plasma_internal_norm, \
        ind_plasma, \
        rmajor, \
        rminor, \
        f_nd_beam_electron, \
        f_nd_plasma_carbon_electron, \
        rndfuel, \
        f_nd_plasma_iron_argon_electron, \
        f_nd_plasma_oxygen_electron, \
        f_res_plasma_neo, \
        res_plasma, \
        rho_plasma_spitzer_classical_profile, \
        t_plasma_res_diffusion, \
        a_plasma_surface, \
        a_plasma_surface_outboard, \
        i_single_null, \
        f_sync_reflect, \
        t_electron_energy_confinement, \
        tauee_in, \
        t_energy_confinement, \
        t_ion_energy_confinement, \
        t_alpha_confinement, \
        f_alpha_energy_confinement, \
        temp_plasma_electron_vol_avg_kev, \
        temp_plasma_electron_on_axis_kev, \
        temp_plasma_electron_density_weighted_kev, \
        temp_plasma_ion_vol_avg_kev, \
        temp_plasma_ion_on_axis_kev, \
        temp_plasma_ion_density_weighted_kev, \
        f_temp_plasma_ion_electron, \
        triang, \
        triang95, \
        vol_plasma, \
        vs_plasma_burn_required, \
        vs_plasma_ramp_required, \
        v_plasma_loop_burn, \
        vs_plasma_ind_ramp, \
        vs_plasma_res_ramp, \
        vs_plasma_total_required, \
        pflux_fw_neutron_mw, \
        pflux_plasma_surface_neutron_avg_mw, \
        wtgpd, \
        a_plasma_poloidal, \
        n_charge_plasma_effective_vol_avg, \
        n_charge_plasma_effective_profile, \
        n_charge_plasma_effective_mass_weighted_vol_avg, \
        j_plasma_bootstrap_sauter_profile, \
        len_plasma_debye_electron_profile, \
        radius_plasma_deuteron_toroidal_larmor_isotropic_profile, \
        radius_plasma_deuteron_toroidal_larmor_isotropic_vol_avg, \
        radius_plasma_triton_toroidal_larmor_isotropic_profile, \
        radius_plasma_triton_toroidal_larmor_isotropic_vol_avg, \
        len_plasma_debye_electron_vol_avg, \
        vel_plasma_electron_profile, \
        vel_plasma_deuteron_vol_avg, \
        vel_plasma_electron_vol_avg, \
        vel_plasma_deuteron_profile, \
        vel_plasma_triton_profile, \
        vel_plasma_triton_vol_avg, \
        vel_plasma_alpha_thermal_profile, \
        vel_plasma_alpha_thermal_vol_avg, \
        vel_plasma_alpha_birth, \
        plasma_coulomb_log_electron_electron_profile, \
        plasma_coulomb_log_electron_electron_vol_avg, \
        plasma_coulomb_log_electron_deuteron_profile, \
        plasma_coulomb_log_electron_deuteron_vol_avg, \
        plasma_coulomb_log_electron_triton_profile, \
        plasma_coulomb_log_electron_triton_vol_avg, \
        plasma_coulomb_log_deuteron_triton_profile, \
        plasma_coulomb_log_deuteron_triton_vol_avg, \
        plasma_coulomb_log_electron_alpha_thermal_profile, \
        plasma_coulomb_log_electron_alpha_thermal_vol_avg, \
        t_plasma_electron_alpha_spitzer_slow_profile, \
        t_plasma_electron_alpha_spitzer_slow_vol_avg, \
        freq_plasma_electron_profile, \
        freq_plasma_electron_vol_avg, \
        freq_plasma_deuteron_profile, \
        freq_plasma_larmor_toroidal_electron_profile, \
        freq_plasma_larmor_toroidal_deuteron_profile, \
        freq_plasma_larmor_toroidal_triton_profile, \
        freq_plasma_upper_hybrid_profile, \
        freq_plasma_upper_hybrid_vol_avg, \
        t_plasma_electron_electron_collision_profile, \
        t_plasma_electron_electron_collision_vol_avg, \
        t_plasma_electron_deuteron_collision_profile, \
        t_plasma_electron_deuteron_collision_vol_avg, \
        t_plasma_electron_triton_collision_profile, \
        t_plasma_electron_triton_collision_vol_avg, \
        t_plasma_electron_alpha_thermal_collision_profile, \
        t_plasma_electron_alpha_thermal_collision_vol_avg, \
        freq_plasma_electron_electron_collision_profile, \
        freq_plasma_electron_electron_collision_vol_avg, \
        freq_plasma_electron_deuteron_collision_profile, \
        freq_plasma_electron_deuteron_collision_vol_avg, \
        freq_plasma_electron_triton_collision_profile, \
        freq_plasma_electron_triton_collision_vol_avg, \
        freq_plasma_electron_alpha_thermal_collision_profile, \
        freq_plasma_electron_alpha_thermal_collision_vol_avg, \
        len_plasma_electron_electron_mean_free_path_profile, \
        len_plasma_electron_electron_mean_free_path_vol_avg, \
        len_plasma_electron_deuteron_mean_free_path_profile, \
        len_plasma_electron_deuteron_mean_free_path_vol_avg, \
        len_plasma_electron_triton_mean_free_path_profile, \
        len_plasma_electron_triton_mean_free_path_vol_avg, \
        len_plasma_electron_alpha_thermal_mean_free_path_profile, \
        len_plasma_electron_alpha_thermal_mean_free_path_vol_avg, \
        res_plasma_fuel_spitzer_profile, \
        res_plasma_fuel_spitzer_vol_avg

    m_beam_amu = 0.0
    m_fuel_amu = 0.0
    m_ions_total_amu = 0.0
    m_plasma_fuel_ions = 0.0
    m_plasma_ions_total = 0.0
    m_plasma_alpha = 0.0
    m_plasma_electron = 0.0
    m_plasma = 0.0
    alphaj = 1.0
    i_alphaj = 0
    alphan = 0.25
    alphap = 0.0
    fusden_alpha_total = 0.0
    fusden_plasma_alpha = 0.0
    alphat = 0.5
    aspect = 2.907
    beamfus0 = 1.0
    beta_total_vol_avg = 0.042
    beta_fast_alpha = 0.0
    beta_vol_avg_max = 0.0
    beta_vol_avg_min = 0.0
    beta_beam = 0.0
    beta_poloidal_vol_avg = 0.0
    beta_poloidal_eps = 0.0
    beta_toroidal_vol_avg = 0.0
    beta_thermal_toroidal_profile = []
    beta_thermal_vol_avg = 0.0
    beta_thermal_poloidal_vol_avg = 0.0
    beta_thermal_toroidal_vol_avg = 0.0
    beta_norm_total = 0.0
    beta_norm_thermal = 0.0
    beta_norm_poloidal = 0.0
    e_plasma_beta_thermal = 0.0
    beta_norm_toroidal = 0.0
    betbm0 = 1.5
    b_plasma_surface_poloidal_average = 0.0
    b_plasma_toroidal_on_axis = 5.68
    b_plasma_toroidal_inboard = 0.0
    b_plasma_toroidal_outboard = 0.0
    b_plasma_toroidal_profile = []
    b_plasma_total = 0.0
    e_plasma_magnetic_stored = 0.0
    burnup = 0.0
    burnup_in = 0.0
    b_plasma_vertical_required = 0.0
    c_beta = 0.5
    csawth = 1.0
    f_vol_plasma = 1.0
    f_r_conducting_wall = 1.35
    nd_plasma_electrons_vol_avg = 9.8e19
    nd_plasma_fuel_ions_vol_avg = 0.0
    dlamee = 0.0
    dlamie = 0.0
    nd_plasma_electron_max_array = np.zeros(8, dtype=np.float64)
    nd_plasma_alphas_vol_avg = 0.0
    nd_beam_ions = 0.0
    nd_beam_ions_out = 0.0
    beta_norm_max = 3.5
    beta_norm_max_wesson = 0.0
    beta_norm_max_menard = 0.0
    beta_norm_max_original_scaling = 0.0
    beta_norm_max_tholerus = 0.0
    beta_norm_max_stambaugh = 0.0
    nd_plasma_electrons_max = 0.0
    nd_plasma_ions_total_vol_avg = 0.0
    nd_plasma_electron_line = 0.0
    nd_plasma_protons_vol_avg = 0.0
    ntau = 0.0
    nTtau = 0.0
    nd_plasma_impurities_vol_avg = 0.0
    beta_poloidal_eps_max = 1.38
    eps = 0.34399724802
    f_c_plasma_auxiliary = 0.0
    f_c_plasma_inductive = 0.0
    f_alpha_electron = 0.0
    f_p_alpha_plasma_deposited = 0.95
    f_alpha_ion = 0.0
    f_plasma_fuel_deuterium = 0.5
    f_p_div_lower = 1.0
    ffwal = 0.92
    f_nd_plasma_pedestal_greenwald = 0.85
    f_nd_plasma_separatrix_greenwald = 0.50
    f_plasma_fuel_helium3 = 0.0
    figmer = 0.0
    fkzohm = 1.0
    f_plasma_fuel_tritium = 0.5
    fusden_total = 0.0
    fusrat_total = 0.0
    fusrat_plasma_dt_profile = []
    fusrat_plasma_dd_triton_profile = []
    fusrat_plasma_dd_helion_profile = []
    fusrat_plasma_dhe3_profile = []
    fusden_plasma = 0.0
    f_c_plasma_non_inductive = 1.0
    ejima_coeff = 0.4
    f_beta_alpha_beam_thermal = 0.0
    hfac = np.zeros(N_CONFINEMENT_SCALINGS, dtype=np.float64)
    hfact = 1.0
    hstar = 1.0
    t_plasma_energy_confinement_max = 10.0
    i_bootstrap_current = 3
    i_beta_component = 0
    i_plasma_current = 4
    i_diamagnetic_current = 0
    i_density_limit = 8
    i_beta_fast_alpha = 1
    i_plasma_ignited = 0
    i_plasma_pedestal = 1
    i_pfirsch_schluter_current = 0
    nd_plasma_pedestal_electron = 4.0e19
    nd_plasma_separatrix_electron = 3.0e19
    alpha_crit = 0.0
    nd_plasma_separatrix_electron_eich_max = 0.0
    plasma_res_factor = 1.0
    radius_plasma_pedestal_density_norm = 1.0
    radius_plasma_pedestal_temp_norm = 1.0
    tbeta = 2.0
    temp_plasma_pedestal_kev = 1.0
    temp_plasma_separatrix_kev = 0.1
    i_beta_norm_max = 1
    i_rad_loss = 1
    i_confinement_time = 34
    i_plasma_wall_gap = 1
    i_plasma_geometry = 0
    i_plasma_shape = 0
    itart = 0
    itartpf = 0
    i_pflux_fw_neutron = 1
    plasma_square = 0.0
    kappa = 1.792
    kappa95 = 1.6
    kappa_ipb = 0.0
    nd_plasma_electron_on_axis = 0.0
    nd_plasma_ions_on_axis = 0.0
    m_s_limit = 0.3
    pres_plasma_thermal_on_axis = 0.0
    pres_plasma_thermal_total_profile = []
    pres_plasma_electron_profile = []
    pres_plasma_ion_total_profile = []
    pres_plasma_fuel_profile = []
    j_plasma_on_axis = 0.0
    j_plasma_bootstrap_sauter_profile = []
    n_plasma_profile_elements = 500
    f_dd_branching_trit = 0.0
    pden_plasma_alpha_mw = 0.0
    pden_alpha_total_mw = 0.0
    f_pden_alpha_electron_mw = 0.0
    p_fw_alpha_mw = 0.0
    f_pden_alpha_ions_mw = 0.0
    p_alpha_total_mw = 0.0
    p_plasma_alpha_mw = 0.0
    p_beam_alpha_mw = 0.0
    p_beam_neutron_mw = 0.0
    p_beam_dt_mw = 0.0
    p_non_alpha_charged_mw = 0.0
    pden_non_alpha_charged_mw = 0.0
    f_temp_plasma_electron_density_vol_avg = 0.0
    p_plasma_inner_rad_mw = 0.0
    pden_plasma_core_rad_mw = 0.0
    p_dd_total_mw = 0.0
    p_dhe3_total_mw = 0.0
    p_plasma_separatrix_mw = 0.0
    p_plasma_separatrix_rmajor_mw = 0.0
    p_div_bt_q_aspect_rmajor_mw = 0.0
    p_div_lower_separatrix_mw = 0.0
    p_div_lower_separatrix_mw = 0.0
    p_div_upper_separatrix_mw = 0.0
    p_div_separatrix_max_mw = 0.0
    p_dt_total_mw = 0.0
    p_plasma_dt_mw = 0.0
    p_plasma_outer_rad_mw = 0.0
    pden_plasma_outer_rad_mw = 0.0
    p_charged_particle_mw = 0.0
    vs_plasma_internal = 0.0
    pflux_fw_rad_mw = 0.0
    pden_ion_electron_equilibration_mw = 0.0
    plasma_current = 0.0
    c_plasma_peng_analytic = 0.0
    c_plasma_peng_double_null = 0.0
    c_plasma_cyclindrical = 0.0
    c_plasma_ipdg89 = 0.0
    c_plasma_todd_empirical_i = 0.0
    c_plasma_todd_empirical_ii = 0.0
    c_plasma_connor_hastie = 0.0
    c_plasma_sauter = 0.0
    c_plasma_fiesta_st = 0.0
    p_plasma_neutron_mw = 0.0
    p_neutron_total_mw = 0.0
    pden_neutron_total_mw = 0.0
    pden_plasma_neutron_mw = 0.0
    p_plasma_ohmic_mw = 0.0
    pden_plasma_ohmic_mw = 0.0
    p_plasma_loss_mw = 0.0
    p_fusion_total_mw = 0.0
    len_plasma_poloidal = 0.0
    p_plasma_rad_mw = 0.0
    pden_plasma_rad_mw = 0.0
    pradsolmw = 0.0
    proton_rate_density = 0.0
    psolradmw = 0.0
    pden_plasma_sync_mw = 0.0
    p_plasma_sync_mw = 0.0
    i_l_h_threshold = 19
    p_l_h_threshold_mw = 0.0
    l_h_threshold_powers = np.zeros(21, dtype=np.float64)
    p_electron_transport_loss_mw = 0.0
    pden_electron_transport_loss_mw = 0.0
    p_ion_transport_loss_mw = 0.0
    pscalingmw = 0.0
    pden_ion_transport_loss_mw = 0.0
    q0 = 1.0
    q95 = 0.0
    molflow_plasma_fuelling_required = 0.0
    tauratio = 1.0
    q95_min = 0.0
    qstar = 0.0
    rad_fraction_sol = 0.8
    rad_fraction_total = 0.0
    f_nd_alpha_electron = 0.1
    f_nd_protium_electrons = 0.0
    ind_plasma_internal_norm = 0.9
    ind_plasma_internal_norm_wesson = 0.0
    ind_plasma_internal_norm_menard = 0.0
    ind_plasma_internal_norm_iter_3 = 0.0
    i_ind_plasma_internal_norm = 0
    ind_plasma = 0.0
    rmajor = 8.14
    rminor = 0.0
    f_nd_beam_electron = 0.005
    f_nd_plasma_carbon_electron = 0.0
    rndfuel = 0.0
    f_nd_plasma_iron_argon_electron = 0.0
    f_nd_plasma_oxygen_electron = 0.0
    f_res_plasma_neo = 0.0
    res_plasma = 0.0
    rho_plasma_spitzer_classical_profile = []
    t_plasma_res_diffusion = 0.0
    a_plasma_surface = 0.0
    a_plasma_surface_outboard = 0.0
    i_single_null = 1
    f_sync_reflect = 0.6
    t_electron_energy_confinement = 0.0
    tauee_in = 0.0
    t_energy_confinement = 0.0
    t_ion_energy_confinement = 0.0
    t_alpha_confinement = 0.0
    f_alpha_energy_confinement = 0.0
    temp_plasma_electron_vol_avg_kev = 12.9
    temp_plasma_electron_on_axis_kev = 0.0
    temp_plasma_electron_density_weighted_kev = 0.0
    temp_plasma_ion_vol_avg_kev = 12.9
    temp_plasma_ion_on_axis_kev = 0.0
    temp_plasma_ion_density_weighted_kev = 0.0
    f_temp_plasma_ion_electron = 1.0
    triang = 0.36
    triang95 = 0.24
    vol_plasma = 0.0
    vs_plasma_burn_required = 0.0
    vs_plasma_ramp_required = 0.0
    v_plasma_loop_burn = 0.0
    vs_plasma_ind_ramp = 0.0
    vs_plasma_res_ramp = 0.0
    vs_plasma_total_required = 0.0
    pflux_fw_neutron_mw = 0.0
    pflux_plasma_surface_neutron_avg_mw = 0.0
    wtgpd = 0.0
    a_plasma_poloidal = 0.0
    n_charge_plasma_effective_vol_avg = 0.0
    n_charge_plasma_effective_profile = []
    n_charge_plasma_effective_mass_weighted_vol_avg = 0.0
    len_plasma_debye_electron_profile = []
    radius_plasma_deuteron_toroidal_larmor_isotropic_profile = []
    radius_plasma_deuteron_toroidal_larmor_isotropic_vol_avg = 0.0
    radius_plasma_triton_toroidal_larmor_isotropic_profile = []
    radius_plasma_triton_toroidal_larmor_isotropic_vol_avg = 0.0
    len_plasma_debye_electron_vol_avg = 0.0
    vel_plasma_electron_profile = []
    vel_plasma_electron_vol_avg = 0.0
    vel_plasma_deuteron_profile = []
    vel_plasma_deuteron_vol_avg = 0.0
    vel_plasma_triton_profile = []
    vel_plasma_triton_vol_avg = 0.0
    vel_plasma_alpha_thermal_profile = []
    vel_plasma_alpha_thermal_vol_avg = 0.0
    vel_plasma_alpha_birth = 0.0
    plasma_coulomb_log_electron_electron_profile = []
    plasma_coulomb_log_electron_electron_vol_avg = 0.0
    plasma_coulomb_log_electron_deuteron_profile = []
    plasma_coulomb_log_electron_deuteron_vol_avg = 0.0
    plasma_coulomb_log_electron_triton_profile = []
    plasma_coulomb_log_electron_triton_vol_avg = 0.0
    plasma_coulomb_log_deuteron_triton_profile = []
    plasma_coulomb_log_deuteron_triton_vol_avg = 0.0
    plasma_coulomb_log_electron_alpha_thermal_profile = []
    plasma_coulomb_log_electron_alpha_thermal_vol_avg = 0.0
    freq_plasma_electron_profile = []
    freq_plasma_electron_vol_avg = 0.0
    freq_plasma_deuteron_profile = []
    freq_plasma_larmor_toroidal_electron_profile = []
    freq_plasma_larmor_toroidal_deuteron_profile = []
    freq_plasma_larmor_toroidal_triton_profile = []
    freq_plasma_upper_hybrid_profile = []
    freq_plasma_upper_hybrid_vol_avg = 0.0
    t_plasma_electron_electron_collision_profile = []
    t_plasma_electron_electron_collision_vol_avg = 0.0
    t_plasma_electron_deuteron_collision_profile = []
    t_plasma_electron_deuteron_collision_vol_avg = 0.0
    t_plasma_electron_triton_collision_profile = []
    t_plasma_electron_triton_collision_vol_avg = 0.0
    t_plasma_electron_alpha_thermal_collision_profile = []
    t_plasma_electron_alpha_thermal_collision_vol_avg = 0.0
    t_plasma_electron_alpha_spitzer_slow_profile = []
    t_plasma_electron_alpha_spitzer_slow_vol_avg = 0.0
    freq_plasma_electron_electron_collision_profile = []
    freq_plasma_electron_electron_collision_vol_avg = 0.0
    freq_plasma_electron_deuteron_collision_profile = []
    freq_plasma_electron_deuteron_collision_vol_avg = 0.0
    freq_plasma_electron_triton_collision_profile = []
    freq_plasma_electron_triton_collision_vol_avg = 0.0
    freq_plasma_electron_alpha_thermal_collision_profile = []
    freq_plasma_electron_alpha_thermal_collision_vol_avg = 0.0
    len_plasma_electron_electron_mean_free_path_profile = []
    len_plasma_electron_electron_mean_free_path_vol_avg = 0.0
    len_plasma_electron_deuteron_mean_free_path_profile = []
    len_plasma_electron_deuteron_mean_free_path_vol_avg = 0.0
    len_plasma_electron_triton_mean_free_path_profile = []
    len_plasma_electron_triton_mean_free_path_vol_avg = 0.0
    len_plasma_electron_alpha_thermal_mean_free_path_profile = []
    len_plasma_electron_alpha_thermal_mean_free_path_vol_avg = 0.0
    res_plasma_fuel_spitzer_profile = []
    res_plasma_fuel_spitzer_vol_avg = 0.0