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

LABELS_CONFINEMENT_SCALINGS = ['User input electron confinement ', 'Neo-Alcator (Ohmic)', 'Mirnov (H)', 'Merezkhin-Muhkovatov (Ohmic)(L)', 'Shimomura (H)', 'Kaye-Goldston (L)', 'ITER 89-P (L)', 'ITER 89-O (L)', 'Rebut-Lallia (L)', 'Goldston (L)', 'T10 (L)', 'JAERI / Odajima-Shimomura (L)', 'Kaye-Big Complex (L)', 'ITER H90-P (H)', 'ITER 89-P & 89-O min (L)', 'Riedel (L)', 'Christiansen (L)', 'Lackner-Gottardi (L)', 'Neo-Kaye (L)', 'Riedel (H)', 'ITER H90-P amended (H)', 'LHD (Stell)', 'Gyro-reduced Bohm (Stell)', 'Lackner-Gottardi (Stell)', 'ITER-93H ELM-free (H)', 'TITAN RFP OBSOLETE ', 'ITER H-97P ELM-free (H)', 'ITER H-97P ELMy (H)', 'ITER-96P (ITER-97L) (L)', 'Valovic modified ELMy (H)', 'Kaye 98 modified (L)', 'ITERH-PB98P(y) (H)', 'IPB98(y) (H)', 'IPB98(y,1) (H)', 'IPB98(y,2) (H)', 'IPB98(y,3) (H)', 'IPB98(y,4) (H)', 'ISS95 (Stell)', 'ISS04 (Stell)', 'DS03 beta-independent (H)', 'Murari "Non-power law" (H)', 'Petty 2008 (ST)(H)', 'Lang high density (H)', 'Hubbard 2017 - nominal (I)', 'Hubbard 2017 - lower (I)', 'Hubbard 2017 - upper (I)', 'Menard NSTX (ST)(H)', 'Menard NSTX-Petty08 hybrid (ST)(H)', 'Buxton NSTX gyro-Bohm (ST)(H)', 'ITPA20 (H)', 'ITPA20-IL (H)'] module-attribute

labels describing energy confinement 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_poloidal_average = None module-attribute

Plasma 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 (see description in LABELS_CONFINEMENT_SCALINGS)

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

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)

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

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

freq_plasma_electron_profile = None module-attribute

Electron 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_deuteron_profile = None module-attribute

Profile of 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)

init_physics_module()

Initialise the physics module

Source code in process/data_structure/physics_variables.py

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

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_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_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, \
        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, \
        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, \
        plfux_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, \
        len_plasma_debye_electron_vol_avg, \
        vel_plasma_electron_profile, \
        vel_plasma_deuteron_profile, \
        vel_plasma_triton_profile, \
        plasma_coulomb_log_electron_electron_profile, \
        plasma_coulomb_log_electron_deuteron_profile, \
        plasma_coulomb_log_electron_triton_profile, \
        freq_plasma_electron_profile, \
        freq_plasma_larmor_toroidal_electron_profile, \
        freq_plasma_larmor_toroidal_deuteron_profile, \
        freq_plasma_larmor_toroidal_triton_profile

    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_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_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
    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
    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
    plfux_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 = []
    len_plasma_debye_electron_vol_avg = 0.0
    vel_plasma_electron_profile = []
    vel_plasma_deuteron_profile = []
    vel_plasma_triton_profile = []
    plasma_coulomb_log_electron_electron_profile = []
    plasma_coulomb_log_electron_deuteron_profile = []
    plasma_coulomb_log_electron_triton_profile = []
    freq_plasma_electron_profile = []
    freq_plasma_larmor_toroidal_electron_profile = []
    freq_plasma_larmor_toroidal_deuteron_profile = []
    freq_plasma_larmor_toroidal_triton_profile = []