PROCESS Variable Descriptions
Introduction
Variables marked with an * are private variables and cannot be accessed outside of their Fortran module scope.
The PROCESS convention on inputs dictates that module variables which can be set from the
input file should be initialised in a routine called init_mod
where mod
is replaced with the
Fortran module name. However, some variables can also be initialised here too and will register as inputs
when they are not.
Output types signify that these module variables are set by models. This does necessarily mean an "output" will appear in the MFile/Outfile.
blanket_library
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
volshldi | Output | real | - | Volume of inboard and outboard shield (m3) |
volshldo | Output | real | - | Volume of inboard and outboard shield (m3) |
volvvi | Output | real | - | Volume of inboard and outboard Vacuum Vessel (m3) |
volvvo | Output | real | - | Volume of inboard and outboard Vacuum Vessel (m3) |
hcryopf | Output | real | - | Clearance between uppermost PF coil and cryostat lid (m) |
vfblkti | Output | real | - | Inboard/outboard void fraction of blanket |
vfblkto | Output | real | - | Inboard/outboard void fraction of blanket |
bldepti | Output | real | - | Inboard/outboard blanket coolant channel length (radial direction) (m) |
bldepto | Output | real | - | Inboard/outboard blanket coolant channel length (radial direction) (m) |
blwidti | Output | real | - | Inboard/outboard blanket mid-plan toroidal circumference for segment (m) |
blwidto | Output | real | - | Inboard/outboard blanket mid-plan toroidal circumference for segment (m) |
bllengi | Output | real | - | Inboard/outboard blanket segment poloidal length (m) |
bllengo | Output | real | - | Inboard/outboard blanket segment poloidal length (m) |
bzfllengi | Output | real | - | Inboard/outboard primary blanket flow lengths (m) |
bzfllengo | Output | real | - | Inboard/outboard primary blanket flow lengths (m) |
bzfllengi_liq | Output | real | - | Inboard/outboard secondary blanket flow lengths (m) |
bzfllengo_liq | Output | real | - | Inboard/outboard secondary blanket flow lengths (m) |
pnucfwi | Output | real | - | Inboard/outboard first wall nuclear heating (MW) |
pnucfwo | Output | real | - | Inboard/outboard first wall nuclear heating (MW) |
tpeakfwi | Output | real | - | Inboard/outboard first wall peak temperature (K) |
tpeakfwo | Output | real | - | Inboard/outboard first wall peak temperature (K) |
mffwi | Output | real | - | Inboard/outboard total mass flow rate to remove inboard FW power (kg/s) |
mffwo | Output | real | - | Inboard/outboard total mass flow rate to remove inboard FW power (kg/s) |
mffw | Output | real | - | Inboard/outboard total mass flow rate to remove inboard FW power (kg/s) |
npfwi | Output | real | - | Inboard/utboard total number of pipes |
npfwo | Output | real | - | Inboard/utboard total number of pipes |
mffwpi | Output | real | - | Inboard/outboard mass flow rate per coolant pipe (kg/s) |
mffwpo | Output | real | - | Inboard/outboard mass flow rate per coolant pipe (kg/s) |
pnucblkti | Output | real | - | Neutron power deposited inboard/outboard blanket blanket (MW) |
pnucblkto | Output | real | - | Neutron power deposited inboard/outboard blanket blanket (MW) |
mfblkti | Output | real | - | Inboard/outboard blanket mass flow rate for coolant (kg/s) |
mfblkto | Output | real | - | Inboard/outboard blanket mass flow rate for coolant (kg/s) |
mfblkt | Output | real | - | Inboard/outboard blanket mass flow rate for coolant (kg/s) |
mfblkti_liq | Output | real | - | Inboard/outboard blanket mass flow rate for liquid breeder (kg/s) |
mfblkto_liq | Output | real | - | Inboard/outboard blanket mass flow rate for liquid breeder (kg/s) |
mfblkt_liq | Output | real | - | Inboard/outboard blanket mass flow rate for liquid breeder (kg/s) |
mftotal | Output | real | - | Total mass flow rate for coolant (kg/s) |
npblkti | Output | real | - | Inboard/outboard total num of pipes |
npblkto | Output | real | - | Inboard/outboard total num of pipes |
mfblktpi | Output | real | - | Inboard/outboard mass flow rate per coolant pipe (kg/s) |
mfblktpo | Output | real | - | Inboard/outboard mass flow rate per coolant pipe (kg/s) |
velblkti | Output | real | - | Inboard/outboard coolant velocity in blanket (m/s) |
velblkto | Output | real | - | Inboard/outboard coolant velocity in blanket (m/s) |
htpmw_fwi | Output | real | - | Inboard/outboard first wall pumping power (MW) |
htpmw_fwo | Output | real | - | Inboard/outboard first wall pumping power (MW) |
htpmw_blkti | Output | real | - | Inboard/outboard blanket pumping power (MW) |
htpmw_blkto | Output | real | - | Inboard/outboard blanket pumping power (MW) |
htpmw_fw_blkti | Output | real | - | Inboard/outboard fw and blanket pumping power (MW) |
htpmw_fw_blkto | Output | real | - | Inboard/outboard fw and blanket pumping power (MW) |
hblnkt | Output | real | - | Blanket internal half-height (m) |
hshld | Output | real | - | Shield internal half-height (m) |
hvv | Output | real | - | Vacuum vessel internal half-height (m) |
icomponent | Output | integer | - | Switch used to specify selected component: blanket=0, shield=1, vacuum vessel=2 |
maths_library
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
double | Parameter | integer | 8 |
physics_module
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
iscz | Output | integer | - | |
err242 | Output | integer | - | |
err243 | Output | integer | - | |
rad_fraction_lcfs | Output | real | - | |
total_plasma_internal_energy | Output | real | - | |
total_loss_power | Output | real | - | |
total_energy_conf_time | Output | real | - | |
ptarmw | Output | real | - | |
lambdaio | Output | real | - | |
drsep | Output | real | - | |
fio | Output | real | - | |
fli | Output | real | - | |
flo | Output | real | - | |
fui | Output | real | - | |
fuo | Output | real | - | |
plimw | Output | real | - | |
plomw | Output | real | - | |
puimw | Output | real | - | |
puomw | Output | real | - | |
rho_star | Output | real | - | |
nu_star | Output | real | - | |
beta_mcdonald | Output | real | - | |
itart_r | Output | real | - | |
first_call | Input | integer | 1 |
pfcoil_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
ngrpmx | Parameter | integer | 10 | maximum number of groups of PF coils |
nclsmx | Parameter | integer | 2 | maximum number of PF coils in a given group |
nptsmx | Parameter | integer | 32 | maximum number of points across the midplane of the plasma at which the field from the PF coils is fixed |
nfixmx | Parameter | integer | 64 | maximum number of fixed current PF coils |
ngc | Parameter | integer | ngrpmx*nclsmx | maximum total number of coils across all groups |
ngc2 | Parameter | integer | ngc+2 | new variable to include 2 additional circuits: plasma and central solenoid |
alfapf | Input | real | 5e-10 | smoothing parameter used in PF coil current calculation at the beginning of pulse (BoP) |
alstroh | Input | real | 400000000.0 | allowable hoop stress in Central Solenoid structural material (Pa) |
i_cs_stress | Output | integer | - | Switch for CS stress calculation:
|
areaoh | Output | real | - | Central solenoid vertical cross-sectional area (m2) |
a_oh_turn | Output | real | - | Central solenoid (OH) trun cross-sectional area (m2) |
awpoh | Output | real | - | central solenoid conductor+void area with area of steel subtracted (m2) |
bmaxoh | Output | real | - | maximum field in central solenoid at end of flat-top (EoF) (T) |
bmaxoh0 | Output | real | - | maximum field in central solenoid at beginning of pulse (T) |
bpf | Output | real | - | peak field at coil i (T) |
ccl0_ma | Output | real | - | PF group current array, flux-swing cancellation current (MA) Input if i_pf_current=0, computed otherwise |
ccls_ma | Output | real | - | PF group current array, equilibrium current (MA) Input if i_pf_current=0, computed otherwise |
cohbop | Output | real | - | Central solenoid overall current density at beginning of pulse (A/m2) |
coheof | Input | real | 18500000.0 | Central solenoid overall current density at end of flat-top (A/m2) ( |
cpt | Output | real | - | current per turn in coil i at time j (A) |
cptdin | Input | real | [40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000. 40000.] | peak current per turn input for PF coil i (A) |
curpfb | Output | real | - | PF coil current array, at beginning of pulse (MA) Indexed by coil number, not group number |
curpff | Output | real | - | PF coil current array, at flat top (MA) Indexed by coil number, not group number |
curpfs | Output | real | - | PF coil current array, at end of pulse (MA) Indexed by coil number, not group number |
etapsu | Input | real | 0.9 | Efficiency of transfer of PF stored energy into or out of storage. |
fcohbof | Output | real | - | ratio of central solenoid overall current density at beginning of flat-top / end of flat-top |
fcohbop | Input | real | 0.9 | ratio of central solenoid overall current density at beginning of pulse / end of flat-top
( |
fcuohsu | Input | real | 0.7 | copper fraction of strand in central solenoid |
fcupfsu | Input | real | 0.69 | copper fraction of cable conductor (PF coils) |
fvssu | Input | real | 1.0 | F-value for |
ipfloc | Input | integer | [2 2 3 0 0 0 0 0 0 0] | Switch for location of PF coil group i:
|
ipfres | Output | integer | - | switch for PF & CS coil conductor type:
|
itr_sum | Output | real | - | total sum of I x turns x radius for all PF coils and CS (Am) |
isumatoh | Input | integer | 1 | switch for superconductor material in central solenoid:
|
isumatpf | Input | integer | 1 | switch for superconductor material in PF coils:
|
j_crit_str_cs | Output | real | - | superconductor strand critical current density under operating conditions in central solenoid (A/m2). Necessary for the cost calculation in $/kA m |
j_crit_str_pf | Output | real | - | superconductor strand critical current density under operating conditions in PF coils (A/m2). Necessary for the cost calculation in $/kA m |
i_pf_current | Input | integer | 1 | Switch for controlling the current of the PF coils:
|
i_sup_pf_shape | Output | integer | - | Switch for the placement of Location 3 (outboard) PF coils when the TF coils are superconducting (i_tf_sup = 1)
|
jscoh_bop | Output | real | - | central solenoid superconductor critical current density (A/m2) at beginning-of-pulse |
jscoh_eof | Output | real | - | central solenoid superconductor critical current density (A/m2) at end-of-flattop |
jcableoh_bop | Output | real | - | central solenoid cable critical current density (A/m2) at beginning-of-pulse |
jcableoh_eof | Output | real | - | central solenoid cable critical current density (A/m2) at end-of-flattop |
ncirt | Output | integer | - | number of PF circuits (including central solenoid and plasma) |
ncls | Input | integer | [1 1 2 0 0 0 0 0 0 0 0 0] | number of PF coils in group j |
nfxfh | Input | integer | 7 | number of filaments the top and bottom of the central solenoid should be broken into during scaling (5 - 10 is good) |
ngrp | Input | integer | 3 | number of groups of PF coils. Symmetric coil pairs should all be in the same group |
nohc | Output | integer | - | number of PF coils (excluding the central solenoid) + 1 |
ohhghf | Input | real | 0.71 | Central solenoid height / TF coil internal height |
oh_steel_frac | Input | real | 0.5 | central solenoid steel fraction ( |
pf_current_safety_factor | Input | real | 1.0 | Ratio of permissible PF coil conductor current density to critical conductor current density based on short-sample DC measurements |
pfcaseth | Output | real | - | steel case thickness for PF coil i (m) |
pfclres | Input | real | 2.5e-08 | PF coil resistivity (if ipfres=1) (Ohm-m) |
rhopfbus | Input | real | 3.93e-08 | Resistivity of CS and PF coil bus bars (irrespective of whether the coils themselves are superconducting or resistive) (Ohm-m) |
pfmmax | Output | real | - | mass of heaviest PF coil (tonnes) |
pfrmax | Output | real | - | radius of largest PF coil (m) |
pfwpmw | Output | real | - | Total mean wall plug power dissipated in PFC and CS power supplies (MW) (issue #713) |
powohres | Output | real | - | central solenoid resistive power during flattop (W) |
powpfres | Output | real | - | total PF coil resistive losses during flattop (W) |
ra | Output | real | - | inner radius of coil i (m) |
rb | Output | real | - | outer radius of coil i (m) |
ric | Output | real | - | peak current in coil i (MA-turns) |
rjconpf | Input | real | [30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000. 30000000.] | average winding pack current density of PF coil i (A/m2) at time of peak
current in that coil (calculated for |
rjohc | Output | real | - | allowable central solenoid current density at end of flat-top (A/m2) |
rjohc0 | Output | real | - | allowable central solenoid current density at beginning of pulse (A/m2) |
rjpfalw | Output | real | - | allowable winding pack current density of PF coil i (A/m2) |
rohc | Output | real | - | radius to the centre of the central solenoid (m) |
routr | Input | real | 1.5 | radial distance (m) from outboard TF coil leg to centre of |
rpf | Output | real | - | radius of PF coil i (m) |
rpf1 | Output | real | - | offset (m) of radial position of |
rpf2 | Input | real | -1.63 | offset (m) of radial position of |
rref | Input | real | [7. 7. 7. 7. 7. 7. 7. 7. 7. 7.] | PF coil radial positioning adjuster:
|
s_tresca_oh | Output | real | - | Maximum shear stress (Tresca criterion) coils/central solenoid [MPa] |
sigpfcalw | Input | real | 500.0 | maximum permissible tensile stress (MPa) in steel coil cases for superconducting
PF coils ( |
sigpfcf | Input | real | 1.0 | fraction of JxB hoop force supported by steel case for superconducting PF coils ( |
sxlg | Output | real | - | mutual inductance matrix (H) |
tmargoh | Output | real | - | Central solenoid temperature margin (K) |
turns | Output | real | - | number of turns in PF coil i |
vf | Input | real | [0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3 0.3] | winding pack void fraction of PF coil i for coolant |
vfohc | Input | real | 0.3 | void fraction of central solenoid conductor for coolant |
vsbn | Output | real | - | total flux swing available for burn (Wb) |
vsefbn | Output | real | - | flux swing from PF coils for burn (Wb) |
vsefsu | Output | real | - | flux swing from PF coils for startup (Wb) |
vseft | Output | real | - | total flux swing from PF coils (Wb) |
vsoh | Output | real | - | total flux swing from the central solenoid (Wb) |
vsohbn | Output | real | - | central solenoid flux swing for burn (Wb) |
vsohsu | Output | real | - | central solenoid flux swing for startup (Wb) |
vssu | Output | real | - | total flux swing for startup ( |
vstot | Output | real | - | total flux swing for pulse (Wb) |
waves | Output | real | - | used in current waveform of PF coils/central solenoid |
whtpf | Output | real | - | total mass of the PF coil conductor (kg) |
whtpfs | Output | real | - | total mass of the PF coil structure (kg) |
wtc | Output | real | - | conductor mass for PF coil i (kg) |
wts | Output | real | - | structure mass for PF coil i (kg) |
zh | Output | real | - | upper point of PF coil i (m) |
zl | Output | real | - | lower point of PF coil i (m) |
zpf | Output | real | - | z (height) location of PF coil i (m) |
zref | Input | real | [3.6 1.2 2.5 1. 1. 1. 1. 1. 1. 1. ] | PF coil vertical positioning adjuster:
|
bmaxcs_lim | Input | real | 13.0 | Central solenoid max field limit [T] |
fbmaxcs | Input | real | 1.0 | F-value for CS mmax field ( |
ld_ratio_cst | Input | real | 3.0 | Ratio of CS coil turn conduit length to depth |
l_cond_cst | Output | real | - | Length of CS of CS coil turn conduit |
d_cond_cst | Output | real | - | Depth/width of CS of CS coil turn conduit |
r_out_cst | Input | real | 0.003 | Length of CS of CS coil turn conduit length |
r_in_cst | Output | real | - | Length of CS of CS coil turn conduit length |
global_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
icase | Input | character | b'Steady-state tokamak model ' | power plant type |
runtitle | Input | character | b"Run Title (change this line using input variable 'runtitle') " | short descriptive title for the run |
verbose | Output | integer | - | switch for turning on/off diagnostic messages
|
run_tests | Output | integer | - | turns on built-in tests if set to 1 |
maxcal | Input | integer | 200 | maximum number of VMCON iterations |
fileprefix | Input | character | b' ' | input file prefix |
output_prefix | Input | character | b' ' | output file prefix |
xlabel | Input | character | b' ' | scan parameter description label |
vlabel | Input | character | b' ' | scan value name label |
xlabel_2 | Input | character | b' ' | scan parameter description label (2nd dimension) |
vlabel_2 | Input | character | b' ' | scan value name label (2nd dimension) |
iscan_global | Output | integer | - | Makes iscan available globally. |
convergence_parameter | Output | real | - | VMCON convergence parameter "sum" |
build_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
aplasmin | Input | real | 0.25 | minimum minor radius (m) |
available_radial_space | Output | real | - | Minimal radial space between plasma and coils (m) |
blarea | Output | real | - | blanket total surface area (m2) |
blareaib | Output | real | - | inboard blanket surface area (m2) |
blareaob | Output | real | - | outboard blanket surface area (m2) |
blbmith | Input | real | 0.17 | inboard blanket box manifold thickness (m) ( |
blbmoth | Input | real | 0.27 | outboard blanket box manifold thickness (m) ( |
blbpith | Input | real | 0.3 | inboard blanket base plate thickness (m) ( |
blbpoth | Input | real | 0.35 | outboard blanket base plate thickness (m) ( |
blbuith | Input | real | 0.365 | inboard blanket breeding zone thickness (m) ( |
blbuoth | Input | real | 0.465 | outboard blanket breeding zone thickness (m) ( |
blnkith | Input | real | 0.115 | inboard blanket thickness (m); (calculated if |
blnkoth | Input | real | 0.235 | outboard blanket thickness (m); calculated if |
blnktth | Output | real | - | top blanket thickness (m), = mean of inboard and outboard blanket thicknesses |
bore | Input | real | 1.42 | central solenoid inboard radius (m) ( |
clhsf | Input | real | 4.268 | cryostat lid height scaling factor (tokamaks) |
ddwex | Input | real | 0.07 | cryostat thickness (m) |
d_vv_in | Input | real | 0.07 | vacuum vessel inboard thickness (TF coil / shield) (m) |
d_vv_out | Input | real | 0.07 | vacuum vessel outboard thickness (TF coil / shield) (m) |
d_vv_top | Input | real | 0.07 | vacuum vessel topside thickness (TF coil / shield) (m) (= d_vv_bot if double-null) |
d_vv_bot | Input | real | 0.07 | vacuum vessel underside thickness (TF coil / shield) (m) |
f_avspace | Input | real | 1.0 | F-value for stellarator radial space check ( |
fcspc | Input | real | 0.6 | Fraction of space occupied by CS pre-compression structure |
fseppc | Input | real | 350000000.0 | Separation force in CS coil pre-compression structure |
fwarea | Output | real | - | first wall total surface area (m2) |
fwareaib | Output | real | - | inboard first wall surface area (m2) |
fwareaob | Output | real | - | outboard first wall surface area (m2) |
fwith | Output | real | - | inboard first wall thickness, initial estimate as calculated (m) |
fwoth | Output | real | - | outboard first wall thickness, initial estimate as calculated (m) |
gapds | Input | real | 0.155 | gap between inboard vacuum vessel and thermal shield (m) ( |
gapoh | Input | real | 0.08 | gap between central solenoid and TF coil (m) ( |
gapomin | Input | real | 0.234 | minimum gap between outboard vacuum vessel and TF coil (m) ( |
gapsto | Output | real | - | gap between outboard vacuum vessel and TF coil (m) |
hmax | Output | real | - | maximum (half-)height of TF coil (inside edge) (m) |
hpfdif | Output | real | - | difference in distance from midplane of upper and lower portions of TF legs (non-zero for single-null devices) (m) |
hpfu | Output | real | - | height to top of (upper) TF coil leg (m) |
hr1 | Output | real | - | half-height of TF coil inboard leg straight section (m) |
iohcl | Input | integer | 1 | Switch for existence of central solenoid:
|
iprecomp | Input | integer | 1 | Switch for existence of central solenoid pre-compression structure:
|
tf_in_cs | Output | integer | - | Switch for placing the TF coil inside the CS
|
ohcth | Input | real | 0.811 | Central solenoid thickness (m) ( |
precomp | Output | real | - | CS coil precompression structure thickness (m) |
rbld | Output | real | - | sum of thicknesses to the major radius (m) |
required_radial_space | Output | real | - | Required space between coil and plasma for blanket shield wall etc (m) |
rinboard | Input | real | 0.651 | plasma inboard radius (m) ( |
rsldi | Output | real | - | radius to inboard shield (inside point) (m) |
rsldo | Output | real | - | radius to outboard shield (outside point) (m) |
r_vv_inboard_out | Output | real | - | Radial plasma facing side position of inboard vacuum vessel [m] |
r_sh_inboard_in | Output | real | - | Radial inner side position of inboard neutronic shield [m] |
r_sh_inboard_out | Output | real | - | Radial plasma facing side position of inboard neutronic shield [m] |
r_tf_inboard_in | Output | real | - | Mid-plane inboard TF coil leg radius at the centre-machine side [m] |
r_tf_inboard_mid | Output | real | - | Mid-plane inboard TF coil leg radius at middle of the coil [m] |
r_tf_inboard_out | Output | real | - | Mid-plane inboard TF coil leg radius at the plasma side [m] |
r_tf_outboard_mid | Output | real | - | Mid-plane outboard TF coil leg radius at the middle of the coil [m] |
i_r_cp_top | Output | integer | - | Switch selecting the he parametrization of the outer radius of the top of the CP part of the TF coil
0 : |
r_cp_top | Output | real | - | Top outer radius of the centropost (ST only) (m) |
f_r_cp | Input | real | 1.4 | Ratio between the top and the midplane TF CP outer radius [-] Not used by default (-1) must be larger than 1 otherwise |
dr_tf_inner_bore | Output | real | - | TF coil horizontal inner bore (m) |
dh_tf_inner_bore | Output | real | - | TF coil vertical inner bore (m) |
scrapli | Input | real | 0.14 | Gap between plasma and first wall, inboard side (m) (if |
scraplo | Input | real | 0.15 | Gap between plasma and first wall, outboard side (m) (if |
sharea | Output | real | - | shield total surface area (m2) |
shareaib | Output | real | - | inboard shield surface area (m2) |
shareaob | Output | real | - | outboard shield surface area (m2) |
shldith | Input | real | 0.69 | inboard shield thickness (m) ( |
shldlth | Input | real | 0.7 | lower (under divertor) shield thickness (m) |
shldoth | Input | real | 1.05 | outboard shield thickness (m) ( |
shldtth | Input | real | 0.6 | upper/lower shield thickness (m); calculated if |
sigallpc | Input | real | 300000000.0 | allowable stress in CSpre-compression structure (Pa) |
tfcth | Output | real | - | inboard TF coil thickness, (centrepost for ST) (m)
(input, calculated or |
tfoffset | Output | real | - | vertical distance between centre of TF coils and centre of plasma (m) |
tfootfi | Input | real | 1.19 | TF coil outboard leg / inboard leg radial thickness
ratio ( |
tfthko | Output | real | - | Outboard TF coil thickness (m) |
tftsgap | Input | real | 0.05 | Minimum metal-to-metal gap between TF coil and thermal shield (m) |
thshield_ib | Input | real | 0.05 | TF-VV thermal shield thickness, inboard (m) |
thshield_ob | Input | real | 0.05 | TF-VV thermal shield thickness, outboard (m) |
thshield_vb | Input | real | 0.05 | TF-VV thermal shield thickness, vertical build (m) |
vgap_vv_thermalshield | Input | real | 0.163 | vertical gap between vacuum vessel and thermal shields (m) |
vgap_xpoint_divertor | Output | real | - | vertical gap between x-point and divertor (m) (if = 0, it is calculated) |
vgaptop | Input | real | 0.6 | vertical gap between top of plasma and first wall (m) (= vgap_xpoint_divertor if double-null) |
vvblgap | Input | real | 0.05 | gap between vacuum vessel and blanket (m) |
plleni | Input | real | 1.0 | length of inboard divertor plate (m) |
plleno | Input | real | 1.0 | length of outboard divertor plate (m) |
plsepi | Input | real | 1.0 | poloidal length, x-point to inboard strike point (m) |
plsepo | Input | real | 1.5 | poloidal length, x-point to outboard strike point (m) |
rspo | Output | real | - | outboard strike point radius (m) |
heat_transport_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
baseel | Input | real | 5000000.0 | base plant electric load (W) |
crypmw | Output | real | - | cryogenic plant power (MW) |
crypmw_max | Input | real | 50.0 | Maximum cryogenic plant power (MW) Constraint equation icc = 87 Scan variable nwseep = 56 |
f_crypmw | Input | real | 1.0 | f-value for maximum cryogenic plant power Iteration variable ixc = 164 Constraint equation icc = 87 |
etatf | Input | real | 0.9 | AC to resistive power conversion for TF coils |
etath | Input | real | 0.35 | thermal to electric conversion efficiency if |
etath_liq | Input | real | 0.35 | |
fachtmw | Output | real | - | facility heat removal (MW) |
fcsht | Output | real | - | total baseline power required at all times (MW) |
fgrosbop | Output | real | - | scaled fraction of gross power to balance-of-plant |
fmgdmw | Output | real | - | power to mgf (motor-generator flywheel) units (MW) (ignored if |
fpumpblkt | Input | real | 0.005 | fraction of total blanket thermal power required to drive the blanket
coolant pumps (default assumes water coolant) ( |
fpumpdiv | Input | real | 0.005 | fraction of total divertor thermal power required to drive the divertor coolant pumps (default assumes water coolant) |
fpumpfw | Input | real | 0.005 | fraction of total first wall thermal power required to drive the FW coolant
pumps (default assumes water coolant) ( |
fpumpshld | Input | real | 0.005 | fraction of total shield thermal power required to drive the shield coolant pumps (default assumes water coolant) |
htpmw_min | Output | real | - | Minimum total electrical power for primary coolant pumps (MW) (NOT RECOMMENDED) |
helpow | Output | real | - | Heat removal at cryogenic temperature tmpcry (W) |
helpow_cryal | Output | real | - | Heat removal at cryogenic temperature tcoolin (W) |
htpmw | Output | real | - | heat transport system electrical pump power (MW) |
htpmw_blkt | Output | real | - | blanket primary coolant mechanical pumping power (MW) |
htpmw_blkt_liq | Output | real | - | blanket secondary coolant mechanical pumping power (MW) |
htpmw_blkt_tot | Output | real | - | blanket primary + secondary coolant mechanical pumping power (MW) |
htpmw_div | Output | real | - | divertor coolant mechanical pumping power (MW) |
htpmw_fw | Output | real | - | first wall coolant mechanical pumping power (MW) |
htpmw_shld | Output | real | - | shield and vacuum vessel coolant mechanical pumping power (MW) |
htpsecmw | Output | real | - | Waste power lost from primary coolant pumps (MW) |
ipowerflow | Input | integer | 1 | switch for power flow model:
|
iprimshld | Input | integer | 1 | Switch for shield thermal power destiny:
|
nphx | Output | integer | - | number of primary heat exchangers |
pacpmw | Output | real | - | total pulsed power system load (MW) |
peakmva | Output | real | - | peak MVA requirement |
pfwdiv | Output | real | - | heat removal from first wall/divertor (MW) |
pgrossmw | Output | real | - | gross electric power (MW) |
pinjht | Output | real | - | power dissipated in heating and current drive system (MW) |
pinjmax | Input | real | 120.0 | maximum injector power during pulse (heating and ramp-up/down phase) (MW) |
pinjwp | Output | real | - | injector wall plug power (MW) |
pinjwpfix | Output | real | - | secondary injector wall plug power (MW) |
pnetelmw | Output | real | - | net electric power (MW) |
precircmw | Output | real | - | recirculating electric power (MW) |
priheat | Output | real | - | total thermal power removed from fusion core (MW) |
psecdiv | Output | real | - | Low-grade heat lost in divertor (MW) |
psechcd | Output | real | - | Low-grade heat lost into HCD apparatus (MW) |
psechtmw | Output | real | - | Low-grade heat (MW) |
pseclossmw | Output | real | - | Low-grade heat (VV + lost)(MW) |
psecshld | Output | real | - | Low-grade heat deposited in shield (MW) |
pthermmw | Output | real | - | High-grade heat useful for electric production (MW) |
pwpm2 | Input | real | 150.0 | base AC power requirement per unit floor area (W/m2) |
tfacpd | Output | real | - | total steady state TF coil AC power demand (MW) |
tlvpmw | Output | real | - | estimate of total low voltage power (MW) |
trithtmw | Input | real | 15.0 | power required for tritium processing (MW) |
tturb | Output | real | - | coolant temperature at turbine inlet (K) ( |
vachtmw | Input | real | 0.5 | vacuum pump power (MW) |
vacuum_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
vacuum_model | Input | character | b'old ' | switch for vacuum pumping model:
|
niterpump | Output | real | - | number of high vacuum pumps (real number), each with the throughput of one
ITER cryopump (50 Pa m3 s-1), all operating at the same time ( |
ntype | Input | integer | 1 | switch for vacuum pump type:
|
nvduct | Output | integer | - | number of ducts (torus to pumps) |
dlscal | Output | real | - | vacuum system duct length scaling |
pbase | Input | real | 0.0005 | base pressure during dwell before gas pre-fill(Pa) |
prdiv | Input | real | 0.36 | divertor chamber pressure during burn (Pa) |
pumptp | Input | real | 1.2155e+22 | Pump throughput (molecules/s) (default is ITER value) |
rat | Input | real | 1.3e-08 | plasma chamber wall outgassing rate (Pa-m/s) |
tn | Input | real | 300.0 | neutral gas temperature in chamber (K) |
vacdshm | Output | real | - | mass of vacuum duct shield (kg) |
vcdimax | Output | real | - | diameter of duct passage (m) |
vpumpn | Output | integer | - | number of high vacuum pumps |
dwell_pump | Output | integer | - | switch for dwell pumping options:
|
pumpareafraction | Input | real | 0.0203 | area of one pumping port as a fraction of plasma surface area |
pumpspeedmax | Input | real | 27.3 | maximum pumping speed per unit area for deuterium & tritium, molecular flow |
pumpspeedfactor | Input | real | 0.167 | effective pumping speed reduction factor due to duct impedance |
initialpressure | Input | real | 1.0 | initial neutral pressure at the beginning of the dwell phase (Pa) |
outgasindex | Input | real | 1.0 | outgassing decay index |
outgasfactor | Input | real | 0.0235 | outgassing prefactor kw: outgassing rate at 1 s per unit area (Pa m s-1) |
dcll_module
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
r_fci | Output | real | - | Radial BZ thickness [m] |
r_backwall | Output | real | - | Radial BZ thickness [m] |
bz_r_ib | Output | real | - | Structure/coolant compositional fractions |
bz_r_ob | Output | real | - | Structure/coolant compositional fractions |
f_vol_stff_plates | Output | real | - | MF/BSS compositional fractions |
f_vol_stl_bz_struct | Output | real | - | MF/BSS compositional fractions |
f_vol_stl_back_wall | Output | real | - | MF/BSS compositional fractions |
f_vol_stl_fw | Output | real | - | MF/BSS compositional fractions |
f_vol_mfbss_stl | Output | real | - | Volume of FCIs, other BZ structure, liquid channels, backwall and MF/BSS [m^3] |
f_vol_mfbss_he | Output | real | - | Volume of FCIs, other BZ structure, liquid channels, backwall and MF/BSS [m^3] |
f_vol_mfbss_pbli | Output | real | - | Volume of FCIs, other BZ structure, liquid channels, backwall and MF/BSS [m^3] |
vol_fci | Output | real | - | BZ masses by composition [kg] |
vol_bz_struct | Output | real | - | BZ masses by composition [kg] |
vol_bz_liq | Output | real | - | BZ masses by composition [kg] |
vol_bz_liq_ib | Output | real | - | BZ masses by composition [kg] |
vol_bz_liq_ob | Output | real | - | BZ masses by composition [kg] |
vol_bw | Output | real | - | BZ masses by composition [kg] |
vol_bss | Output | real | - | BZ masses by composition [kg] |
wht_cer | Output | real | - | Backwall masses by composition [kg] |
wht_stl_struct | Output | real | - | Backwall masses by composition [kg] |
wht_cool_struct | Output | real | - | Backwall masses by composition [kg] |
wht_bw_stl | Output | real | - | MF/BSS masses by composition [kg] |
wht_bw_cool | Output | real | - | MF/BSS masses by composition [kg] |
wht_mfbss_stl | Output | real | - | FW masses by composition [kg] |
wht_mfbss_cool | Output | real | - | FW masses by composition [kg] |
wht_mfbss_pbli | Output | real | - | FW masses by composition [kg] |
fwmass_stl | Output | real | - | Total masses of material in blanket [kg] |
fwmass_cool | Output | real | - | Total masses of material in blanket [kg] |
mass_cool_blanket | Output | real | - | Total mass for an inboard/outboard reactor segment [kg] |
mass_liq_blanket | Output | real | - | Total mass for an inboard/outboard reactor segment [kg] |
mass_stl_blanket | Output | real | - | Total mass for an inboard/outboard reactor segment [kg] |
mass_segm_ib | Output | real | - | |
mass_segm_ob | Output | real | - |
sctfcoil_module
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
tf_fit_t | Output | real | - | Dimensionless winding pack width |
tf_fit_z | Output | real | - | Dimensionless winding pack radial thickness |
tf_fit_y | Output | real | - | Ratio of peak field with ripple to nominal axisymmetric peak field |
tfc_current | Output | real | - | Current in each TF coil |
awpc | Output | real | - | Total cross-sectional area of winding pack including GW insulation and insertion gap [m2] |
awptf | Output | real | - | Total cross-sectional area of winding pack without ground insulation and insertion gap [m2] |
a_tf_steel | Output | real | - | Inboard coil steel coil cross-sectional area [m2] |
a_tf_ins | Output | real | - | Inboard coil insulation cross-section per coil [m2] |
f_tf_steel | Output | real | - | Inboard coil steel fraction [-] |
f_tf_ins | Output | real | - | Inboard coil insulation fraction [-] |
h_cp_top | Output | real | - | Vertical distance from the midplane to the top of the tapered section [m] |
r_tf_outboard_in | Output | real | - | Radial position of plasma-facing edge of TF coil outboard leg [m] |
r_tf_outboard_out | Output | real | - | Radial position of outer edge of TF coil inboard leg [m] |
r_wp_inner | Output | real | - | Radial position of inner edge and centre of winding pack [m] |
r_wp_outer | Output | real | - | Radial position of outer edge and centre of winding pack [m] |
r_wp_centre | Output | real | - | Radial position of centre and centre of winding pack [m] |
dr_tf_wp_top | Output | real | - | Conductor layer radial thickness at centercollumn top [m] Ground insulation layer included, only defined for itart = 1 |
vol_ins_cp | Output | real | - | CP turn insulation volume [m3] |
vol_gr_ins_cp | Output | real | - | CP ground insulation volume [m3] |
vol_case_cp | Output | real | - | Volume of the CP outer casing cylinder |
t_wp_toroidal | Output | real | - | Minimal toroidal thickness of of winding pack [m] |
t_wp_toroidal_av | Output | real | - | Averaged toroidal thickness of of winding pack [m] |
t_lat_case_av | Output | real | - | Average lateral casing thickness [m] |
a_case_front | Output | real | - | Front casing area [m2] |
a_case_nose | Output | real | - | Nose casing area [m2] |
a_ground_ins | Output | real | - | Inboard mid-plane cross-section area of the WP ground insulation [m2] |
a_leg_ins | Output | real | - | TF ouboard leg turn insulation area per coil [m2] |
a_leg_gr_ins | Output | real | - | TF outboard leg ground insulation area per coil [m2] |
a_leg_cond | Output | real | - | Exact TF ouboard leg conductor area [m2] |
theta_coil | Output | real | - | Half toroidal angular extent of a single TF coil inboard leg |
tan_theta_coil | Output | real | - | Tan half toroidal angular extent of a single TF coil inboard leg |
t_conductor_radial | Output | real | - | Conductor area radial and toroidal dimension (integer turn only) [m] |
t_conductor_toroidal | Output | real | - | Conductor area radial and toroidal dimension (integer turn only) [m] |
t_cable_radial | Output | real | - | Cable area radial and toroidal dimension (integer turn only) [m] |
t_cable_toroidal | Output | real | - | Cable area radial and toroidal dimension (integer turn only) [m] |
t_turn_radial | Output | real | - | Turn radial and toroidal dimension (integer turn only) [m] |
t_turn_toroidal | Output | real | - | Turn radial and toroidal dimension (integer turn only) [m] |
t_cable | Output | real | - | Cable area averaged dimension (square shape) [m] |
vforce_inboard_tot | Output | real | - | Total inboard vertical tension (all coils) [N] |
vv_stress_quench | Output | real | - | The Tresca stress experienced by the Vacuum Vessel when the SCTF coil quenches [Pa] |
copper* | Variable | type | - | |
hastelloy* | Variable | type | - | |
solder* | Variable | type | - | |
jacket* | Variable | type | - | |
helium* | Variable | type | - | |
croco_strand_area | Output | real | - | |
croco_strand_critical_current | Output | real | - | |
conductor_copper_area | Output | real | - | |
conductor_copper_fraction | Output | real | - | |
conductor_copper_bar_area | Output | real | - | |
conductor_hastelloy_area | Output | real | - | |
conductor_hastelloy_fraction | Output | real | - | |
conductor_helium_area | Output | real | - | |
conductor_helium_fraction | Output | real | - | |
conductor_solder_area | Output | real | - | |
conductor_solder_fraction | Output | real | - | |
conductor_jacket_area | Output | real | - | |
conductor_jacket_fraction | Output | real | - | |
conductor_rebco_area | Output | real | - | |
conductor_rebco_fraction | Output | real | - | |
conductor_critical_current | Output | real | - | |
conductor_acs | Output | real | - | |
conductor_area | Output | real | - | Area of cable space inside jacket |
t1 | Output | real | - | |
time2 | Output | real | - | |
tau2 | Output | real | - | |
estotft | Output | real | - | |
is_leg_cp_temp_same | Output | integer | - |
stellarator_module
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
f_n | Output | real | - | |
f_r | Output | real | - | |
f_aspect | Output | real | - | |
f_b | Output | real | - | |
f_i | Output | real | - | |
f_a | Output | real | - | |
first_call | Variable | logical | .true. | |
first_call_stfwbs | Variable | logical | .true. |
buildings_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
admv | Input | real | 100000.0 | administration building volume (m3) |
admvol | Output | real | - | volume of administration buildings (m3) |
aux_build_l | Input | real | 60.0 | aux building supporting tokamak processes length, width, height (m) |
aux_build_w | Input | real | 30.0 | aux building supporting tokamak processes length, width, height (m) |
aux_build_h | Input | real | 5.0 | aux building supporting tokamak processes length, width, height (m) |
auxcool_l | Input | real | 20.0 | Site-Wide Auxiliary Cooling Water facility length, width, height (m) |
auxcool_w | Input | real | 20.0 | Site-Wide Auxiliary Cooling Water facility length, width, height (m) |
auxcool_h | Input | real | 5.0 | Site-Wide Auxiliary Cooling Water facility length, width, height (m) |
bioshld_thk | Input | real | 2.5 | Radial thickness of bio-shield around reactor (m) |
chemlab_l | Input | real | 50.0 | Chemistry labs and treatment buldings length, width, height (m) |
chemlab_w | Input | real | 30.0 | Chemistry labs and treatment buldings length, width, height (m) |
chemlab_h | Input | real | 6.0 | Chemistry labs and treatment buldings length, width, height (m) |
clh1 | Input | real | 2.5 | vertical clearance from TF coil to cryostat (m) (calculated for tokamaks) |
clh2 | Input | real | 15.0 | clearance beneath TF coil to foundation (including basement) (m) |
control_buildings_l | Input | real | 80.0 | control building length, width, height (m) |
control_buildings_w | Input | real | 60.0 | control building length, width, height (m) |
control_buildings_h | Input | real | 6.0 | control building length, width, height (m) |
conv | Input | real | 60000.0 | control building volume (m3) |
convol | Output | real | - | volume of control, protection and i&c building (m3) |
crane_arm_h | Input | real | 10.0 | vertical dimension of crane arm, operating over reactor (m) |
crane_clrnc_h | Input | real | 4.0 | horizontal clearance to building wall for crane operation (m) |
crane_clrnc_v | Input | real | 3.0 | vertical clearance for crane operation (m) |
cryomag_l | Input | real | 120.0 | Cryogenic Buildings for Magnet and Fuel Cycle length, width, height (m) |
cryomag_w | Input | real | 90.0 | Cryogenic Buildings for Magnet and Fuel Cycle length, width, height (m) |
cryomag_h | Input | real | 5.0 | Cryogenic Buildings for Magnet and Fuel Cycle length, width, height (m) |
cryostore_l | Input | real | 160.0 | Magnet Cryo Storage Tanks length, width, height (m) |
cryostore_w | Input | real | 30.0 | Magnet Cryo Storage Tanks length, width, height (m) |
cryostore_h | Input | real | 20.0 | Magnet Cryo Storage Tanks length, width, height (m) |
cryostat_clrnc | Input | real | 2.5 | vertical clearance from TF coil to cryostat (m) |
cryvol | Output | real | - | volume of cryoplant building (m3) |
efloor | Output | real | - | effective total floor space (m2) |
elecdist_l | Input | real | 380.0 | Transformers and electrical distribution facilities length, width, height (m) |
elecdist_w | Input | real | 350.0 | Transformers and electrical distribution facilities length, width, height (m) |
elecdist_h | Input | real | 5.0 | Transformers and electrical distribution facilities length, width, height (m) |
elecload_l | Input | real | 100.0 | Electric (eesential and non-essential) load centres length, width, height (m) |
elecload_w | Input | real | 90.0 | Electric (eesential and non-essential) load centres length, width, height (m) |
elecload_h | Input | real | 3.0 | Electric (eesential and non-essential) load centres length, width, height (m) |
elecstore_l | Input | real | 100.0 | Energy Storage facilities length, width, height (m) |
elecstore_w | Input | real | 60.0 | Energy Storage facilities length, width, height (m) |
elecstore_h | Input | real | 12.0 | Energy Storage facilities length, width, height (m) |
elevol | Output | real | - | volume of electrical equipment building (m3) |
esbldgm3 | Input | real | 1000.0 | volume of energy storage equipment building (m3) (not used if |
fc_building_l | Input | real | 60.0 | Fuel Cycle facilities length, width (m) |
fc_building_w | Input | real | 60.0 | Fuel Cycle facilities length, width (m) |
fndt | Input | real | 2.0 | foundation thickness (m) |
gas_buildings_l | Input | real | 25.0 | air & gas supply (amalgamated) buildings length, width, height (m) |
gas_buildings_w | Input | real | 15.0 | air & gas supply (amalgamated) buildings length, width, height (m) |
gas_buildings_h | Input | real | 5.0 | air & gas supply (amalgamated) buildings length, width, height (m) |
ground_clrnc | Input | real | 5.0 | clearance beneath TF coil (m) |
hcd_building_l | Input | real | 70.0 | HCD building length, width, height (m) |
hcd_building_w | Input | real | 40.0 | HCD building length, width, height (m) |
hcd_building_h | Input | real | 25.0 | HCD building length, width, height (m) |
hccl | Input | real | 5.0 | clearance around components in hot cell (m) |
hcwt | Input | real | 1.5 | hot cell wall thickness (m) |
heat_sink_l | Input | real | 160.0 | heat sinks length, width, height (m) |
heat_sink_w | Input | real | 80.0 | heat sinks length, width, height (m) |
heat_sink_h | Input | real | 12.0 | heat sinks length, width, height (m) |
hot_sepdist | Input | real | 2.0 | hot cell storage component separation distance (m) |
hotcell_h | Input | real | 12.0 | hot cell storage and maintenance facility height (m) |
hw_storage_l | Input | real | 20.0 | hazardous waste storage building length, width, height (m) |
hw_storage_w | Input | real | 10.0 | hazardous waste storage building length, width, height (m) |
hw_storage_h | Input | real | 5.0 | hazardous waste storage building length, width, height (m) |
i_bldgs_size | Output | integer | - | switch between routines estimating building sizes (0 = default; 1 = updated) |
i_bldgs_v | Output | integer | - | switch to select verbose output for buildings (1 = verbose) |
ilw_smelter_l | Input | real | 50.0 | radioactive waste smelting facility length, width, height (m) |
ilw_smelter_w | Input | real | 30.0 | radioactive waste smelting facility length, width, height (m) |
ilw_smelter_h | Input | real | 30.0 | radioactive waste smelting facility length, width, height (m) |
ilw_storage_l | Input | real | 120.0 | ILW waste storage building length, width, height (m) |
ilw_storage_w | Input | real | 100.0 | ILW waste storage building length, width, height (m) |
ilw_storage_h | Input | real | 8.0 | ILW waste storage building length, width, height (m) |
llw_storage_l | Input | real | 45.0 | LLW waste storage building length, width, height (m) |
llw_storage_w | Input | real | 20.0 | LLW waste storage building length, width, height (m) |
llw_storage_h | Input | real | 5.0 | LLW waste storage building length, width, height (m) |
magnet_pulse_l | Input | real | 105.0 | pulsed magnet power building length, width, height (m) |
magnet_pulse_w | Input | real | 40.0 | pulsed magnet power building length, width, height (m) |
magnet_pulse_h | Input | real | 5.0 | pulsed magnet power building length, width, height (m) |
magnet_trains_l | Input | real | 120.0 | steady state magnet power trains building length, width, height (m) |
magnet_trains_w | Input | real | 90.0 | steady state magnet power trains building length, width, height (m) |
magnet_trains_h | Input | real | 5.0 | steady state magnet power trains building length, width, height (m) |
maint_cont_l | Input | real | 125.0 | maintenance control building length, width, height (m) |
maint_cont_w | Input | real | 100.0 | maintenance control building length, width, height (m) |
maint_cont_h | Input | real | 6.0 | maintenance control building length, width, height (m) |
mbvfac | Input | real | 2.8 | maintenance building volume multiplication factor |
nbi_sys_l | Input | real | 225.0 | NBI system length, width (m) |
nbi_sys_w | Input | real | 185.0 | NBI system length, width (m) |
pfbldgm3 | Input | real | 20000.0 | volume of PF coil power supply building (m3) |
pibv | Input | real | 20000.0 | power injection building volume (m3) |
qnty_sfty_fac | Input | real | 2.0 | quantity safety factor for component use during plant lifetime |
rbvfac | Input | real | 1.6 | reactor building volume multiplication factor |
rbrt | Input | real | 1.0 | reactor building roof thickness (m) |
rbvol | Output | real | - | reactor building volume (m3) |
rbwt | Input | real | 2.0 | reactor building wall thickness (m) |
reactor_clrnc | Input | real | 4.0 | clearance around reactor (m) |
reactor_fndtn_thk | Input | real | 2.0 | reactor building foundation thickness (m) |
reactor_hall_l | Output | real | - | reactor building length, width, height (m) |
reactor_hall_w | Output | real | - | reactor building length, width, height (m) |
reactor_hall_h | Output | real | - | reactor building length, width, height (m) |
reactor_roof_thk | Input | real | 1.0 | reactor building roof thickness (m) |
reactor_wall_thk | Input | real | 2.0 | reactor building wall thickness (m) |
rmbvol | Output | real | - | volume of maintenance and assembly building (m3) |
robotics_l | Input | real | 50.0 | robotics buildings length, width, height (m) |
robotics_w | Input | real | 30.0 | robotics buildings length, width, height (m) |
robotics_h | Input | real | 30.0 | robotics buildings length, width, height (m) |
row | Input | real | 4.0 | clearance to building wall for crane operation (m) |
rxcl | Input | real | 4.0 | clearance around reactor (m) |
sec_buildings_l | Input | real | 30.0 | security & safety buildings length, width, height (m) |
sec_buildings_w | Input | real | 25.0 | security & safety buildings length, width, height (m) |
sec_buildings_h | Input | real | 6.0 | security & safety buildings length, width, height (m) |
shmf | Input | real | 0.5 | fraction of shield mass per TF coil to be moved in the maximum shield lift |
shov | Input | real | 100000.0 | shops and warehouse volume (m3) |
shovol | Output | real | - | volume of shops and buildings for plant auxiliaries (m3) |
staff_buildings_area | Input | real | 480000.0 | footprint of staff buildings (m2) |
staff_buildings_h | Input | real | 5.0 | staff buildings height (m) |
stcl | Input | real | 3.0 | clearance above crane to roof (m) |
tfcbv | Input | real | 20000.0 | volume of TF coil power supply building (m3) (calculated if TF coils are superconducting) |
transp_clrnc | Input | real | 1.0 | transportation clearance between components (m) |
trcl | Input | real | 1.0 | transportation clearance between components (m) |
triv | Input | real | 40000.0 | volume of tritium, fuel handling and health physics buildings (m3) |
turbine_hall_l | Input | real | 109.0 | turbine hall length, width, height (m) |
turbine_hall_w | Input | real | 62.0 | turbine hall length, width, height (m) |
turbine_hall_h | Input | real | 15.0 | turbine hall length, width, height (m) |
tw_storage_l | Input | real | 90.0 | tritiated waste storage building length, width, height (m) |
tw_storage_w | Input | real | 30.0 | tritiated waste storage building length, width, height (m) |
tw_storage_h | Input | real | 5.0 | tritiated waste storage building length, width, height (m) |
volrci | Output | real | - | internal volume of reactor building (m3) |
volnucb | Output | real | - | sum of nuclear buildings volumes (m3) |
warm_shop_l | Input | real | 100.0 | warm shop length, width, height (m) |
warm_shop_w | Input | real | 50.0 | warm shop length, width, height (m) |
warm_shop_h | Input | real | 10.0 | warm shop length, width, height (m) |
water_buildings_l | Input | real | 110.0 | water, laundry & drainage buildings length, width, height (m) |
water_buildings_w | Input | real | 10.0 | water, laundry & drainage buildings length, width, height (m) |
water_buildings_h | Input | real | 5.0 | water, laundry & drainage buildings length, width, height (m) |
wgt | Input | real | 500000.0 | reactor building crane capacity (kg) (calculated if 0 is input) |
wgt2 | Input | real | 100000.0 | hot cell crane capacity (kg) (calculated if 0 is input) |
workshop_l | Input | real | 150.0 | [cold] workshop buildings length, width, height (m) |
workshop_w | Input | real | 125.0 | [cold] workshop buildings length, width, height (m) |
workshop_h | Input | real | 10.0 | [cold] workshop buildings length, width, height (m) |
wrbi | Output | real | - | distance from centre of machine to building wall (m) |
wsvol | Output | real | - | volume of warm shop building (m3) |
wsvfac | Input | real | 1.9 | warm shop building volume multiplication factor |
a_reactor_bldg | Input | real | 8320.0 | Floor area of reactor building in m^2 |
a_ee_ps_bldg | Input | real | 21330.0 | Floor area of electrical equipment and power supply building in m^2 |
a_aux_services_bldg | Input | real | 1000.0 | Floor area of auxiliary services building in m^2 |
a_hot_cell_bldg | Input | real | 8430.0 | Floor area of hot cell building in m^2 |
a_reactor_service_bldg | Input | real | 2440.0 | Floor area of reactor service building in m^2 |
a_service_water_bldg | Input | real | 1567.0 | Floor area of service water building in m^2 |
a_fuel_handling_bldg | Input | real | 1670.0 | Floor area of fuel handling and storage building in m^2 |
a_control_room_bldg | Input | real | 2880.0 | Floor area of controlroom building in m^2 |
a_ac_ps_bldg | Input | real | 6423.0 | Floor area of AC power supply building in m^2 |
a_admin_bldg | Input | real | 25674.0 | Floor area of admin building in m^2 |
a_site_service_bldg | Input | real | 8300.0 | Floor area of site service building in m^2 |
a_cryo_inert_gas_bldg | Input | real | 18380.0 | Floor area of cryogenics and inert gas storage building in m^2 |
a_security_bldg | Input | real | 4552.0 | Floor area of security building in m^2 |
times_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
pulsetimings | Input | real | 1.0 | Switch for pulse timings (if lpulse=1):
|
tburn | Input | real | 1000.0 | burn time (s) (calculated if |
tburn0 | Output | real | - | burn time (s) - used for internal consistency |
tcycle | Output | real | - | full cycle time (s) |
tdown | Output | real | - | down time (s) |
tdwell | Input | real | 1800.0 | time between pulses in a pulsed reactor (s) ( |
t_fusion_ramp | Input | real | 10.0 | heating time, after current ramp up (s) |
tim | Output | real | - | array of time points during plasma pulse (s) |
timelabel | Input | character | [b'Start ' b'BOP ' b'EOR ' b'BOF ' b'EOF ' b'EOP '] | array of time labels during plasma pulse (s) |
intervallabel | Input | character | [b'tramp ' b'tohs ' b't_fusion_ra' b'tburn ' b'tqnch '] | time intervals - as strings (s) |
tohs | Input | real | 30.0 | plasma current ramp-up time for current initiation (s) (calculated if |
tohsin | Output | real | - | Switch for plasma current ramp-up time (if lpulse=0):
|
tpulse | Output | real | - | pulse length = tohs + t_fusion_ramp + tburn + tqnch |
tqnch | Input | real | 15.0 | shut down time for PF coils (s); if pulsed, = tohs |
tramp | Input | real | 15.0 | initial PF coil charge time (s); if pulsed, = tohs |
physics_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
ipnlaws | Parameter | integer | 50 | number of energy confinement time scaling laws |
abeam | Output | real | - | beam ion mass (amu) |
afuel | Output | real | - | average mass of fuel portion of ions (amu) |
aion | Output | real | - | average mass of all ions (amu) |
alphaj | Input | real | 1.0 | current profile index (calculated from q_0 and q if |
alphan | Input | real | 0.25 | density profile index |
alphap | Output | real | - | pressure profile index |
alpharate | Output | real | - | alpha particle production rate (particles/m3/sec) |
alphat | Input | real | 0.5 | temperature profile index |
aspect | Input | real | 2.907 | aspect ratio ( |
beamfus0 | Input | real | 1.0 | multiplier for beam-background fusion calculation |
beta | Input | real | 0.042 | total plasma beta ( |
betaft | Output | real | - | fast alpha beta component |
betalim | Output | real | - | allowable beta |
betalim_lower | Output | real | - | allowable lower beta |
betanb | Output | real | - | neutral beam beta component |
betap | Output | real | - | poloidal beta |
normalised_total_beta | Output | real | - | normaised total beta |
betbm0 | Input | real | 1.5 | leading coefficient for NB beta fraction |
bp | Output | real | - | poloidal field (T) |
bt | Input | real | 5.68 | toroidal field on axis (T) ( |
btot | Output | real | - | total toroidal + poloidal field (T) |
burnup | Output | real | - | fractional plasma burnup |
burnup_in | Output | real | - | fractional plasma burnup user input |
bvert | Output | real | - | vertical field at plasma (T) |
c_beta | Input | real | 0.5 | Destabalisation parameter for iprofile=6 beta limit |
csawth | Input | real | 1.0 | coeff. for sawteeth effects on burn V-s requirement |
cvol | Input | real | 1.0 | multiplying factor times plasma volume (normally=1) |
cwrmax | Input | real | 1.35 | maximum ratio of conducting wall distance to plasma minor radius for
vertical stability ( |
dene | Input | real | 9.8e+19 | electron density (/m3) ( |
deni | Output | real | - | fuel ion density (/m3) |
dlamee | Output | real | - | electron-electron coulomb logarithm |
dlamie | Output | real | - | ion-electron coulomb logarithm |
dlimit | Output | real | - | density limit (/m3) as calculated using various models |
dnalp | Output | real | - | thermal alpha density (/m3) |
dnbeam | Output | real | - | hot beam ion density, variable (/m3) |
dnbeam2 | Output | real | - | hot beam ion density from calculation (/m3) |
dnbeta | Input | real | 3.5 | Troyon-like coefficient for beta scaling |
dnelimt | Output | real | - | density limit (/m3) |
dnitot | Output | real | - | total ion density (/m3) |
dnla | Output | real | - | line averaged electron density (/m3) |
dnprot | Output | real | - | proton ash density (/m3) |
dntau | Output | real | - | plasma average "n-tau" (seconds/m3) |
dnz | Output | real | - | high Z ion density (/m3) |
gradient_length_ne | Output | real | - | Max. normalized gradient length in el. density (ipedestal==0 only) |
gradient_length_te | Output | real | - | Max. normalized gradient length in el. temperature (ipedestal==0 only) |
ealphadt | Parameter | real | 3520.0D0 | alpha birth energy in D-T reaction (keV) |
epbetmax | Input | real | 1.38 | maximum (eps*beta_poloidal) ( |
eps | Input | real | 0.34399724802 | inverse aspect ratio |
aux_current_fraction | Output | real | - | fraction of plasma current produced by auxiliary current drive |
inductive_current_fraction | Output | real | - | fraction of plasma current produced inductively |
falpe | Output | real | - | fraction of alpha energy to electrons |
falpha | Input | real | 0.95 | fraction of alpha power deposited in plasma (Physics of Energetic Ions, p.2489) |
falpi | Output | real | - | fraction of alpha power to ions |
fdeut | Input | real | 0.5 | deuterium fuel fraction |
ftar | Input | real | 1.0 | fraction of power to the lower divertor in double null configuration
( |
ffwal | Input | real | 0.92 | factor to convert plasma surface area to first wall area in neutron wall
load calculation ( |
fgwped | Input | real | 0.85 | fraction of Greenwald density to set as pedestal-top density. If |
fgwsep | Input | real | 0.5 | fraction of Greenwald density to set as separatrix density. If |
fhe3 | Output | real | - | helium-3 fuel fraction |
figmer | Output | real | - | physics figure of merit (= plasma_currentaspect*sbar, where |
fkzohm | Input | real | 1.0 | Zohm elongation scaling adjustment factor ( |
fplhsep | Input | real | 1.0 | F-value for Psep >= Plh + Paux ( |
fpdivlim | Input | real | 1.0 | F-value for minimum pdivt ( |
fne0 | Input | real | 1.0 | f-value for the constraint ne(0) > ne(ped) ( |
ftrit | Input | real | 0.5 | tritium fuel fraction |
fusionrate | Output | real | - | fusion reaction rate (reactions/m3/sec) |
fvsbrnni | Input | real | 1.0 | fraction of the plasma current produced by non-inductive means ( |
gamma | Input | real | 0.4 | Ejima coefficient for resistive startup V-s formula |
gammaft | Output | real | - | ratio of (fast alpha + neutral beam beta) to thermal beta |
hfac | Output | real | - | H factors for an ignited plasma for each energy confinement time scaling law |
hfact | Input | real | 1.0 | H factor on energy confinement times, radiation corrected ( |
taumax | Input | real | 10.0 | Maximum allowed energy confinement time (s) |
i_bootstrap_current | Input | integer | 3 | switch for bootstrap current scaling
|
iculbl | Output | integer | - | switch for beta limit scaling (
|
i_plasma_current | Input | integer | 4 | switch for plasma current scaling to use
|
i_diamagnetic_current | Output | integer | - | switch for diamagnetic current scaling
|
idensl | Input | integer | 7 | switch for density limit to enforce (
|
idivrt | Input | integer | 2 | number of divertors (calculated from |
ifalphap | Input | integer | 1 | switch for fast alpha pressure calculation
|
ignite | Output | integer | - | switch for ignition assumption. Obviously, ignite must be zero if current drive is required. If ignite 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.
|
iinvqd | Input | integer | 1 | switch for inverse quadrature in L-mode scaling laws 5 and 9:
|
ipedestal | Input | integer | 1 | switch for pedestal profiles:
|
i_pfirsch_schluter_current | Output | integer | - | switch for Pfirsch-Schlüter current scaling (issue #413):
|
neped | Input | real | 4e+19 | electron density of pedestal [m-3] (`ipedestal==1) |
nesep | Input | real | 3e+19 | electron density at separatrix [m-3] (`ipedestal==1) |
alpha_crit | Output | real | - | critical ballooning parameter value |
nesep_crit | Output | real | - | critical electron density at separatrix [m-3] |
plasma_res_factor | Input | real | 1.0 | plasma resistivity pre-factor |
rhopedn | Input | real | 1.0 | r/a of density pedestal ( |
rhopedt | Input | real | 1.0 | r/a of temperature pedestal ( |
rho_te_max | Output | real | - | r/a where the temperature gradient is largest ( |
rho_ne_max | Output | real | - | r/a where the density gradient is largest ( |
tbeta | Input | real | 2.0 | temperature profile index beta (`ipedestal==1) |
teped | Input | real | 1.0 | electron temperature of pedestal (keV) ( |
tesep | Input | real | 0.1 | electron temperature at separatrix (keV) ( |
iprofile | Input | integer | 1 | switch for current profile consistency:
|
iradloss | Input | integer | 1 | switch for radiation loss term usage in power balance (see User Guide):
|
isc | Input | integer | 34 | switch for energy confinement time scaling law (see description in |
tauscl | Parameter | character | (/'Neo-Alcator (ohmic)', 'Mirnov (H)', 'Merezkhin-Muhkovatov (L)', 'Shimomura (H)', 'Kaye-Goldston (L)', 'ITER 89-P (L)', 'ITER 89-O (L)', 'Rebut-Lallia (L)', 'Goldston (L)', 'T10 (L)', 'JAERI-88 (L)', 'Kaye-Big Complex (L)', 'ITER H90-P (H)', 'ITER Mix (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 (H)', 'TITAN RFP OBSOLETE ', 'ITER H-97P ELM-free (H)', 'ITER H-97P ELMy (H)', 'ITER-96P (L)', 'Valovic modified ELMy(H)', 'Kaye PPPL April 98 (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 (H)', 'Murari et al NPL (H)', 'Petty 2008 (H)', 'Lang et al. 2012 (H)', 'Hubbard 2017 - nom (I)', 'Hubbard 2017 - lower (I)', 'Hubbard 2017 - upper (I)', 'NSTX (Spherical) (H)', 'NSTX-Petty08 Hybrid (H)', 'NSTX gyro-Bohm Buxton(H)', 'Input tauee_in ', 'ITPA20 (H)'/) | tauscl(ipnlaws) : labels describing energy confinement scaling laws:
|
iscrp | Input | integer | 1 | switch for plasma-first wall clearances:
|
ishape | Output | integer | - | switch for plasma cross-sectional shape calculation:
|
itart | Output | integer | - | switch for spherical tokamak (ST) models:
|
itartpf | Output | integer | - | switch for Spherical Tokamak PF models:
|
iwalld | Input | integer | 1 | switch for neutron wall load calculation:
|
kappa | Input | real | 1.792 | plasma separatrix elongation (calculated if |
kappa95 | Input | real | 1.6 | plasma elongation at 95% surface (calculated if |
kappaa | Output | real | - | plasma elongation calculated as xarea/(pi.a^2) |
kappaa_ipb | Output | real | - | Volume measure of plasma elongation |
ne0 | Output | real | - | central electron density (/m3) |
ni0 | Output | real | - | central ion density (/m3) |
m_s_limit | Input | real | 0.3 | margin to vertical stability |
p0 | Output | real | - | central total plasma pressure (Pa) |
palppv | Output | real | - | alpha power per volume (MW/m3) |
palpepv | Output | real | - | alpha power per volume to electrons (MW/m3) |
palpfwmw | Output | real | - | alpha power escaping plasma and reaching first wall (MW) |
palpipv | Output | real | - | alpha power per volume to ions (MW/m3) |
palpmw | Output | real | - | alpha power (MW) |
palpnb | Output | real | - | alpha power from hot neutral beam ions (MW) |
pchargemw | Output | real | - | non-alpha charged particle fusion power (MW) |
pchargepv | Output | real | - | non-alpha charged particle fusion power per volume (MW/m3) |
pcoef | Output | real | - | profile factor (= n-weighted T / average T) |
pinnerzoneradmw | Output | real | - | radiation power from inner zone (MW) |
pcoreradpv | Output | real | - | total core radiation power per volume (MW/m3) |
pdd | Output | real | - | deuterium-deuterium fusion power (MW) |
pdhe3 | Output | real | - | deuterium-helium3 fusion power (MW) |
pdivt | Output | real | - | power to conducted to the divertor region (MW) |
pdivl | Output | real | - | power conducted to the lower divertor region (calculated if |
pdivu | Output | real | - | power conducted to the upper divertor region (calculated if |
pdivmax | Output | real | - | power conducted to the divertor with most load (calculated if |
pdt | Output | real | - | deuterium-tritium fusion power (MW) |
pouterzoneradmw | Output | real | - | radiation power from outer zone (MW) |
pedgeradpv | Output | real | - | edge radiation power per volume (MW/m3) |
pfuscmw | Output | real | - | charged particle fusion power (MW) |
phiint | Output | real | - | internal plasma V-s |
photon_wall | Output | real | - | Nominal mean radiation load on inside surface of reactor (MW/m2) |
piepv | Output | real | - | ion/electron equilibration power per volume (MW/m3) |
plasma_current | Output | real | - | plasma current (A) |
pneutmw | Output | real | - | neutron fusion power (MW) |
pneutpv | Output | real | - | neutron fusion power per volume (MW/m3) |
pohmmw | Output | real | - | ohmic heating power (MW) |
pohmpv | Output | real | - | ohmic heating power per volume (MW/m3) |
powerht | Output | real | - | heating power (= transport loss power) (MW) used in confinement time calculation |
powfmw | Output | real | - | fusion power (MW) |
pperim | Output | real | - | plasma poloidal perimeter (m) |
pradmw | Output | real | - | total radiation power from inside LCFS (MW) |
pradpv | Output | real | - | total radiation power per volume (MW/m3) |
pradsolmw | Output | real | - | radiation power from SoL (MW) |
protonrate | Output | real | - | proton production rate (particles/m3/sec) |
psolradmw | Output | real | - | SOL radiation power (MW) ( |
psyncpv | Output | real | - | synchrotron radiation power per volume (MW/m3) |
ilhthresh | Input | integer | 19 | switch for L-H mode power threshold scaling to use (see pthrmw for list) |
plhthresh | Output | real | - | L-H mode power threshold (MW) (chosen via ilhthresh, and enforced if constraint equation 15 is on) |
pthrmw | Output | real | - | L-H power threshold for various scalings (MW)
|
ptremw | Output | real | - | electron transport power (MW) |
ptrepv | Output | real | - | electron transport power per volume (MW/m3) |
ptrimw | Output | real | - | ion transport power (MW) |
pscalingmw | Output | real | - | Total transport power from scaling law (MW) |
ptripv | Output | real | - | ion transport power per volume (MW/m3) |
q | Input | real | 3.0 | Safety factor 'near' plasma edge ( |
q0 | Input | real | 1.0 | safety factor on axis |
q95 | Output | real | - | safety factor at 95% surface |
qfuel | Output | real | - | plasma fuelling rate (nucleus-pairs/s) |
tauratio | Input | real | 1.0 | tauratio /1.0/ : ratio of He and pellet particle confinement times |
qlim | Output | real | - | lower limit for edge safety factor |
qstar | Output | real | - | cylindrical safety factor |
rad_fraction_sol | Input | real | 0.8 | SoL radiation fraction |
rad_fraction_total | Output | real | - | Radiation fraction total = SoL + LCFS radiation / total power deposited in plasma |
ralpne | Input | real | 0.1 | thermal alpha density/electron density ( |
protium | Output | real | - | Seeded protium density / electron density. |
rli | Input | real | 0.9 | plasma normalised internal inductance (calculated from alphaj if |
rlp | Output | real | - | plasma inductance (H) |
rmajor | Input | real | 8.14 | plasma major radius (m) ( |
rminor | Output | real | - | plasma minor radius (m) |
rnbeam | Input | real | 0.005 | hot beam density / n_e ( |
rncne | Output | real | - | n_carbon / n_e |
rndfuel | Output | real | - | fuel burnup rate (reactions/second) |
rnfene | Output | real | - | n_highZ / n_e |
rnone | Output | real | - | n_oxygen / n_e |
rpfac | Output | real | - | neo-classical correction factor to rplas |
rplas | Output | real | - | plasma resistance (ohm) |
res_time | Output | real | - | plasma current resistive diffusion time (s) |
sarea | Output | real | - | plasma surface area |
sareao | Output | real | - | outboard plasma surface area |
sf | Output | real | - | shape factor = plasma poloidal perimeter / (2.pi.rminor) |
i_single_null | Input | integer | 1 | switch for single null / double null plasma:
|
ssync | Input | real | 0.6 | synchrotron wall reflectivity factor |
tauee | Output | real | - | electron energy confinement time (sec) |
tauee_in | Output | real | - | Input electron energy confinement time (sec) ( |
taueff | Output | real | - | global thermal energy confinement time (sec) |
tauei | Output | real | - | ion energy confinement time (sec) |
taup | Output | real | - | alpha particle confinement time (sec) |
te | Input | real | 12.9 | volume averaged electron temperature (keV) ( |
te0 | Output | real | - | central electron temperature (keV) |
ten | Output | real | - | density weighted average electron temperature (keV) |
ti | Input | real | 12.9 | volume averaged ion temperature (keV). N.B. calculated from te if |
ti0 | Output | real | - | central ion temperature (keV) |
tin | Output | real | - | density weighted average ion temperature (keV) |
tratio | Input | real | 1.0 | ion temperature / electron temperature(used to calculate ti if |
triang | Input | real | 0.36 | plasma separatrix triangularity (calculated if |
triang95 | Input | real | 0.24 | plasma triangularity at 95% surface (calculated if |
vol | Output | real | - | plasma volume (m3) |
vsbrn | Output | real | - | V-s needed during flat-top (heat + burn times) (Wb) |
vshift | Output | real | - | plasma/device midplane vertical shift - single null |
vsind | Output | real | - | internal and external plasma inductance V-s (Wb) |
vsres | Output | real | - | resistive losses in startup V-s (Wb) |
vsstt | Output | real | - | total V-s needed (Wb) |
wallmw | Output | real | - | average neutron wall load (MW/m2) |
wtgpd | Output | real | - | mass of fuel used per day (g) |
xarea | Output | real | - | plasma cross-sectional area (m2) |
zeff | Output | real | - | plasma effective charge |
zeffai | Output | real | - | mass weighted plasma effective charge |
cost_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
abktflnc | Input | real | 5.0 | allowable first wall/blanket neutron fluence (MW-yr/m2) ( |
adivflnc | Input | real | 7.0 | allowable divertor heat fluence (MW-yr/m2) |
blkcst | Output | real | - | blanket direct cost (M$) |
c221 | Output | real | - | total account 221 cost (M$) - first wall, blanket, shield, support structure and div plates |
c222 | Output | real | - | total account 222 cost (M$) - TF coils + PF coils |
capcost | Output | real | - | total capital cost including interest (M$) |
cconfix | Input | real | 80.0 | fixed cost of superconducting cable ($/m) |
cconshpf | Input | real | 70.0 | cost of PF coil steel conduit/sheath ($/m) |
cconshtf | Input | real | 75.0 | cost of TF coil steel conduit/sheath ($/m) |
cdcost | Output | real | - | current drive direct costs (M$) |
cdirt | Output | real | - | total plant direct cost (M$) |
cdrlife | Output | real | - | Full power year lifetime of heating/current drive system (y) |
cdrlife_cal | Output | real | - | Calendar year lifetime of heating/current drive system (y) |
cfactr | Input | real | 0.75 | Total plant availability fraction; input if |
cpfact | Output | real | - | Total plant capacity factor |
cfind | Input | real | [0.244 0.244 0.244 0.29 ] | indirect cost factor (func of lsa) (cost model = 0) |
cland | Input | real | 19.2 | cost of land (M$) |
coe | Output | real | - | cost of electricity ($/MW-hr) |
coecap | Output | real | - | capital cost of electricity (m$/kW-hr) |
coefuelt | Output | real | - | 'fuel' (including replaceable components) contribution to cost of electricity (m$/kW-hr) |
coeoam | Output | real | - | operation and maintenance contribution to cost of electricity (m$/kW-hr) |
concost | Output | real | - | plant construction cost (M$) |
costexp | Input | real | 0.8 | cost exponent for scaling in 2015 costs model |
costexp_pebbles | Input | real | 0.6 | cost exponent for pebbles in 2015 costs model |
cost_factor_buildings | Input | real | 1.0 | cost scaling factor for buildings |
cost_factor_land | Input | real | 1.0 | cost scaling factor for land |
cost_factor_tf_coils | Input | real | 1.0 | cost scaling factor for TF coils |
cost_factor_fwbs | Input | real | 1.0 | cost scaling factor for fwbs |
cost_factor_rh | Input | real | 1.0 | cost scaling factor for remote handling |
cost_factor_vv | Input | real | 1.0 | cost scaling factor for vacuum vessel |
cost_factor_bop | Input | real | 1.0 | cost scaling factor for energy conversion system |
cost_factor_misc | Input | real | 1.0 | cost scaling factor for remaining subsystems |
maintenance_fwbs | Input | real | 0.2 | Maintenance cost factor: first wall, blanket, shield, divertor |
maintenance_gen | Input | real | 0.05 | Maintenance cost factor: All other components except coils, vacuum vessel, thermal shield, cryostat, land |
amortization | Input | real | 13.6 | amortization factor (fixed charge factor) "A" (years) |
cost_model | Input | integer | 1 | Switch for cost model:
|
i_cp_lifetime | Output | integer | - | Switch for the centrepost lifetime constraint 0 : The CP full power year lifetime is set by the user via cplife_input 1 : The CP lifetime is equal to the divertor lifetime 2 : The CP lifetime is equal to the breeding blankets lifetime 3 : The CP lifetime is equal to the plant lifetime |
cowner | Input | real | 0.15 | owner cost factor |
cplife_input | Input | real | 2.0 | User input full power year lifetime of the centrepost (years) (i_cp_lifetime = 0) |
cplife | Output | real | - | Calculated full power year lifetime of centrepost (years) |
cplife_cal | Output | real | - | Calculated calendar year lifetime of centrepost (years) |
cpstcst | Output | real | - | ST centrepost direct cost (M$) |
cpstflnc | Input | real | 10.0 | allowable ST centrepost neutron fluence (MW-yr/m2) |
crctcore | Output | real | - | reactor core costs (categories 221, 222 and 223) |
csi | Input | real | 16.0 | allowance for site costs (M$) |
cturbb | Input | real | 38.0 | cost of turbine building (M$) |
decomf | Input | real | 0.1 | proportion of constructed cost required for decommissioning fund |
dintrt | Output | real | - | diff between borrowing and saving interest rates |
divcst | Output | real | - | divertor direct cost (M$) |
divlife | Output | real | - | Full power lifetime of divertor (y) |
divlife_cal | Output | real | - | Calendar year lifetime of divertor (y) |
dtlife | Output | real | - | period prior to the end of the plant life that the decommissioning fund is used (years) |
fcap0 | Input | real | 1.165 | average cost of money for construction of plant assuming design/construction time of six years |
fcap0cp | Input | real | 1.08 | average cost of money for replaceable components assuming lead time for these of two years |
fcdfuel | Input | real | 0.1 | fraction of current drive cost treated as fuel (if |
fcontng | Input | real | 0.195 | project contingency factor |
fcr0 | Input | real | 0.0966 | fixed charge rate during construction |
fkind | Input | real | 1.0 | multiplier for Nth of a kind costs |
fwallcst | Output | real | - | first wall cost (M$) |
iavail | Input | integer | 2 | Switch for plant availability model:
|
ibkt_life | Output | integer | - | Switch for fw/blanket lifetime calculation in availability module:
|
life_dpa | Input | real | 50.0 | Allowable DPA from DEMO fw/blanket lifetime calculation in availability module |
bktcycles | Input | real | 1000.0 | Number of fusion cycles to reach allowable DPA from DEMO fw/blanket lifetime calculation |
avail_min | Input | real | 0.75 | Minimum availability ( |
tok_build_cost_per_vol | Input | real | 1283.0 | Unit cost for tokamak complex buildings, including building and site services ($/m3) |
light_build_cost_per_vol | Input | real | 270.0 | Unit cost for unshielded non-active buildings ($/m3) |
favail | Input | real | 1.0 | F-value for minimum availability ( |
num_rh_systems | Input | integer | 4 | Number of remote handling systems (1-10) |
conf_mag | Input | real | 0.99 | c parameter, which determines the temperature margin at which magnet lifetime starts to decline |
div_prob_fail | Input | real | 0.0002 | Divertor probability of failure (per op day) |
div_umain_time | Input | real | 0.25 | Divertor unplanned maintenance time (years) |
div_nref | Input | real | 7000.0 | Reference value for cycle cycle life of divertor |
div_nu | Input | real | 14000.0 | The cycle when the divertor fails with 100% probability |
fwbs_nref | Input | real | 20000.0 | Reference value for cycle life of blanket |
fwbs_nu | Input | real | 40000.0 | The cycle when the blanket fails with 100% probability |
fwbs_prob_fail | Input | real | 0.0002 | Fwbs probability of failure (per op day) |
fwbs_umain_time | Input | real | 0.25 | Fwbs unplanned maintenance time (years) |
redun_vacp | Input | real | 25.0 | Vacuum system pump redundancy level (%) |
redun_vac | Output | integer | - | Number of redundant vacuum pumps |
t_operation | Output | real | - | Operational time (yrs) |
tbktrepl | Input | real | 0.5 | time taken to replace blanket (y) ( |
tcomrepl | Input | real | 0.5 | time taken to replace both blanket and divertor (y) ( |
tdivrepl | Input | real | 0.25 | time taken to replace divertor (y) ( |
uubop | Input | real | 0.02 | unplanned unavailability factor for balance of plant ( |
uucd | Input | real | 0.02 | unplanned unavailability factor for current drive ( |
uudiv | Input | real | 0.04 | unplanned unavailability factor for divertor ( |
uufuel | Input | real | 0.02 | unplanned unavailability factor for fuel system ( |
uufw | Input | real | 0.04 | unplanned unavailability factor for first wall ( |
uumag | Input | real | 0.02 | unplanned unavailability factor for magnets ( |
uuves | Input | real | 0.04 | unplanned unavailability factor for vessel ( |
ifueltyp | Output | integer | - | Switch for fuel type:
|
ipnet | Output | integer | - | Switch for net electric power calculation:
|
ireactor | Input | integer | 1 | Switch for net electric power and cost of electricity calculations:
|
lsa | Input | integer | 4 | Level of safety assurance switch (generally, use 3 or 4):
|
moneyint | Output | real | - | interest portion of capital cost (M$) |
output_costs | Input | integer | 1 | Switch for costs output:
|
discount_rate | Input | real | 0.0435 | effective cost of money in constant dollars |
startupratio | Input | real | 1.0 | ratio of additional HCD power for start-up to flat-top operational requirements |
startuppwr | Output | real | - | cost associated with additional HCD system power required on start-up ($) |
supercond_cost_model | Output | integer | - | Switch for superconductor cost model:
|
tlife | Input | real | 30.0 | Full power year plant lifetime (years) |
tmain | Output | real | - | Maintenance time for replacing CP (years) (iavail = 3) |
u_unplanned_cp | Output | real | - | User-input CP unplanned unavailability (iavail = 3) |
ucad | Parameter | real | 180.0D0 | unit cost for administration buildings (M$/m3) |
ucaf | Parameter | real | 1.5D6 | unit cost for aux facility power equipment ($) |
ucahts | Parameter | real | 31.0D0 | unit cost for aux heat transport equipment ($/W**exphts) |
ucap | Parameter | real | 17.0D0 | unit cost of auxiliary transformer ($/kVA) |
ucblbe | Input | real | 260.0 | unit cost for blanket beryllium ($/kg) |
ucblbreed | Input | real | 875.0 | unit cost for breeder material ($/kg) ( |
ucblli | Input | real | 875.0 | unit cost for blanket lithium ($/kg) (30% Li6) |
ucblli2o | Input | real | 600.0 | unit cost for blanket Li_2O ($/kg) |
ucbllipb | Input | real | 10.3 | unit cost for blanket Li-Pb ($/kg) (30% Li6) |
ucblss | Input | real | 90.0 | unit cost for blanket stainless steel ($/kg) |
ucblvd | Input | real | 200.0 | unit cost for blanket vanadium ($/kg) |
ucbpmp | Parameter | real | 2.925D5 | vacuum system backing pump cost ($) |
ucbus | Input | real | 0.123 | cost of aluminium bus for TF coil ($/A-m) |
uccase | Input | real | 50.0 | cost of superconductor case ($/kg) |
ucco | Parameter | real | 350.0D0 | unit cost for control buildings (M$/m3) |
uccpcl1 | Input | real | 250.0 | cost of high strength tapered copper ($/kg) |
uccpclb | Input | real | 150.0 | cost of TF outboard leg plate coils ($/kg) |
uccpmp | Parameter | real | 3.9D5 | vacuum system cryopump cost ($) |
uccr | Parameter | real | 460.0D0 | unit cost for cryogenic building (M$/vol) |
uccry | Input | real | 93000.0 | heat transport system cryoplant costs ($/W**expcry) |
uccryo | Input | real | 32.0 | unit cost for vacuum vessel ($/kg) |
uccu | Input | real | 75.0 | unit cost for copper in superconducting cable ($/kg) |
ucdgen | Parameter | real | 1.7D6 | cost per 8 MW diesel generator ($) |
ucdiv | Input | real | 280000.0 | cost of divertor blade ($) |
ucdtc | Parameter | real | 13.0D0 | detritiation, air cleanup cost ($/10000m3/hr) |
ucduct | Parameter | real | 4.225D4 | vacuum system duct cost ($/m) |
ucech | Input | real | 3.0 | ECH system cost ($/W) |
ucel | Parameter | real | 380.0D0 | unit cost for electrical equipment building (M$/m3) |
uces1 | Parameter | real | 3.2D4 | MGF (motor-generator flywheel) cost factor ($/MVA**0.8) |
uces2 | Parameter | real | 8.8D3 | MGF (motor-generator flywheel) cost factor ($/MJ**0.8) |
ucf1 | Input | real | 22300000.0 | cost of fuelling system ($) |
ucfnc | Input | real | 35.0 | outer PF coil fence support cost ($/kg) |
ucfpr | Parameter | real | 4.4D7 | cost of 60g/day tritium processing unit ($) |
ucfuel | Input | real | 3.45 | unit cost of D-T fuel (M$/year/1200MW) |
ucfwa | Parameter | real | 6.0D4 | first wall armour cost ($/m2) |
ucfwps | Parameter | real | 1.0D7 | first wall passive stabiliser cost ($) |
ucfws | Parameter | real | 5.3D4 | first wall structure cost ($/m2) |
ucgss | Parameter | real | 35.0D0 | cost of reactor structure ($/kg) |
uche3 | Input | real | 1000000.0 | cost of helium-3 ($/kg) |
uchrs | Input | real | 87900000.0 | cost of heat rejection system ($) |
uchts | Input | real | [15.3 19.1] | cost of heat transport system equipment per loop ($/W); dependent on coolant type (coolwh) |
uciac | Input | real | 150000000.0 | cost of instrumentation, control & diagnostics ($) |
ucich | Input | real | 3.0 | ICH system cost ($/W) |
ucint | Parameter | real | 35.0D0 | superconductor intercoil structure cost ($/kg) |
uclh | Input | real | 3.3 | lower hybrid system cost ($/W) |
uclv | Parameter | real | 16.0D0 | low voltage system cost ($/kVA) |
ucmb | Parameter | real | 260.0D0 | unit cost for reactor maintenance building (M$/m3) |
ucme | Input | real | 125000000.0 | cost of maintenance equipment ($) |
ucmisc | Input | real | 25000000.0 | miscellaneous plant allowance ($) |
ucnbi | Input | real | 3.3 | NBI system cost ($/W) |
ucnbv | Parameter | real | 1000.0D0 | cost of nuclear building ventilation ($/m3) |
ucoam | Input | real | [68.8 68.8 68.8 74.4] | annual cost of operation and maintenance (M$/year/1200MW**0.5) |
ucpens | Input | real | 32.0 | penetration shield cost ($/kg) |
ucpfb | Input | real | 210.0 | cost of PF coil buses ($/kA-m) |
ucpfbk | Input | real | 16600.0 | cost of PF coil DC breakers ($/MVA**0.7) |
ucpfbs | Input | real | 4900.0 | cost of PF burn power supplies ($/kW**0.7) |
ucpfcb | Input | real | 75000.0 | cost of PF coil AC breakers ($/circuit) |
ucpfdr1 | Input | real | 150.0 | cost factor for dump resistors ($/MJ) |
ucpfic | Input | real | 10000.0 | cost of PF instrumentation and control ($/channel) |
ucpfps | Input | real | 35000.0 | cost of PF coil pulsed power supplies ($/MVA) |
ucphx | Parameter | real | 15.0D0 | primary heat transport cost ($/W**exphts) |
ucpp | Parameter | real | 48.0D0 | cost of primary power transformers ($/kVA**0.9) |
ucrb | Input | real | 400.0 | cost of reactor building (M$/m3) |
ucsc | Input | real | [ 600. 600. 300. 600. 600. 600. 300. 1200. 1200.] | cost of superconductor ($/kg) |
sc_mat_cost_0 | Input | real | [ 4.8 2. 1. 4.8 4.8 47.4 1. 47.4 47.4] | cost of superconductor ($/kA m) at 6.4 T, 4.2 K |
ucsh | Parameter | real | 115.0D0 | cost of shops and warehouses (M$/m3) |
ucshld | Input | real | 32.0 | cost of shield structural steel ($/kg) |
ucswyd | Parameter | real | 1.84D7 | switchyard equipment costs ($) |
uctfbr | Input | real | 1.22 | cost of TF coil breakers ($/W**0.7) |
uctfbus | Input | real | 100.0 | cost of TF coil bus ($/kg) |
uctfdr | Parameter | real | 1.75D-4 | cost of TF coil dump resistors ($/J) |
uctfgr | Parameter | real | 5000.0D0 | additional cost of TF coil dump resistors ($/coil) |
uctfic | Parameter | real | 1.0D4 | cost of TF coil instrumentation and control ($/coil/30) |
uctfps | Input | real | 24.0 | cost of TF coil power supplies ($/W**0.7) |
uctfsw | Input | real | 1.0 | cost of TF coil slow dump switches ($/A) |
uctpmp | Parameter | real | 1.105D5 | cost of turbomolecular pump ($) |
uctr | Parameter | real | 370.0D0 | cost of tritium building ($/m3) |
ucturb | Input | real | [2.30e+08 2.45e+08] | cost of turbine plant equipment ($) (dependent on coolant type coolwh) |
ucvalv | Parameter | real | 3.9D5 | vacuum system valve cost ($) |
ucvdsh | Parameter | real | 26.0D0 | vacuum duct shield cost ($/kg) |
ucviac | Parameter | real | 1.3D6 | vacuum system instrumentation and control cost ($) |
ucwindpf | Input | real | 465.0 | cost of PF coil superconductor windings ($/m) |
ucwindtf | Input | real | 480.0 | cost of TF coil superconductor windings ($/m) |
ucws | Parameter | real | 460.0D0 | cost of active assembly shop ($/m3) |
ucwst | Input | real | [0. 3.94 5.91 7.88] | cost of waste disposal (M$/y/1200MW) |
process_input
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
nin | Parameter | integer | 10 | |
maxlen* | Parameter | integer | 2000 | |
line* | Variable | character | - | |
linelen* | Variable | integer | - | |
lineno* | Variable | integer | - | |
iptr* | Variable | integer | - | |
infile* | Variable | integer | - | |
outfile* | Variable | integer | - | |
report_changes* | Variable | integer | - | |
icode* | Variable | integer | - | |
subscript_present* | Variable | logical | - | |
error* | Variable | logical | - | |
error_message* | Variable | character | - | |
show_changes* | Variable | integer | - | |
constraints_exist* | Variable | logical | - |
reinke_module
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
vcritx | Output | real | - |
error_handling
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
errors_on | Output | logical | - | |
error_okay* | Parameter | integer | 0 | |
error_info* | Parameter | integer | 1 | |
error_warn* | Parameter | integer | 2 | |
error_severe* | Parameter | integer | 3 | |
error_id* | Variable | integer | - | error_id : identifier for final message encountered |
error_status | Output | integer | - | error_status : overall status flag for a run; on exit:
|
int_default* | Parameter | integer | -999999 | |
flt_default* | Parameter | real | real(INT_DEFAULT, kind(1.0D0)) | |
idiags | Input | integer | [-999999 -999999 -999999 -999999 -999999 -999999 -999999 -999999] | |
fdiags | Input | real | [-999999. -999999. -999999. -999999. -999999. -999999. -999999. -999999.] | |
error_head* | Variable | type | - | |
error_tail* | Variable | type | - | |
error_type* | Variable | type | - |
rebco_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
rebco_thickness | Input | real | 1e-06 | thickness of REBCO layer in tape (m) ( |
copper_thick | Input | real | 0.0001 | thickness of copper layer in tape (m) ( |
hastelloy_thickness | Input | real | 5e-05 | thickness of Hastelloy layer in tape (m) |
tape_width | Input | real | 0.004 | Mean width of tape (m) |
tape_thickness | Input | real | 6.5e-05 | thickness of tape, inc. all layers (hts, copper, substrate, etc.) (m) |
croco_od | Output | real | - | Outer diameter of CroCo strand (m) |
croco_id | Output | real | - | Inner diameter of CroCo copper tube (m) |
croco_thick | Input | real | 0.0025 | Thickness of CroCo copper tube (m) ( |
copper_rrr | Input | real | 100.0 | residual resistivity ratio copper in TF superconducting cable |
coppera_m2 | Output | real | - | TF coil current / copper area (A/m2) |
coppera_m2_max | Input | real | 100000000.0 | Maximum TF coil current / copper area (A/m2) |
f_coppera_m2 | Input | real | 1.0 | f-value for constraint 75: TF coil current / copper area < copperA_m2_max |
copperaoh_m2 | Output | real | - | CS coil current / copper area (A/m2) ( |
copperaoh_m2_max | Input | real | 100000000.0 | Maximum CS coil current / copper area (A/m2) |
f_copperaoh_m2 | Input | real | 1.0 | f-value for constraint 88: CS coil current / copper area < copperA_m2_max |
stack_thickness | Output | real | - | |
tapes | Output | real | - | |
rebco_area | Output | real | - | |
copper_area | Output | real | - | |
hastelloy_area | Output | real | - | |
solder_area | Output | real | - | |
croco_area | Output | real | - |
reinke_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
impvardiv | Input | integer | 9 | Index of impurity to be iterated for Reinke divertor detachment criterion |
lhat | Input | real | 4.33 | Connection length factor L|| = lhat qstar R for Reinke criterion, default value from Post et al. 1995 J. Nucl. Mat. 220-2 1014 |
fzmin | Output | real | - | Minimum impurity fraction necessary for detachment. This is the impurity at the SOL/Div. |
fzactual | Input | real | 0.001 | Actual impurity fraction of divertor impurity (impvardiv) in the SoL (taking
impurity_enrichment into account) ( |
reinke_mode | Output | integer | - | Switch for Reinke criterion H/I mode:
|
water_usage_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
airtemp | Input | real | 15.0 | ambient air temperature (degrees Celsius) |
watertemp | Input | real | 5.0 | water temperature (degrees Celsius) |
windspeed | Input | real | 4.0 | wind speed (m/s) |
waterdens | Input | real | 998.02 | density of water (kg/m3) for simplicity, set to static value applicable to water at 21 degC |
latentheat | Input | real | 2257000.0 | latent heat of vaporization (J/kg) for simplicity, set to static value applicable at 1 atm (100 kPa) air pressure |
volheat | Output | real | - | volumetric heat of vaporization (J/m3) |
evapratio | Output | real | - | evaporation ratio: ratio of the heat used to evaporate water to the total heat discharged through the tower |
evapvol | Output | real | - | evaporated volume of water (m3) |
energypervol | Output | real | - | input waste (heat) energy cooled per evaporated volume (J/m3) |
volperenergy | Output | real | - | volume evaporated by units of heat energy (m3/MJ) |
waterusetower | Output | real | - | total volume of water used in cooling tower (m3) |
wateruserecirc | Output | real | - | total volume of water used in recirculating system (m3) |
wateruseonethru | Output | real | - | total volume of water used in once-through system (m3) |
structure_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
aintmass | Output | real | - | intercoil structure mass (kg) |
clgsmass | Output | real | - | gravity support structure for TF coil, PF coil and intercoil support systems (kg) |
coldmass | Output | real | - | total mass of components at cryogenic temperatures (kg) |
fncmass | Output | real | - | PF coil outer support fence mass (kg) |
gsmass | Output | real | - | reactor core gravity support mass (kg) |
fwbs_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
bktlife | Output | real | - | Full power blanket lifetime (years) |
bktlife_cal | Output | real | - | Calendar year blanket lifetime (years) |
coolmass | Output | real | - | mass of water coolant (in shield, blanket, first wall, divertor) [kg] |
vvmass | Output | real | - | vacuum vessel mass [kg] |
denstl | Input | real | 7800.0 | density of steel [kg m^-3] |
denw | Input | real | 19250.0 | density of tungsten [kg m^-3] |
denwc | Input | real | 15630.0 | density of tungsten carbide [kg m^-3] |
dewmkg | Output | real | - | total mass of vacuum vessel + cryostat [kg] (calculated if blktmodel>0) |
emult | Input | real | 1.269 | energy multiplication in blanket and shield |
emultmw | Output | real | - | power due to energy multiplication in blanket and shield [MW] |
fblss | Input | real | 0.09705 | KIT blanket model: steel fraction of breeding zone |
fdiv | Input | real | 0.115 | Solid angle fraction taken by one divertor |
fhcd | Output | real | - | area fraction covered by heating/current drive apparatus plus diagnostics |
fhole | Output | real | - | area fraction taken up by other holes (IFE) |
fwbsshape | Input | integer | 2 | switch for first wall, blanket, shield and vacuum vessel shape:
|
fwlife | Output | real | - | first wall full-power year lifetime (y) |
fwmass | Output | real | - | first wall mass [kg] |
fw_armour_mass | Output | real | - | first wall armour mass [kg] |
fw_armour_thickness | Input | real | 0.005 | first wall armour thickness [m] |
fw_armour_vol | Output | real | - | first wall armour volume [m^3] |
iblanket | Input | integer | 1 | switch for blanket model:
|
iblnkith | Input | integer | 1 | switch for inboard blanket:
|
inuclear | Output | integer | - | switch for nuclear heating in the coils:
|
qnuc | Output | real | - | nuclear heating in the coils (W) ( |
li6enrich | Input | real | 30.0 | lithium-6 enrichment of breeding material (%) |
pnucblkt | Output | real | - | nuclear heating in the blanket [MW] |
pnuc_cp | Output | real | - | Total nuclear heating in the ST centrepost [MW] |
pnuc_cp_sh | Output | real | - | Neutronic shield nuclear heating in the ST centrepost [MW] |
pnuc_cp_tf | Output | real | - | TF neutronic nuclear heating in the ST centrepost [MW] |
pnucdiv | Output | real | - | nuclear heating in the divertor [MW] |
pnucfw | Output | real | - | nuclear heating in the first wall [MW] |
pnuchcd | Output | real | - | nuclear heating in the HCD apparatus and diagnostics [MW] |
pnucloss | Output | real | - | nuclear heating lost via holes [MW] |
pnucvvplus | Output | real | - | nuclear heating to vacuum vessel and beyond [MW] |
pnucshld | Output | real | - | nuclear heating in the shield [MW] |
whtblkt | Output | real | - | mass of blanket [kg] |
whtblss | Output | real | - | mass of blanket - steel part [kg] |
armour_fw_bl_mass | Output | real | - | Total mass of armour, first wall and blanket [kg] |
breeder_f | Input | real | 0.5 | Volume ratio: Li4SiO4/(Be12Ti+Li4SiO4) ( |
breeder_multiplier | Input | real | 0.75 | combined breeder/multipler fraction of blanket by volume |
vfcblkt | Input | real | 0.05295 | He coolant fraction of blanket by volume ( |
vfpblkt | Input | real | 0.1 | He purge gas fraction of blanket by volume ( |
whtblli4sio4 | Output | real | - | mass of lithium orthosilicate in blanket [kg] ( |
whtbltibe12 | Output | real | - | mass of titanium beryllide in blanket [kg] ( |
neut_flux_cp | Output | real | - | Centrepost TF fast neutron flux (E > 0.1 MeV) [m^(-2).^(-1)] This variable is only calculated for superconducting (i_tf_sup = 1 ) spherical tokamal magnet designs (itart = 0) |
f_neut_shield | Input | real | -1.0 | Fraction of nuclear power shielded before the CP magnet (ST) ( neut_absorb = -1 --> a fit on simplified MCNP neutronic calculation is used assuming water cooled (13%) tungesten carbyde ) |
vffwi | Output | real | - | Inboard/outboard FW coolant void fraction |
vffwo | Output | real | - | Inboard/outboard FW coolant void fraction |
psurffwi | Output | real | - | Surface heat flux on first wall [MW] (sum = pradfw) |
psurffwo | Output | real | - | Surface heat flux on first wall [MW] (sum = pradfw) |
volfw | Output | real | - | First wall volume [m3] |
fblss_ccfe | Output | real | - | Fractions of blanket by volume: steel, lithium orthosilicate, titanium beryllide |
fblli2sio4 | Output | real | - | Fractions of blanket by volume: steel, lithium orthosilicate, titanium beryllide |
fbltibe12 | Output | real | - | Fractions of blanket by volume: steel, lithium orthosilicate, titanium beryllide |
breedmat | Input | integer | 1 | breeder material switch (iblanket=2 (KIT HCPB)):
|
densbreed | Output | real | - | density of breeder material [kg m^-3] ( |
fblbe | Input | real | 0.6 | beryllium fraction of blanket by volume (if |
fblbreed | Input | real | 0.154 | breeder fraction of blanket breeding zone by volume ( |
fblhebmi | Input | real | 0.4 | helium fraction of inboard blanket box manifold by volume ( |
fblhebmo | Input | real | 0.4 | helium fraction of outboard blanket box manifold by volume ( |
fblhebpi | Input | real | 0.6595 | helium fraction of inboard blanket back plate by volume ( |
fblhebpo | Input | real | 0.6713 | helium fraction of outboard blanket back plate by volume ( |
hcdportsize | Input | integer | 1 | switch for size of heating/current drive ports (
|
nflutf | Output | real | - | peak fast neutron fluence on TF coil superconductor [n m^-2] ( |
npdiv | Input | integer | 2 | number of divertor ports ( |
nphcdin | Input | integer | 2 | number of inboard ports for heating/current drive ( |
nphcdout | Input | integer | 2 | number of outboard ports for heating/current drive ( |
tbr | Output | real | - | tritium breeding ratio ( |
tritprate | Output | real | - | tritium production rate [g day^-1] ( |
wallpf | Input | real | 1.21 | neutron wall load peaking factor ( |
whtblbreed | Output | real | - | mass of blanket - breeder part [kg] ( |
whtblbe | Output | real | - | mass of blanket - beryllium part [kg] |
iblanket_thickness | Input | integer | 2 | Blanket thickness switch (Do not set blnkith, blnkoth, fwith or fwoth when
|
primary_pumping | Input | integer | 2 | Switch for pumping power for primary coolant (mechanical power only and peak first wall
temperature is only calculated if
|
i_shield_mat | Output | integer | - | Switch for shield material - currently only applied in costing routines
|
secondary_cycle | Output | integer | - | Switch for power conversion cycle:
|
secondary_cycle_liq | Input | integer | 4 | Switch for power conversion cycle for the liquid breeder component of the blanket:
|
coolwh | Input | integer | 1 | Switch for blanket coolant (set via blkttype):
|
afwi | Input | real | 0.008 | inner radius of inboard first wall/blanket coolant channels (stellarator only) [m] |
afwo | Input | real | 0.008 | inner radius of outboard first wall/blanket coolant channels (stellarator only) [m] |
fwcoolant | Input | character | b'helium' | switch for first wall coolant (can be different from blanket coolant):
|
fw_wall | Input | real | 0.003 | wall thickness of first wall coolant channels [m] |
afw | Input | real | 0.006 | radius of first wall cooling channels [m] |
pitch | Input | real | 0.02 | pitch of first wall cooling channels [m] |
fwinlet | Input | real | 573.0 | inlet temperature of first wall coolant [K] |
fwoutlet | Input | real | 823.0 | outlet temperature of first wall coolant [K] |
fwpressure | Input | real | 15500000.0 | first wall coolant pressure [Pa] ( |
tpeak | Input | real | 873.0 | peak first wall temperature [K] |
roughness | Input | real | 1e-06 | first wall channel roughness epsilon [m] |
fw_channel_length | Input | real | 4.0 | Length of a single first wall channel (all in parallel) [m]
( |
peaking_factor | Input | real | 1.0 | peaking factor for first wall heat loads. (Applied separately to inboard and outboard loads. Applies to both neutron and surface loads. Only used to calculate peak temperature - not the coolant flow rate.) |
blpressure | Input | real | 15500000.0 | blanket coolant pressure [Pa] ( |
inlet_temp | Input | real | 573.0 | inlet temperature of blanket coolant [K] ( |
outlet_temp | Input | real | 823.0 | Outlet temperature of blanket coolant [K] (
|
coolp | Input | real | 15500000.0 | blanket coolant pressure [Pa] (stellarator only) |
nblktmodpo | Input | integer | 8 | number of outboard blanket modules in poloidal direction ( |
nblktmodpi | Input | integer | 7 | number of inboard blanket modules in poloidal direction ( |
nblktmodto | Input | integer | 48 | number of outboard blanket modules in toroidal direction ( |
nblktmodti | Input | integer | 32 | number of inboard blanket modules in toroidal direction ( |
tfwmatmax | Input | real | 823.0 | maximum temperature of first wall material [K] ( |
fw_th_conductivity | Input | real | 28.34 | thermal conductivity of first wall material at 293 K (W/m/K) (Temperature dependence is as for unirradiated Eurofer) |
fvoldw | Input | real | 1.74 | area coverage factor for vacuum vessel volume |
fvolsi | Input | real | 1.0 | area coverage factor for inboard shield volume |
fvolso | Input | real | 0.64 | area coverage factor for outboard shield volume |
fwclfr | Input | real | 0.15 | first wall coolant fraction (calculated if |
praddiv | Output | real | - | Radiation power incident on the divertor (MW) |
pradfw | Output | real | - | Radiation power incident on the first wall (MW) |
pradhcd | Output | real | - | Radiation power incident on the heating and current drive system (MW) |
pradloss | Output | real | - | Radiation power lost through holes (eventually hits shield) (MW) Only used for stellarator |
ptfnuc | Output | real | - | nuclear heating in the TF coil (MW) |
ptfnucpm3 | Output | real | - | nuclear heating in the TF coil (MW/m3) ( |
rdewex | Output | real | - | cryostat radius [m] |
zdewex | Output | real | - | cryostat height [m] |
rpf2dewar | Input | real | 0.5 | radial distance between outer edge of largest ( |
vdewex | Output | real | - | cryostat volume [m^3] |
vdewin | Output | real | - | vacuum vessel volume [m^3] |
vfshld | Input | real | 0.25 | coolant void fraction in shield |
volblkt | Output | real | - | volume of blanket [m^3] |
volblkti | Output | real | - | volume of inboard blanket [m^3] |
volblkto | Output | real | - | volume of outboard blanket [m^3] |
volshld | Output | real | - | volume of shield [m^3] |
whtshld | Output | real | - | mass of shield [kg] |
wpenshld | Output | real | - | mass of the penetration shield [kg] |
wtshldi | Output | real | - | mass of inboard shield [kg] |
wtshldo | Output | real | - | mass of outboard shield [kg] |
irefprop | Input | integer | 1 | Switch to use REFPROP routines (stellarator only) |
fblli | Output | real | - | lithium fraction of blanket by volume (stellarator only) |
fblli2o | Input | real | 0.08 | lithium oxide fraction of blanket by volume (stellarator only) |
fbllipb | Input | real | 0.68 | lithium lead fraction of blanket by volume (stellarator only) |
fblvd | Output | real | - | vanadium fraction of blanket by volume (stellarator only) |
wtblli2o | Output | real | - | mass of blanket - Li_2O part [kg] |
wtbllipb | Output | real | - | mass of blanket - Li-Pb part [kg] |
whtblvd | Output | real | - | mass of blanket - vanadium part [kg] |
whtblli | Output | real | - | mass of blanket - lithium part [kg] |
vfblkt | Input | real | 0.25 | coolant void fraction in blanket. |
blktmodel | Output | integer | - | switch for blanket/tritium breeding model (see iblanket):
|
declblkt | Input | real | 0.075 | neutron power deposition decay length of blanket structural material [m] (stellarators only) |
declfw | Input | real | 0.075 | neutron power deposition decay length of first wall structural material [m] (stellarators only) |
declshld | Input | real | 0.075 | neutron power deposition decay length of shield structural material [m] (stellarators only) |
blkttype | Input | integer | 3 | Switch for blanket type:
|
etaiso | Input | real | 0.85 | isentropic efficiency of FW and blanket coolant pumps |
etahtp | Input | real | 0.95 | electrical efficiency of primary coolant pumps BLANKET REFACTOR For DCLL, but to be used by all mods that share blanket library after testing. Thermodynamic Model for primary_pumping == 2 |
ipump | Output | integer | - | Switch for whether the FW and BB are on the same pump system i.e. do they have the same primary coolant or not - =0 FW and BB have the same primary coolant, flow = FWin->FWout->BBin->BBout - =1 FW and BB have the different primary coolant and are on different pump systems |
i_bb_liq | Output | integer | - | Switch for Liquid Metal Breeder Material - =0 PbLi - =1 Li |
icooldual | Output | integer | - | Switch to specify whether breeding blanket is single-cooled or dual-coolant. - =0 Single coolant used for FW and Blanket (H2O or He). Solid Breeder. - =1 Single coolant used for FW and Blanket (H2O or He). Liquid metal breeder circulted for tritium extraction. - =2 Dual coolant: primary coolant (H2O or He) for FW and blanket structure; secondary coolant is self-cooled liquid metal breeder. |
ifci | Output | integer | - | Switch for Flow Channel Insert (FCI) type if liquid metal breeder blanket. - =0 Thin conducting walls, default electrical conductivity (bz_channel_conduct_liq) is Eurofer - =1 Insulating Material, assumed perfect electrical insulator, default density (den_ceramic) is for SiC - =2 Insulating Material, electrical conductivity (bz_channel_conduct_liq) is input (default Eurofer), default density (den_ceramic) is for SiC |
ims | Output | integer | - | Switch for Multi Module Segment (MMS) or Single Modle Segment (SMS) - =0 MMS - =1 SMS |
n_liq_recirc | Input | integer | 10 | Number of liquid metal breeder recirculations per day, for use with icooldual=1 |
r_f_liq_ib | Input | real | 0.5 | Radial fraction of BZ liquid channels |
r_f_liq_ob | Input | real | 0.5 | Radial fraction of BZ liquid channels |
w_f_liq_ib | Input | real | 0.5 | Toroidal fraction of BZ liquid channels |
w_f_liq_ob | Input | real | 0.5 | Toroidal fraction of BZ liquid channels |
den_ceramic | Input | real | 3210.0 | FCI material density |
th_wall_secondary | Input | real | 0.0125 | Liquid metal coolant/breeder wall thickness thin conductor or FCI [m] |
bz_channel_conduct_liq | Input | real | 833000.0 | Liquid metal coolant/breeder thin conductor or FCI wall conductance [A V^-1 m^-1] |
a_bz_liq | Input | real | 0.2 | Toroidal width of the rectangular cooling channel [m] for long poloidal sections of blanket breeding zone |
b_bz_liq | Input | real | 0.2 | Radial width of the rectangular cooling channel [m] for long poloidal sections of blanket breeding zone |
nopol | Input | integer | 2 | Number of poloidal sections in a liquid metal breeder/coolant channel for module/segment |
nopipes | Input | integer | 4 | Number of Liquid metal breeder/coolant channels per module/segment |
den_liq | Input | real | 9500.0 | Liquid metal breeder/coolant density [kg m^-3] |
wht_liq | Output | real | - | Liquid metal |
wht_liq_ib | Output | real | - | Liquid metal |
wht_liq_ob | Output | real | - | Liquid metal |
specific_heat_liq | Input | real | 190.0 | Liquid metal breeder/coolant specific heat [J kg^-1 K^-1] |
thermal_conductivity_liq | Input | real | 30.0 | Liquid metal breeder/coolant thermal conductivity [W m^-1 K^-1] |
dynamic_viscosity_liq | Output | real | - | Liquid metal breeder/coolant dynamic viscosity [Pa s] |
electrical_conductivity_liq | Output | real | - | Liquid metal breeder/coolant electrical conductivity [Ohm m] |
hartmann_liq | Output | real | - | Hartmann number |
b_mag_blkt | Input | real | [5. 5.] | Toroidal Magnetic feild strength for IB/OB blanket [T] |
etaiso_liq | Input | real | 0.85 | Isentropic efficiency of blanket liquid breeder/coolant pumps |
blpressure_liq | Input | real | 1700000.0 | blanket liquid metal breeder/coolant pressure [Pa] |
inlet_temp_liq | Input | real | 570.0 | Inlet (scan var 68) and Outlet (scan var 69) temperature of the liquid breeder/coolant [K] |
outlet_temp_liq | Input | real | 720.0 | Inlet (scan var 68) and Outlet (scan var 69) temperature of the liquid breeder/coolant [K] |
rhof_fw | Output | real | - | Density of the FW primary coolant |
visc_fw | Output | real | - | Viscosity of the FW primary coolant |
rhof_bl | Output | real | - | Density of the blanket primary coolant |
visc_bl | Output | real | - | Viscosity of the blanket primary coolant |
cp_fw | Output | real | - | Spesific heat for FW and blanket primary coolant(s) |
cv_fw | Output | real | - | Spesific heat for FW and blanket primary coolant(s) |
cp_bl | Output | real | - | Spesific heat for FW and blanket primary coolant(s) |
cv_bl | Output | real | - | Spesific heat for FW and blanket primary coolant(s) |
f_nuc_pow_bz_struct | Input | real | 0.34 | For a dual-coolant blanket, fraction of BZ power cooled by primary coolant |
f_nuc_pow_bz_liq | Input | real | 0.66 | For a dual-coolant blanket, fraction of BZ self-cooled power (secondary coolant) |
pnuc_fw_ratio_dcll | Input | real | 0.14 | For a dual-coolant blanket, ratio of FW/Blanket nuclear power as fraction of total |
pnuc_blkt_ratio_dcll | Input | real | 0.86 | For a dual-coolant blanket, ratio of FW/Blanket nuclear power as fraction of total |
bzfllengi_n_rad | Input | integer | 4 | Number of radial and poloidal sections that make up the total primary coolant flow length in a blanket module (IB and OB) |
bzfllengi_n_pol | Input | integer | 2 | Number of radial and poloidal sections that make up the total primary coolant flow length in a blanket module (IB and OB) |
bzfllengo_n_rad | Input | integer | 4 | Number of radial and poloidal sections that make up the total primary coolant flow length in a blanket module (IB and OB) |
bzfllengo_n_pol | Input | integer | 2 | Number of radial and poloidal sections that make up the total primary coolant flow length in a blanket module (IB and OB) |
bzfllengi_n_rad_liq | Input | integer | 2 | Number of radial and poloidal sections that make up the total secondary coolant/breeder flow length in a blanket module (IB and OB) |
bzfllengi_n_pol_liq | Input | integer | 2 | Number of radial and poloidal sections that make up the total secondary coolant/breeder flow length in a blanket module (IB and OB) |
bzfllengo_n_rad_liq | Input | integer | 2 | Number of radial and poloidal sections that make up the total secondary coolant/breeder flow length in a blanket module (IB and OB) |
bzfllengo_n_pol_liq | Input | integer | 2 | Number of radial and poloidal sections that make up the total secondary coolant/breeder flow length in a blanket module (IB and OB) |
constraint_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
auxmin | Input | real | 0.1 | minimum auxiliary power (MW) ( |
betpmx | Input | real | 0.19 | maximum poloidal beta ( |
bigqmin | Input | real | 10.0 | minimum fusion gain Q ( |
bmxlim | Input | real | 12.0 | maximum peak toroidal field (T) ( |
fauxmn | Input | real | 1.0 | f-value for minimum auxiliary power ( |
fbeta | Input | real | 1.0 | f-value for epsilon beta-poloidal ( |
fbetap | Input | real | 1.0 | f-value for poloidal beta ( |
fbetatry | Input | real | 1.0 | f-value for beta limit ( |
fbetatry_lower | Input | real | 1.0 | f-value for (lower) beta limit ( |
fcpttf | Input | real | 1.0 | f-value for TF coil current per turn upper limit
( |
fcwr | Input | real | 1.0 | f-value for conducting wall radius / rminor limit
( |
fdene | Input | real | 1.0 | f-value for density limit ( |
fdivcol | Input | real | 1.0 | f-value for divertor collisionality ( |
fdtmp | Input | real | 1.0 | f-value for first wall coolant temperature rise
( |
fecrh_ignition | Input | real | 1.0 | f-value for ecrh ignition constraint
( |
fflutf | Input | real | 1.0 | f-value for neutron fluence on TF coil ( |
ffuspow | Input | real | 1.0 | f-value for maximum fusion power ( |
fgamcd | Input | real | 1.0 | f-value for current drive gamma ( |
fhldiv | Input | real | 1.0 | f-value for divertor heat load ( |
fiooic | Input | real | 0.5 | f-value for TF coil operating current / critical current ratio
( |
fipir | Input | real | 1.0 | f-value for Ip/Irod upper limit constraint equation icc = 46 iteration variable ixc = 72 |
fjohc | Input | real | 1.0 | f-value for central solenoid current at end-of-flattop
( |
fjohc0 | Input | real | 1.0 | f-value for central solenoid current at beginning of pulse
( |
fjprot | Input | real | 1.0 | f-value for TF coil winding pack current density
( |
flhthresh | Input | real | 1.0 | f-value for L-H power threshold ( |
fmva | Input | real | 1.0 | f-value for maximum MVA ( |
fnbshinef | Input | real | 1.0 | f-value for maximum neutral beam shine-through fraction
( |
fncycle | Input | real | 1.0 | f-value for minimum CS coil stress load cycles
( |
fnesep | Input | real | 1.0 | f-value for Eich critical separatrix density
( |
foh_stress | Input | real | 1.0 | f-value for Tresca yield criterion in Central Solenoid
( |
fpeakb | Input | real | 1.0 | f-value for maximum toroidal field ( |
fpinj | Input | real | 1.0 | f-value for injection power ( |
fpnetel | Input | real | 1.0 | f-value for net electric power ( |
fportsz | Input | real | 1.0 | f-value for neutral beam tangency radius limit
( |
fpsepbqar | Input | real | 1.0 | f-value for maximum Psep*Bt/qAR limit ( |
fpsepr | Input | real | 1.0 | f-value for maximum Psep/R limit ( |
fptemp | Input | real | 1.0 | f-value for peak centrepost temperature ( |
fptfnuc | Input | real | 1.0 | f-value for maximum TF coil nuclear heating ( |
fq | Input | real | 1.0 | f-value for edge safety factor ( |
fqval | Input | real | 1.0 | f-value for Q ( |
fradpwr | Input | real | 0.99 | f-value for core radiation power limit ( |
fradwall | Input | real | 1.0 | f-value for upper limit on radiation wall load ( |
freinke | Input | real | 1.0 | f-value for Reinke detachment criterion ( |
frminor | Input | real | 1.0 | f-value for minor radius limit ( |
fstrcase | Input | real | 1.0 | f-value for maximum TF coil case Tresca yield criterion
( |
fstrcond | Input | real | 1.0 | f-value for maxiumum TF coil conduit Tresca yield criterion
( |
fstr_wp | Input | real | 1.0 | f-value for maxiumum TF coil strain absolute value
( |
fmaxvvstress | Input | real | 1.0 | f-value for maximum permitted stress of the VV
( |
ftbr | Input | real | 1.0 | f-value for minimum tritium breeding ratio ( |
ftburn | Input | real | 1.0 | f-value for minimum burn time ( |
ftcycl | Input | real | 1.0 | f-value for cycle time ( |
ftmargoh | Input | real | 1.0 | f-value for central solenoid temperature margin
( |
ftmargtf | Input | real | 1.0 | f-value for TF coil temperature margin ( |
ftohs | Input | real | 1.0 | f-value for plasma current ramp-up time ( |
ftpeak | Input | real | 1.0 | f-value for first wall peak temperature ( |
fvdump | Input | real | 1.0 | f-value for dump voltage ( |
fvs | Input | real | 1.0 | f-value for flux-swing (V-s) requirement (STEADY STATE)
( |
fvvhe | Input | real | 1.0 | f-value for vacuum vessel He concentration limit ( |
fwalld | Input | real | 1.0 | f-value for maximum wall load ( |
fzeffmax | Input | real | 1.0 | f-value for maximum zeff ( |
gammax | Input | real | 2.0 | maximum current drive gamma ( |
maxradwallload | Input | real | 1.0 | Maximum permitted radiation wall load (MW/m^2) ( |
mvalim | Input | real | 40.0 | maximum MVA limit ( |
nbshinefmax | Input | real | 0.001 | maximum neutral beam shine-through fraction ( |
nflutfmax | Input | real | 1e+23 | max fast neutron fluence on TF coil (n/m2) ( |
pdivtlim | Input | real | 150.0 | Minimum pdivt [MW] ( |
peakfactrad | Input | real | 3.33 | peaking factor for radiation wall load ( |
peakradwallload | Output | real | - | Peak radiation wall load (MW/m^2) ( |
pnetelin | Input | real | 1000.0 | required net electric power (MW) ( |
powfmax | Input | real | 1500.0 | maximum fusion power (MW) ( |
psepbqarmax | Input | real | 9.5 | maximum ratio of Psep*Bt/qAR (MWT/m) ( |
pseprmax | Input | real | 25.0 | maximum ratio of power crossing the separatrix to plasma major radius (Psep/R) (MW/m)
( |
ptfnucmax | Input | real | 0.001 | maximum nuclear heating in TF coil (MW/m3) ( |
tbrmin | Input | real | 1.1 | minimum tritium breeding ratio ( |
tbrnmn | Input | real | 1.0 | minimum burn time (s) (KE - no longer itv., see issue #706) |
tcycmn | Output | real | - | minimum cycle time (s) ( |
tohsmn | Input | real | 1.0 | minimum plasma current ramp-up time (s) ( |
vvhealw | Input | real | 1.0 | allowed maximum helium concentration in vacuum vessel at end of plant life (appm)
( |
walalw | Input | real | 1.0 | allowable neutron wall-load (MW/m2) ( |
taulimit | Input | real | 5.0 | Lower limit on taup/taueff the ratio of alpha particle to energy confinement
times ( |
ftaulimit | Input | real | 1.0 | f-value for lower limit on taup/taueff the ratio of alpha particle to energy
confinement times ( |
fniterpump | Input | real | 1.0 | f-value for constraint that number of pumps < tfno
( |
zeffmax | Input | real | 3.6 | maximum value for Zeff ( |
fpoloidalpower | Input | real | 1.0 | f-value for constraint on rate of change of energy in poloidal field
( |
fpsep | Input | real | 1.0 | f-value to ensure separatrix power is less than value from Kallenbach divertor (Not required as constraint 69 is an equality) |
fcqt | Input | real | 1.0 | TF coil quench temparature remains below tmax_croco
( |
cs_fatigue_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
residual_sig_hoop | Input | real | 240000000.0 | residual hoop stress in strucutal material (Pa) |
n_cycle | Output | real | - | Allowable number of cycles for CS stress model |
n_cycle_min | Input | real | 20000.0 | Minimum llowable number of cycles for CS stress model |
t_crack_radial | Input | real | 0.006 | Initial depth of crack in thickness of conduit (m) |
t_crack_vertical | Input | real | 0.00089 | Inital vertical crack size (m) |
t_structural_radial | Input | real | 0.07 | Thickness of CS conductor conduit (m) |
t_structural_vertical | Input | real | 0.022 | Vertical thickness of CS conductor conduit (m) |
bkt_life_csf | Output | real | - | Switch to pass bkt_life cycles to n_cycle_min |
sf_vertical_crack | Input | real | 2.0 | Safety factor for vertical crack size (-) |
sf_radial_crack | Input | real | 2.0 | Safety factor for radial crack size (-) |
sf_fast_fracture | Input | real | 1.5 | safety factor for stress intensity factor (-) |
paris_coefficient | Input | real | 6.5e-13 | Paris equation material coefficient (-) |
paris_power_law | Input | real | 3.5 | Paris equation material power law (-) |
walker_coefficient | Input | real | 0.436 | walker coefficent (-) |
fracture_toughness | Input | real | 200.0 | fracture toughness (MPa m^1/2) |
pf_power_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
acptmax | Output | real | - | average of currents in PF circuits (kA) |
ensxpfm | Output | real | - | maximum stored energy in the PF circuits (MJ) |
iscenr | Input | integer | 2 | Switch for PF coil energy storage option:
|
pfckts | Output | real | - | number of PF coil circuits |
spfbusl | Output | real | - | total PF coil circuit bus length (m) |
spsmva | Output | real | - | sum of PF power supply ratings (MVA) |
srcktpm | Output | real | - | sum of resistive PF coil power (kW) |
vpfskv | Output | real | - | PF coil voltage (kV) |
peakpoloidalpower | Output | real | - | Peak absolute rate of change of stored energy in poloidal field (MW) |
maxpoloidalpower | Input | real | 1000.0 | Maximum permitted absolute rate of change of stored energy in poloidal field (MW) |
poloidalpower | Output | real | - | Poloidal power usage at time t (MW) |
define_iteration_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
dummy | Output | real | - |
tfcoil_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
acasetf | Output | real | - | external case area per coil (inboard leg) (m2) |
acasetfo | Output | real | - | external case area per coil (outboard leg) (m2) |
acndttf | Output | real | - | area of the cable conduit (m2) |
acond | Output | real | - | Winding pack conductor area [m2] Does not include the area of voids and central helium channel |
acstf | Output | real | - | Cable space area (per turn) [m2] Includes the area of voids and central helium channel |
insulation_area | Output | real | - | single turn insulation area (m2) |
aiwp | Output | real | - | winding pack turn insulation area per coil (m2) |
sig_tf_case_max | Input | real | 600000000.0 | Allowable maximum shear stress (Tresca criterion) in TF coil case (Pa) |
sig_tf_wp_max | Input | real | 600000000.0 | Allowable maximum shear stress (Tresca criterion) in TF coil conduit (Pa) Allowable Tresca stress in TF coil structural material (Pa) |
arealeg | Output | real | - | outboard TF leg area (m2) |
aswp | Output | real | - | winding pack structure area (m2) |
avwp | Output | real | - | winding pack void (He coolant) area (m2) |
awphec | Output | real | - | winding pack He coil area (m2) |
bcritsc | Input | real | 24.0 | upper critical field (T) for Nb3Sn superconductor at zero temperature and
strain ( |
bmaxtf | Output | real | - | mean peak field at TF coil (T) |
bmaxtfrp | Output | real | - | peak field at TF conductor with ripple (T) |
casestr | Output | real | - | case strain |
casthi | Output | real | - | inboard TF coil case plasma side thickness (m) (calculated for stellarators) |
casthi_fraction | Input | real | 0.05 | inboard TF coil case plasma side thickness as a fraction of tfcth |
casthi_is_fraction | Output | logical | - | logical switch to make casthi a fraction of TF coil thickness ( |
casths | Output | real | - | inboard TF coil sidewall case thickness (m) (calculated for stellarators) |
casths_fraction | Input | real | 0.06 | inboard TF coil sidewall case thickness as a fraction of tftort |
tfc_sidewall_is_fraction | Output | logical | - | logical switch to make casths a fraction of TF coil thickness ( |
t_conductor | Output | real | - | Conductor (cable + steel conduit) area averaged dimension [m] |
t_turn_tf | Output | real | - | TF coil turn edge length including turn insulation [m] If the turn is not a square (i_tf_turns_integer = 1) a squared turn of equivelent size is use to calculated this quantity If the t_turn_tf is non zero, cpttf is calculated |
t_turn_tf_is_input | Output | logical | - | Boolean switch to activated when the user set the TF coil turn dimensions Not an input |
f_t_turn_tf | Input | real | 1.0 | f-value for TF turn edge length constraint If the turn is not a square (i_tf_turns_integer = 1) a squared turn of equivelent size is use for this constraint iteration variable ixc = 175 constraint equation icc = 86 |
t_turn_tf_max | Input | real | 0.05000000074505806 | TF turn edge length including turn insulation upper limit [m] If the turn is not a square (i_tf_turns_integer = 1) a squared turn of equivelent size is use for this constraint constraint equation icc = 86 |
t_cable_tf | Output | real | - | TF coil superconducting cable squared/rounded dimensions [m] If the turn is not a square (i_tf_turns_integer = 1) a squared cable of equivelent size is use to calculated this quantity If the t_cable_tf is non zero, cpttf is calculated |
t_cable_tf_is_input | Output | logical | - | Boolean switch to activated when the user set the TF coil cable dimensions Not an input |
acs | Output | real | - | Area of space inside conductor (m2) |
cdtfleg | Output | real | - | TF outboard leg current density (A/m2) (resistive coils only) |
cforce | Output | real | - | centering force on inboard leg (per coil) (N/m) |
cplen | Output | real | - | length of TF coil inboard leg ('centrepost') ( |
cpttf | Input | real | 70000.0 | TF coil current per turn (A). (calculated for stellarators) (calculated for
integer-turn TF coils |
cpttf_max | Input | real | 90000.0 | Max TF coil current per turn [A]. (for stellarators and |
dcase | Input | real | 8000.0 | density of coil case (kg/m3) |
dcond | Input | real | [6080. 6080. 6070. 6080. 6080. 8500. 6070. 8500. 8500.] | density of superconductor type given by i_tf_sc_mat/isumatoh/isumatpf (kg/m3) |
dcondins | Input | real | 1800.0 | density of conduit + ground-wall insulation (kg/m3) |
dhecoil | Input | real | 0.005 | diameter of central helium channel in TF winding (m) |
estotftgj | Output | real | - | total stored energy in the toroidal field (GJ) |
b_crit_upper_nbti | Input | real | 14.86 | upper critical field of GL_nbti |
t_crit_nbti | Input | real | 9.04 | critical temperature of GL_nbti |
max_force_density | Output | real | - | Maximal (WP averaged) force density in TF coils at 1 point. (MN/m3) |
fcutfsu | Input | real | 0.69 | copper fraction of cable conductor (TF coils) (iteration variable 59) |
fhts | Input | real | 0.5 | technology adjustment factor for critical current density fit for isumat..=2 Bi-2212 superconductor, to describe the level of technology assumed (i.e. to account for stress, fatigue, radiation, AC losses, joints or manufacturing variations; 1.0 would be very optimistic) |
insstrain | Output | real | - | Radial strain in insulator |
i_tf_stress_model | Input | integer | 1 | Switch for the TF coil stress model 0 : Generalized plane strain formulation, Issues #977 and #991, O(n^3) 1 : Old plane stress model (only for SC) 2 : Axisymmetric extended plane strain, Issues #1414 and #998, O(n) |
i_tf_tresca | Output | integer | - | Switch for TF coil conduit Tresca stress criterion: 0 : Tresca (no adjustment); 1 : Tresca with CEA adjustment factors (radial+2%, vertical+60%) |
i_tf_wp_geom | Input | integer | -1 | Switch for TF WP geometry selection 0 : Rectangular geometry 1 : Double rectangular geometry 2 : Trapezoidal geometry (constant lateral casing thickness) Default setting for backward compatibility if i_tf_turns_integer = 0 : Double rectangular if i_tf_turns_integer = 1 : Rectangular |
i_tf_case_geom | Output | integer | - | Switch for TF case geometry selection 0 : Circular front case (ITER design) 1 : Straight front case |
i_tf_turns_integer | Output | integer | - | Switch for TF coil integer/non-integer turns: 0 : non-integer turns 1 : integer turns |
i_tf_sc_mat | Input | integer | 1 | Switch for superconductor material in TF coils:
|
i_tf_sup | Input | integer | 1 | Switch for TF coil conductor model:
|
i_tf_shape | Output | integer | - | Switch for TF coil toroidal shape:
|
i_tf_cond_eyoung_axial | Output | integer | - | Switch for the behavior of the TF coil conductor elastic axial properties
|
i_tf_cond_eyoung_trans | Input | integer | 1 | Switch for the behavior of the elastic properties of the TF coil
conductorin the transverse direction. Only active if
|
n_pancake | Input | integer | 10 | Number of pancakes in TF coil. Only used if |
n_layer | Input | integer | 20 | Number of layers in TF coil. Only used if |
n_rad_per_layer | Input | integer | 100 | Size of the arrays per layers storing the radial dependent stress quantities (stresses, strain displacement etc..) |
i_tf_bucking | Input | integer | -1 | Switch for TF inboard suport structure design: Default setting for backward compatibility - if copper resistive TF (i_tf_sup = 0) : Free standing TF without bucking structure - if Superconducting TF (i_tf_sup = 1) : Free standing TF with a steel casing - if aluminium TF (i_tf_sup = 2) : Free standing TF with a bucking structure Rem : the case is a bucking structure - =0 : Free standing TF without case/bucking cyliner (only a conductor layer) - =1 : Free standing TF with a case/bucking cylinder made of - if copper resistive TF (i_tf_sup = 0) : used defined bucking cylinder - if Superconducting TF (i_tf_sup = 1) : Steel casing - if aluminium resisitive TF (i_tf_sup = 2) : used defined bucking cylinder - =2 : The TF is in contact with the CS : "bucked and wedged design" Fast version : thin TF-CS interface neglected in the stress calculations (3 layers) The CS is frictionally decoupled from the TF, does not carry axial tension - =3 : The TF is in contact with the CS : "bucked and wedged design" Full version : thin TF-CS Kapton interface introduced in the stress calculations (4 layers) The CS and kaptop are frictionally decoupled from the TF, do not carry axial tension |
n_tf_graded_layers | Input | integer | 1 | Number of layers of different stress properties in the WP. If |
n_tf_stress_layers | Output | integer | - | Number of layers considered for the inboard TF stress calculations set in initial.f90 from i_tf_bucking and n_tf_graded_layers |
n_tf_wp_layers | Input | integer | 5 | Maximum number of layers that can be considered in the TF coil composited/smeared stress analysis. This is the layers of one turn, not the entire WP. Default: 5. void, conductor, copper, conduit, insulation. |
jbus | Input | real | 1250000.0 | bussing current density (A/m2) |
j_crit_str_tf | Output | real | - | j_crit_str : superconductor strand critical current density under operating conditions (A/m2). Necessary for the cost calculation in $/kAm |
j_crit_str_0 | Input | real | [5.96905476e+08 1.92550153e+09 7.24544683e+08 5.49858624e+08 6.69284510e+08 0.00000000e+00 8.98964415e+08 1.15875300e+09 8.65652123e+08] | j_crit_str_pf_0 : superconductor strand critical current density at 6 T and 4.2 K (A/m2) Necessary for the cost calculation in $/kAm |
jwdgcrt | Output | real | - | critical current density for winding pack (A/m2) |
jwdgpro | Output | real | - | allowable TF coil winding pack current density, for dump temperature rise protection (A/m2) |
jwptf | Output | real | - | winding pack engineering current density (A/m2) |
oacdcp | Output | real | - | Overall current density in TF coil inboard legs midplane (A/m2) Rem SK : Not used in tfcoil to set the current any more. Should not be used as iteration variable 12 any more. It is now calculated. |
eyoung_ins | Input | real | 100000000.0 | Insulator Young's modulus [Pa]. Default value (1.0D8) setup the following values - SC TF, eyoung_ins = 20 Gpa (default value from DDD11-2 v2 2 (2009)) - Al TF, eyoung_ins = 2.5 GPa (Kapton polymer) |
eyoung_steel | Input | real | 205000000000.0 | Steel case Young's modulus (Pa) (default value from DDD11-2 v2 2 (2009)) |
eyoung_cond_axial | Input | real | 660000000.0 | SC TF coil conductor Young's modulus in the parallel (along the wire/tape)
direction [Pa]
Set by user input only if |
eyoung_cond_trans | Output | real | - | SC TF coil conductor Young's modulus in the transverse direction [Pa]
Set by user input only if |
eyoung_res_tf_buck | Input | real | 150000000000.0 | Resistive TF magnets bucking cylinder young modulus (Pa) |
eyoung_copper | Input | real | 117000000000.0 | Copper young modulus. Default value taken from wikipedia |
eyoung_al | Input | real | 69000000000.0 | Aluminium young modulus. Default value taken from wikipedia |
poisson_steel | Input | real | 0.3 | Steel Poisson's ratio, Source : https://www.engineeringtoolbox.com/metals-poissons-ratio-d_1268.html |
poisson_copper | Input | real | 0.35 | Copper Poisson's ratio. Source : https://www.engineeringtoolbox.com/poissons-ratio-d_1224.html |
poisson_al | Input | real | 0.35 | Aluminium Poisson's ratio. Source : https://www.engineeringtoolbox.com/poissons-ratio-d_1224.html |
poisson_ins | Input | real | 0.34 | Insulation Poisson's ratio. Default: Kapton. Source : DuPont™ Kapton® HN datasheet. |
poisson_cond_axial | Input | real | 0.30000001192092896 | SC TF coil conductor Poisson's ratio in the parallel-transverse direction |
poisson_cond_trans | Input | real | 0.30000001192092896 | SC TF coil conductor Poisson's ratio in the transverse-transverse direction |
rbmax | Output | real | - | Radius of maximum TF B-field (m) |
tflegres | Output | real | - | TF coil leg resistance (ohm) |
toroidalgap | Input | real | 1.0 | Minimal distance between two toroidal coils. (m) |
ftoroidalgap | Input | real | 1.0 | F-value for minimum tftort ( |
ripmax | Input | real | 1.0 | aximum allowable toroidal field ripple amplitude at plasma edge (%) |
ripple | Output | real | - | peak/average toroidal field ripple at plasma edge (%) |
ritfc | Output | real | - | total (summed) current in TF coils (A) |
n_radial_array | Parameter | integer | 50 | Size of the radial distribution arrays per layers used for stress, strain and displacement distibution |
radial_array | Output | real | - | Array refining the radii of the stress calculations arrays |
sig_tf_r | Output | real | - | TF Inboard leg radial stress in steel r distribution at mid-plane [Pa] |
sig_tf_t | Output | real | - | TF Inboard leg tangential stress in steel r distribution at mid-plane [Pa] |
deflect | Output | real | - | TF coil radial deflection (displacement) radial distribution [m] |
sig_tf_z | Output | real | - | TF Inboard leg vertical tensile stress in steel at mid-plane [Pa] |
sig_tf_vmises | Output | real | - | TF Inboard leg Von-Mises stress in steel r distribution at mid-plane [Pa] |
sig_tf_tresca | Output | real | - | TF Inboard leg maximum shear stress (Tresca criterion) in steel r distribution at mid-plane [Pa] |
sig_tf_cs_bucked | Output | real | - | Maximum shear stress (Tresca criterion) in CS structures at CS flux swing [Pa]:
Quantity only computed for bucked and wedged design ( |
sig_tf_case | Output | real | - | Maximum shear stress (Tresca criterion) in TF casing steel structures (Pa) |
sig_tf_wp | Output | real | - | |
str_cs_con_res | Input | real | -0.005 | Residual manufacturing strain in CS superconductor material |
str_pf_con_res | Input | real | -0.005 | Residual manufacturing strain in PF superconductor material |
str_tf_con_res | Input | real | -0.005 | Residual manufacturing strain in TF superconductor material
If |
str_wp | Output | real | - | Axial (vertical) strain in the TF coil winding pack found by
self-consistent stress/strain calculation.
if |
str_wp_max | Input | real | 0.007 | Maximum allowed absolute value of the strain in the TF coil
( |
i_str_wp | Input | integer | 1 | Switch for the behavior of the TF strain used to compute the strain-dependent critical surface:
|
quench_model | Input | character | b'exponential ' | switch for TF coil quench model (Only applies to REBCO magnet at present, issue #522):
|
time1 | Output | real | - | Time at which TF quench is detected (s) |
tcritsc | Input | real | 16.0 | critical temperature (K) for superconductor at zero field and strain ( |
tdmptf | Input | real | 10.0 | fast discharge time for TF coil in event of quench (s) ( For REBCO model, meaning depends on quench_model:
|
tfareain | Output | real | - | Area of inboard midplane TF legs (m2) |
tfbusl | Output | real | - | TF coil bus length (m) |
tfbusmas | Output | real | - | TF coil bus mass (kg) |
tfckw | Output | real | - | available DC power for charging the TF coils (kW) |
tfcmw | Output | real | - | Peak power per TF power supply (MW) |
tfcpmw | Output | real | - | Peak resistive TF coil inboard leg power (MW) |
tfjtsmw | Output | real | - | TF joints resistive power losses (MW) |
tfcryoarea | Output | real | - | surface area of toroidal shells covering TF coils (m2) |
tficrn | Output | real | - | TF coil half-width - inner bore (m) |
tfind | Output | real | - | TF coil inductance (H) |
tfinsgap | Input | real | 0.01 | TF coil WP insertion gap (m) |
tflegmw | Output | real | - | TF coil outboard leg resistive power (MW) |
rhocp | Output | real | - | TF coil inboard leg resistivity [Ohm-m]. If |
rhotfleg | Output | real | - | Resistivity of a TF coil leg (Ohm-m) |
rhotfbus | Input | real | 2.62e-08 | Resistivity of a TF coil bus (Ohm-m). Default value takes the same res as the leg one |
frhocp | Input | real | 1.0 | Centrepost resistivity enhancement factor. For |
frholeg | Input | real | 1.0 | Ouboard legs resistivity enhancement factor. Only used for |
i_cp_joints | Input | integer | -1 | Switch for CP demoutable joints type -= 0 : Clampled joints -= 1 : Sliding joints Default value (-1) choses : Sliding joints for resistive magnets (i_tf_sup = 0, 2) Clampled joints for superconducting magents (i_tf_sup = 1) |
rho_tf_joints | Input | real | 2.5e-10 | TF joints surfacic resistivity [ohm.m]. Feldmetal joints assumed. |
n_tf_joints_contact | Input | integer | 6 | Number of contact per turn |
n_tf_joints | Input | integer | 4 | Number of joints Ex: n_tf_joints = 2 for top and bottom CP joints |
th_joint_contact | Input | real | 0.03 | TF sliding joints contact pad width [m] |
pres_joints | Output | real | - | Calculated TF joints resistive power losses [W] |
tfleng | Output | real | - | TF coil circumference (m) |
eff_tf_cryo | Input | real | -1.0 | TF cryoplant efficiency (compared to pefect Carnot cycle). Using -1 set the default value depending on magnet technology:
|
n_tf | Input | real | 16.0 | Number of TF coils (default = 50 for stellarators). Number of TF coils outer legs for ST |
tfocrn | Output | real | - | TF coil half-width - outer bore (m) |
tfsai | Output | real | - | area of the inboard TF coil legs (m2) |
tfsao | Output | real | - | area of the outboard TF coil legs (m2) |
tftmp | Input | real | 4.5 | peak helium coolant temperature in TF coils and PF coils (K) |
tftort | Input | real | 1.0 | TF coil toroidal thickness (m) |
thicndut | Input | real | 0.0008 | conduit insulation thickness (m) |
layer_ins | Output | real | - | Additional insulation thickness between layers (m) |
thkcas | Input | real | 0.3 | inboard TF coil case outer (non-plasma side) thickness (m) ( |
dr_tf_wp | Output | real | - | radial thickness of winding pack (m) ( |
thwcndut | Input | real | 0.008 | TF coil conduit case thickness (m) ( |
tinstf | Input | real | 0.018 | Thickness of the ground insulation layer surrounding (m)
- Superconductor TF ( |
tmargmin_tf | Output | real | - | minimum allowable temperature margin : TF coils (K) |
tmargmin_cs | Output | real | - | minimum allowable temperature margin : CS (K) |
tmargmin | Output | real | - | minimum allowable temperature margin : TFC AND CS (K) |
temp_margin | Output | real | - | temperature margin (K) |
tmargtf | Output | real | - | TF coil temperature margin (K) |
tmaxpro | Input | real | 150.0 | maximum temp rise during a quench for protection (K) |
tmax_croco | Input | real | 200.0 | CroCo strand: maximum permitted temp during a quench (K) |
croco_quench_temperature | Output | real | - | CroCo strand: Actual temp reached during a quench (K) |
tmpcry | Input | real | 4.5 | coil temperature for cryogenic plant power calculation (K) |
n_tf_turn | Output | real | - | number of turns per TF coil |
vdalw | Input | real | 20.0 | max voltage across TF coil during quench (kV) ( |
vforce | Output | real | - | vertical tension on inboard leg/coil (N) |
f_vforce_inboard | Input | real | 0.5 | Fraction of the total vertical force taken by the TF inboard leg tension
Not used for resistive |
vforce_outboard | Output | real | - | Vertical tension on outboard leg/coil (N) |
vftf | Input | real | 0.4 | coolant fraction of TFC 'cable' ( |
voltfleg | Output | real | - | volume of each TF coil outboard leg (m3) |
vtfkv | Output | real | - | TF coil voltage for resistive coil including bus (kV) |
vtfskv | Output | real | - | voltage across a TF coil during quench (kV) |
whtcas | Output | real | - | mass per coil of external case (kg) |
whtcon | Output | real | - | TF coil conductor mass per coil (kg/coil).
For |
whtconcu | Output | real | - | copper mass in TF coil conductor (kg/coil).
For |
whtconal | Output | real | - | Aluminium mass in TF coil conductor (kg/coil).
For |
whtconin | Output | real | - | conduit insulation mass in TF coil conductor (kg/coil) |
whtconsc | Output | real | - | superconductor mass in TF coil cable (kg/coil) |
whtconsh | Output | real | - | steel conduit mass in TF coil conductor (kg/coil) |
whtgw | Output | real | - | mass of ground-wall insulation layer per coil (kg/coil) |
whttf | Output | real | - | total mass of the TF coils (kg) |
wwp1 | Output | real | - | width of first step of winding pack (m) |
wwp2 | Output | real | - | width of second step of winding pack (m) |
dthet | Output | real | - | angle of arc i (rad) |
radctf | Output | real | - | radius of arc i (m) |
xarc | Output | real | - | x location of arc point i on surface (m) |
xctfc | Output | real | - | x location of arc centre i (m) |
yarc | Output | real | - | y location of arc point i on surface (m) |
yctfc | Output | real | - | y location of arc centre i (m) |
tfa | Output | real | - | Horizontal radius of inside edge of TF coil (m) |
tfb | Output | real | - | Vertical radius of inside edge of TF coil (m) |
drtop | Output | real | - | centrepost taper maximum radius adjustment (m) |
dztop | Output | real | - | centrepost taper height adjustment (m) |
etapump | Input | real | 0.8 | centrepost coolant pump efficiency |
fcoolcp | Input | real | 0.3 | coolant fraction of TF coil inboard legs ( |
fcoolleg | Input | real | 0.2 | coolant fraction of TF coil outboard legs |
a_cp_cool | Output | real | - | Centrepost cooling area toroidal cross-section (constant over the whole CP) |
ncool | Output | real | - | number of centrepost coolant tubes |
ppump | Output | real | - | centrepost coolant pump power (W) |
prescp | Output | real | - | resistive power in the centrepost (itart=1) [W].
If |
presleg | Output | real | - | Summed resistive power in the TF coil legs [W]. Remain 0 if |
ptempalw | Input | real | 473.15 | maximum peak centrepost temperature (K) ( |
rcool | Input | real | 0.005 | average radius of coolant channel (m) ( |
tcoolin | Input | real | 313.15 | centrepost coolant inlet temperature (K) |
dtiocool | Output | real | - | inlet / outlet TF coil coolant temperature rise (K) |
tcpav | Input | real | 373.15 | Average temperature of centrepost called CP (K). Only used for resistive coils
to compute the resisitive heating. Must be an iteration variable for
ST ( |
tcpav2 | Output | real | - | Computed centrepost average temperature (K) (for consistency) |
tlegav | Input | real | -1.0 | Average temperature of the TF outboard legs [K]. If |
tcpmax | Output | real | - | peak centrepost temperature (K) |
vcool | Input | real | 20.0 | inlet centrepost coolant flow speed at midplane (m/s) ( |
vol_cond_cp | Output | real | - | Exact conductor volume in the centrepost (m3) |
whtcp | Output | real | - | mass of TF coil inboard legs (kg) |
whttflgs | Output | real | - | mass of the TF coil legs (kg) |
cryo_cool_req | Output | real | - | Cryo cooling requirement at helium temp 4.5K (kW) |
theta1_coil | Input | real | 45.0 | The angle of the outboard arc forming the TF coil current center line [deg] |
theta1_vv | Input | real | 1.0 | The angle of the outboard arc forming the Vacuum Vessel current center line [deg] |
max_vv_stress | Input | real | 143000000.0 | The allowable peak maximum shear stress in the vacuum vessel due to quench and fast discharge of the TF coils [Pa] |
ife_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
maxmat | Parameter | integer | 8 | Total number of materials in IFE device. Material numbers are as follows:
|
bldr | Input | real | 1.0 | radial thickness of IFE blanket (m; calculated |
bldrc | Input | real | 1.0 | radial thickness of IFE curtain (m; |
bldzl | Input | real | 4.0 | vertical thickness of IFE blanket below chamber (m) |
bldzu | Input | real | 4.0 | vertical thickness of IFE blanket above chamber (m) |
blmatf | Input | real | [[0.05 0. 0.45 0. 0.2 0. 0.3 0. 0. ] [0.05 0. 0.45 0. 0.2 0. 0.3 0. 0. ] [0.05 0. 0.45 0. 0.2 0. 0.3 0. 0. ]] | IFE blanket material fractions |
blmatm | Output | real | - | IFE blanket material masses (kg) |
blmatv | Output | real | - | IFE blanket material volumes (m3) |
blvol | Output | real | - | IFE blanket volume (m3) |
cdriv0 | Input | real | 154.3 | IFE generic/laser driver cost at edrive=0 (M$) |
cdriv1 | Input | real | 163.2 | IFE low energy heavy ion beam driver cost extrapolated to |
cdriv2 | Input | real | 244.9 | IFE high energy heavy ion beam driver cost extrapolated to |
cdriv3 | Input | real | 1.463 | IFE driver cost ($/J wall plug) ( |
chdzl | Input | real | 9.0 | vertical thickness of IFE chamber below centre (m) |
chdzu | Input | real | 9.0 | vertical thickness of IFE chamber above centre (m) |
chmatf | Input | real | [1. 0. 0. 0. 0. 0. 0. 0. 0.] | IFE chamber material fractions |
chmatm | Output | real | - | IFE chamber material masses (kg) |
chmatv | Output | real | - | IFE chamber material volumes (m3) |
chrad | Input | real | 6.5 | radius of IFE chamber (m) ( |
chvol | Output | real | - | IFE chamber volume (m3) |
dcdrv0 | Input | real | 111.4 | IFE generic/laser driver cost gradient (M$/MJ) |
dcdrv1 | Input | real | 78.0 | HIB driver cost gradient at low energy (M$/MJ) |
dcdrv2 | Input | real | 59.9 | HIB driver cost gradient at high energy (M$/MJ) |
drveff | Input | real | 0.28 | IFE driver wall plug to target efficiency ( |
edrive | Input | real | 5000000.0 | IFE driver energy (J) ( |
etadrv | Output | real | - | IFE driver wall plug to target efficiency |
etali | Input | real | 0.4 | IFE lithium pump wall plug efficiency ( |
etave | Input | real | [0.082 0.079 0.076 0.073 0.069 0.066 0.062 0.059 0.055 0.051] | IFE driver efficiency vs driver energy ( |
fauxbop | Input | real | 0.06 | fraction of gross electric power to balance-of-plant (IFE) |
fbreed | Input | real | 0.51 | fraction of breeder external to device core |
fburn | Input | real | 0.3333 | IFE burn fraction (fraction of tritium fused/target) |
flirad | Input | real | 0.78 | radius of FLiBe/lithium inlet (m) ( |
frrmax | Input | real | 1.0 | f-value for maximum IFE repetition rate ( |
fwdr | Input | real | 0.01 | radial thickness of IFE first wall (m) |
fwdzl | Input | real | 0.01 | vertical thickness of IFE first wall below chamber (m) |
fwdzu | Input | real | 0.01 | vertical thickness of IFE first wall above chamber (m) |
fwmatf | Input | real | [[0.05 0. 0.95 0. 0. 0. 0. 0. 0. ] [0.05 0. 0.95 0. 0. 0. 0. 0. 0. ] [0.05 0. 0.95 0. 0. 0. 0. 0. 0. ]] | IFE first wall material fractions |
fwmatm | Output | real | - | IFE first wall material masses (kg) |
fwmatv | Output | real | - | IFE first wall material volumes (kg) |
fwvol | Output | real | - | IFE first wall volume (m3) |
gain | Output | real | - | IFE target gain |
gainve | Input | real | [ 60. 95. 115. 125. 133. 141. 152. 160. 165. 170.] | IFE target gain vs driver energy ( |
htpmw_ife | Output | real | - | IFE heat transport system electrical pump power (MW) |
ife | Output | integer | - | Switch for IFE option:
|
ifedrv | Input | integer | 2 | Switch for type of IFE driver:
|
ifetyp | Output | integer | - | Switch for type of IFE device build:
|
lipmw | Output | real | - | IFE lithium pump power (MW; |
mcdriv | Input | real | 1.0 | IFE driver cost multiplier |
mflibe | Output | real | - | total mass of FLiBe (kg) |
pdrive | Input | real | 23000000.0 | IFE driver power reaching target (W) ( |
pfusife | Input | real | 1000.0 | IFE input fusion power (MW) ( |
pifecr | Input | real | 10.0 | IFE cryogenic power requirements (MW) |
ptargf | Input | real | 2.0 | IFE target factory power at 6 Hz repetition rate (MW) |
r1 | Output | real | - | IFE device radial build (m) |
r2 | Output | real | - | IFE device radial build (m) |
r3 | Output | real | - | IFE device radial build (m) |
r4 | Output | real | - | IFE device radial build (m) |
r5 | Output | real | - | IFE device radial build (m) |
r6 | Output | real | - | IFE device radial build (m) |
r7 | Output | real | - | IFE device radial build (m) |
reprat | Output | real | - | IFE driver repetition rate (Hz) |
rrin | Input | real | 6.0 | Input IFE repetition rate (Hz) ( |
rrmax | Input | real | 20.0 | maximum IFE repetition rate (Hz) |
shdr | Input | real | 1.7 | radial thickness of IFE shield (m) |
shdzl | Input | real | 5.0 | vertical thickness of IFE shield below chamber (m) |
shdzu | Input | real | 5.0 | vertical thickness of IFE shield above chamber (m) |
shmatf | Input | real | [[0.05 0.19 0. 0. 0. 0.665 0.095 0. 0. ] [0.05 0.19 0. 0. 0. 0.665 0.095 0. 0. ] [0.05 0.19 0. 0. 0. 0.665 0.095 0. 0. ]] | IFE shield material fractions |
shmatm | Output | real | - | IFE shield material masses (kg) |
shmatv | Output | real | - | IFE shield material volumes (kg) |
shvol | Output | real | - | IFE shield volume (m3) |
sombdr | Input | real | 2.7 | radius of cylindrical blanket section below chamber ( |
somtdr | Input | real | 2.7 | radius of cylindrical blanket section above chamber ( |
taufall | Output | real | - | Lithium Fall Time (s) |
tdspmw | Input | real | 0.01 | IFE target delivery system power (MW) |
tfacmw | Output | real | - | IFE target factory power (MW) |
tgain | Input | real | 85.0 | IFE target gain (if |
uccarb | Input | real | 50.0 | cost of carbon cloth ($/kg) |
ucconc | Input | real | 0.1 | cost of concrete ($/kg) |
ucflib | Input | real | 84.0 | cost of FLiBe ($/kg) |
uctarg | Input | real | 0.3 | cost of IFE target ($/target) |
v1dr | Output | real | - | radial thickness of IFE void between first wall and blanket (m) |
v1dzl | Output | real | - | vertical thickness of IFE void 1 below chamber (m) |
v1dzu | Output | real | - | vertical thickness of IFE void 1 above chamber (m) |
v1matf | Input | real | [[1. 0. 0. 0. 0. 0. 0. 0. 0.] [1. 0. 0. 0. 0. 0. 0. 0. 0.] [1. 0. 0. 0. 0. 0. 0. 0. 0.]] | IFE void 1 material fractions |
v1matm | Output | real | - | IFE void 1 material masses (kg) |
v1matv | Output | real | - | IFE void 1 material volumes (kg) |
v1vol | Output | real | - | IFE void 1 volume (m3) |
v2dr | Input | real | 2.0 | radial thickness of IFE void between blanket and shield (m) |
v2dzl | Input | real | 7.0 | vertical thickness of IFE void 2 below chamber (m) |
v2dzu | Input | real | 7.0 | vertical thickness of IFE void 2 above chamber (m) |
v2matf | Input | real | [[1. 0. 0. 0. 0. 0. 0. 0. 0.] [1. 0. 0. 0. 0. 0. 0. 0. 0.] [1. 0. 0. 0. 0. 0. 0. 0. 0.]] | IFE void 2 material fractions |
v2matm | Output | real | - | IFE void 2 material masses (kg) |
v2matv | Output | real | - | IFE void 2 material volumes (kg) |
v2vol | Output | real | - | IFE void 2 volume (m3) |
v3dr | Input | real | 43.3 | radial thickness of IFE void outside shield (m) |
v3dzl | Input | real | 30.0 | vertical thickness of IFE void 3 below chamber (m) |
v3dzu | Input | real | 20.0 | vertical thickness of IFE void 3 above chamber (m) |
v3matf | Input | real | [[1. 0. 0. 0. 0. 0. 0. 0. 0.] [1. 0. 0. 0. 0. 0. 0. 0. 0.] [1. 0. 0. 0. 0. 0. 0. 0. 0.]] | IFE void 3 material fractions |
v3matm | Output | real | - | IFE void 3 material masses (kg) |
v3matv | Output | real | - | IFE void 3 material volumes (kg) |
v3vol | Output | real | - | IFE void 3 volume (m3) |
zl1 | Output | real | - | IFE vertical build below centre (m) |
zl2 | Output | real | - | IFE vertical build below centre (m) |
zl3 | Output | real | - | IFE vertical build below centre (m) |
zl4 | Output | real | - | IFE vertical build below centre (m) |
zl5 | Output | real | - | IFE vertical build below centre (m) |
zl6 | Output | real | - | IFE vertical build below centre (m) |
zl7 | Output | real | - | IFE vertical build below centre (m) |
zu1 | Output | real | - | IFE vertical build above centre (m) |
zu2 | Output | real | - | IFE vertical build above centre (m) |
zu3 | Output | real | - | IFE vertical build above centre (m) |
zu4 | Output | real | - | IFE vertical build above centre (m) |
zu5 | Output | real | - | IFE vertical build above centre (m) |
zu6 | Output | real | - | IFE vertical build above centre (m) |
zu7 | Output | real | - | IFE vertical build above centre (m) |
stellarator_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
istell | Output | integer | - | Switch for stellarator option (set via
|
bmn | Input | real | 0.001 | relative radial field perturbation |
f_asym | Input | real | 1.0 | divertor heat load peaking factor |
f_rad | Input | real | 0.85 | radiated power fraction in SOL |
f_w | Input | real | 0.5 | island size fraction factor |
fdivwet | Input | real | 0.333333333333333 | wetted fraction of the divertor area |
flpitch | Input | real | 0.001 | field line pitch (rad) |
hportamax | Output | real | - | maximum available area for horizontal ports (m2) |
hportpmax | Output | real | - | maximum available poloidal extent for horizontal ports (m) |
hporttmax | Output | real | - | maximum available toroidal extent for horizontal ports (m) |
iotabar | Input | real | 1.0 | rotational transform (reciprocal of tokamak q) for stellarator confinement time scaling laws |
isthtr | Input | integer | 3 | Switch for stellarator auxiliary heating method:
|
m_res | Input | integer | 5 | poloidal resonance number (1) |
max_gyrotron_frequency | Output | real | - | Maximal available gyrotron frequency (input parameter) (Hz) |
n_res | Input | integer | 5 | toroidal resonance number (1) |
shear | Input | real | 0.5 | magnetic shear, derivative of iotabar (1) |
te0_ecrh_achievable | Output | real | - | maximal central electron temperature as achievable by the ECRH, input. (keV) |
vportamax | Output | real | - | maximum available area for vertical ports (m2) |
vportpmax | Output | real | - | maximum available poloidal extent for vertical ports (m) |
vporttmax | Output | real | - | maximum available toroidal extent for vertical ports (m) |
powerht_constraint | Output | real | - | |
powerscaling_constraint | Output | real | - |
constants
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
iotty | Parameter | integer | 6 | Standard output unit identifier |
nout | Parameter | integer | 11 | Output file unit identifier |
nplot | Parameter | integer | 12 | Plot data file unit identifier |
mfile | Parameter | integer | 13 | Machine-optimised output file unit |
vfile | Parameter | integer | 14 | Verbose diagnostics file |
opt_file | Parameter | integer | 15 | Optimisation information output file number |
sig_file | Parameter | integer | 16 | TF inboard stress radial distributions file number |
degrad | Parameter | real | 0.01745329251D0 | degrees to radians, = pi/180 |
echarge | Parameter | real | 1.60217733D-19 | electron charge [C] |
emass | Parameter | real | 9.10938370D-31 | electron mass [kg] |
mproton | Parameter | real | 1.6726231D-27 | proton mass [kg] |
pi | Parameter | real | 3.1415926535897932D0 | pi |
rmu0 | Parameter | real | 1.256637062D-6 | permeability of free space [H/m] |
twopi | Parameter | real | 6.2831853071795862D0 | 2 pi |
umass | Parameter | real | 1.660538921D-27 | unified atomic mass unit [kg |
epsilon0 | Parameter | real | 8.85418781D-12 | permittivity of free space [Farad/m] |
cph2o | Parameter | real | 4180.0D0 | specific heat capacity of water (J/kg/K) |
dcopper | Input | real | 8900.0 | density of copper (kg/m3) |
dalu | Input | real | 2700.0 | density of aluminium (kg/m3) |
denh2o | Parameter | real | 985.0D0 | density of water (kg/m3) |
k_copper | Parameter | real | 330.0D0 | Copper thermal conductivity (W/m/K) |
kh2o | Parameter | real | 0.651D0 | thermal conductivity of water (W/m/K) |
muh2o | Parameter | real | 4.71D-4 | water dynamic viscosity (kg/m/s) |
n_day_year | Parameter | real | 365.2425D0 | Average number of days in a year |
stellarator_configuration
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
stella_config_name | Output | character | - | |
stella_config_symmetry | Output | integer | - | |
stella_config_coilspermodule | Output | integer | - | |
stella_config_rmajor_ref | Output | real | - | |
stella_config_rminor_ref | Output | real | - | |
stella_config_coil_rmajor | Output | real | - | |
stella_config_coil_rminor | Output | real | - | |
stella_config_aspect_ref | Output | real | - | |
stella_config_bt_ref | Output | real | - | |
stella_config_wp_area | Output | real | - | |
stella_config_wp_bmax | Output | real | - | |
stella_config_i0 | Output | real | - | |
stella_config_a1 | Output | real | - | |
stella_config_a2 | Output | real | - | |
stella_config_dmin | Output | real | - | |
stella_config_inductance | Output | real | - | |
stella_config_coilsurface | Output | real | - | |
stella_config_coillength | Output | real | - | |
stella_config_max_portsize_width | Output | real | - | |
stella_config_maximal_coil_height | Output | real | - | |
stella_config_min_plasma_coil_distance | Output | real | - | |
stella_config_derivative_min_lcfs_coils_dist | Output | real | - | |
stella_config_plasma_volume | Output | real | - | |
stella_config_plasma_surface | Output | real | - | |
stella_config_wp_ratio | Output | real | - | |
stella_config_max_force_density | Output | real | - | |
stella_config_max_force_density_mnm | Output | real | - | |
stella_config_min_bend_radius | Output | real | - | |
stella_config_epseff | Output | real | - | |
stella_config_max_lateral_force_density | Output | real | - | |
stella_config_max_radial_force_density | Output | real | - | |
stella_config_centering_force_max_mn | Output | real | - | |
stella_config_centering_force_min_mn | Output | real | - | |
stella_config_centering_force_avg_mn | Output | real | - | |
stella_config_neutron_peakfactor | Output | real | - | |
sc_d11_star_mono_input | Output | real | - | |
sc_nu_star_mono_input | Output | real | - |
pfcoil_module
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
nef | Output | integer | - | |
nfxf | Output | integer | - | |
ricpf | Output | real | - | |
ssq0 | Output | real | - | |
sig_axial | Output | real | - | |
sig_hoop | Output | real | - | |
axial_force | Output | real | - | |
rfxf | Output | real | - | |
zfxf | Output | real | - | |
cfxf | Output | real | - | |
xind | Output | real | - | |
rcls | Output | real | - | |
zcls | Output | real | - | |
ccls | Output | real | - | |
ccl0 | Output | real | - | |
bpf2 | Output | real | - | |
vsdum | Output | real | - | |
first_call | Input | logical | 1 | |
cslimit | Output | logical | - | |
conductorpf* | Variable | type | - | |
croco_strand* | Variable | type | - |
scan_module
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
ipnscns | Parameter | integer | 1000 | Maximum number of scan points |
ipnscnv | Parameter | integer | 81 | Number of available scan variables |
noutvars | Parameter | integer | 84 | |
width | Parameter | integer | 110 | |
scan_dim | Input | integer | 1 | 1-D or 2-D scan switch (1=1D, 2=2D) |
isweep | Output | integer | - | Number of scan points to calculate |
isweep_2 | Output | integer | - | Number of 2D scan points to calculate |
nsweep | Input | integer | 1 | Switch denoting quantity to scan:
inlet_temp_liq' : Inlet temperature of blanket liquid metal coolant/breeder (K)
<LI> 69 outlet_temp_liq' : Outlet temperature of blanket liquid metal coolant/breeder (K)
blpressure_liq' : Blanket liquid metal breeder/coolant pressure (Pa)
<LI> 71 n_liq_recirc' : Selected number of liquid metal breeder recirculations per day
bz_channel_conduct_liq' : Conductance of liquid metal breeder duct walls (A V-1 m-1)
<LI> 73 pnuc_fw_ratio_dcll' : Ratio of FW nuclear power as fraction of total (FW+BB)
|
nsweep_2 | Input | integer | 3 | nsweep_2 /3/ : switch denoting quantity to scan for 2D scan: |
sweep | Output | real | - | sweep(ipnscns) /../: actual values to use in scan |
sweep_2 | Output | real | - | sweep_2(ipnscns) /../: actual values to use in 2D scan |
first_call_1d | Input | logical | 1 | |
first_call_2d | Input | logical | 1 |
neoclassics_constants
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
no_roots | Parameter | integer | 30 |
neoclassics_module
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
species | Variable | character | (/"e", "D", "T", "a"/) | |
densities | Output | real | - | |
temperatures | Output | real | - | |
dr_densities | Output | real | - | |
dr_temperatures | Output | real | - | |
roots | Variable | real | 0 | |
weights | Variable | real | 0 | |
nu | Variable | real | 0 | |
nu_star | Variable | real | 0 | |
nu_star_averaged | Variable | real | 0 | |
vd | Variable | real | 0 | |
kt | Variable | real | 0 | |
er | Variable | real | 0.0 | |
iota | Variable | real | 1.0d0 | |
d11_mono | Variable | real | 0 | |
d11_plateau | Variable | real | 0 | |
nu_star_mono_input | Output | real | - | |
d11_star_mono_input | Output | real | - | |
d13_star_mono_input | Output | real | - | |
d111 | Variable | real | 0 | |
d112 | Variable | real | 0 | |
d113 | Variable | real | 0 | |
q_flux | Variable | real | 0 | |
gamma_flux | Variable | real | 0 | |
d31_mono | Variable | real | 0 | |
eps_eff | Variable | real | 1d-5 | |
r_eff | Variable | real | 0 |
primary_pumping_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
gamma_he | Input | real | 1.667 | ratio of specific heats for helium ( |
t_in_bb | Input | real | 573.13 | temperature in FW and blanket coolant at blanket entrance ( |
t_out_bb | Input | real | 773.13 | temperature in FW and blanket coolant at blanket exit ( |
p_he | Input | real | 8000000.0 | pressure in FW and blanket coolant at pump exit ( |
dp_he | Input | real | 550000.0 | pressure drop in FW and blanket coolant including heat exchanger and pipes ( |
dp_fw_blkt | Input | real | 150000.0 | pressure drop in FW and blanket coolant including heat exchanger and pipes ( |
dp_fw | Input | real | 150000.0 | pressure drop in FW coolant including heat exchanger and pipes ( |
dp_blkt | Input | real | 3500.0 | pressure drop in blanket coolant including heat exchanger and pipes ( |
dp_liq | Input | real | 10000000.0 | pressure drop in liquid metal blanket coolant including heat exchanger and pipes ( |
htpmw_fw_blkt | Output | real | - | mechanical pumping power for FW and blanket including heat exchanger and
pipes ( |
divertor_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
adas | Output | real | - | area divertor / area main plasma (along separatrix) |
anginc | Input | real | 0.262 | angle of incidence of field line on plate (rad) |
beta_div | Input | real | 1.0 | field line angle wrt divertor target plate (degrees) |
betai | Input | real | 1.0 | poloidal plane angle between divertor plate and leg, inboard (rad) |
betao | Input | real | 1.0 | poloidal plane angle between divertor plate and leg, outboard (rad) |
bpsout | Input | real | 0.6 | reference B_p at outboard divertor strike point (T) |
c1div | Input | real | 0.45 | fitting coefficient to adjust ptpdiv, ppdiv |
c2div | Input | real | -7.0 | fitting coefficient to adjust ptpdiv, ppdiv |
c3div | Input | real | 0.54 | fitting coefficient to adjust ptpdiv, ppdiv |
c4div | Input | real | -3.6 | fitting coefficient to adjust ptpdiv, ppdiv |
c5div | Input | real | 0.7 | fitting coefficient to adjust ptpdiv, ppdiv |
c6div | Output | real | - | fitting coefficient to adjust ptpdiv, ppdiv |
delld | Input | real | 1.0 | coeff for power distribution along main plasma |
dendiv | Output | real | - | plasma density at divertor (10**20 /m3) |
densin | Output | real | - | density at plate (on separatrix) (10**20 /m3) |
divclfr | Input | real | 0.3 | divertor coolant fraction |
divdens | Input | real | 10000.0 | divertor structure density (kg/m3) |
divdum | Output | integer | - | switch for divertor Zeff model:
|
divfix | Input | real | 0.2 | divertor structure vertical thickness (m) |
divmas | Output | real | - | divertor plate mass (kg) |
divplt | Input | real | 0.035 | divertor plate thickness (m) (from Spears, Sept 1990) |
divsur | Output | real | - | divertor surface area (m2) |
fdfs | Input | real | 10.0 | radial gradient ratio |
fdiva | Input | real | 1.11 | divertor area fudge factor (for ITER, Sept 1990) |
fhout | Output | real | - | fraction of power to outboard divertor (for single null) |
fififi | Input | real | 0.004 | coefficient for gamdiv |
flux_exp | Input | real | 2.0 | The plasma flux expansion in the divertor (default 2; Wade 2020) |
frrp | Input | real | 0.4 | fraction of radiated power to plate |
hldiv | Output | real | - | divertor heat load (MW/m2) |
i_hldiv | Output | integer | - | switch for user input hldiv:
|
hldivlim | Input | real | 5.0 | heat load limit (MW/m2) |
ksic | Input | real | 0.8 | power fraction for outboard double-null scrape-off plasma |
lamp | Output | real | - | power flow width (m) |
minstang | Output | real | - | minimum strike angle for heat flux calculation |
omegan | Input | real | 1.0 | pressure ratio (nT)_plasma / (nT)_scrape-off |
omlarg | Output | real | - | power spillage to private flux factor |
ppdivr | Output | real | - | peak heat load at plate (with radiation) (MW/m2) |
prn1 | Input | real | 0.285 | n-scrape-off / n-average plasma; (input for |
ptpdiv | Output | real | - | peak temperature at the plate (eV) |
rconl | Output | real | - | connection length ratio, outboard side |
rlclolcn | Output | real | - | ratio of collision length / connection length |
rlenmax | Input | real | 0.5 | maximum value for length ratio (rlclolcn) ( |
rsrd | Output | real | - | effective separatrix/divertor radius ratio |
tconl | Output | real | - | main plasma connection length (m) |
tdiv | Input | real | 2.0 | temperature at divertor (eV) (input for stellarator only, calculated for tokamaks) |
tsep | Output | real | - | temperature at the separatrix (eV) |
xparain | Input | real | 2100.0 | parallel heat transport coefficient (m2/s) |
xpertin | Input | real | 2.0 | perpendicular heat transport coefficient (m2/s) |
zeffdiv | Input | real | 1.0 | Zeff in the divertor region (if |
current_drive_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
beamwd | Input | real | 0.58 | width of neutral beam duct where it passes between the TF coils (m) T Inoue et al, Design of neutral beam system for ITER-FEAT, Fusion Engineering and Design, Volumes 56-57, October 2001, Pages 517-521) |
bigq | Output | real | - | Fusion gain; P_fusion / (P_injection + P_ohmic) |
bootstrap_current_fraction | Output | real | - | bootstrap current fraction (enforced; see i_bootstrap_current) |
bootstrap_current_fraction_max | Input | real | 0.9 | maximum fraction of plasma current from bootstrap; if |
bscf_iter89 | Output | real | - | bootstrap current fraction, ITER 1989 model |
bscf_nevins | Output | real | - | bootstrap current fraction, Nevins et al model |
bscf_sauter | Output | real | - | bootstrap current fraction, Sauter et al model |
bscf_wilson | Output | real | - | bootstrap current fraction, Wilson et al model |
bscf_sakai | Output | real | - | Bootstrap current fraction, Sakai et al model |
cboot | Input | real | 1.0 | bootstrap current fraction multiplier ( |
cnbeam | Output | real | - | neutral beam current (A) |
diacf_hender | Output | real | - | diamagnetic current fraction, Hender fit |
diacf_scene | Output | real | - | diamagnetic current fraction, SCENE fit |
diamagnetic_current_fraction | Output | real | - | diamagnetic current fraction |
echpwr | Output | real | - | ECH power (MW) |
echwpow | Output | real | - | ECH wall plug power (MW) |
effcd | Output | real | - | current drive efficiency (A/W) |
harnum | Input | real | 2.0 | cyclotron harmonic frequency number, used in cut-off function |
wave_mode | Output | integer | - | Switch for ECRH wave mode :
|
enbeam | Input | real | 1000.0 | neutral beam energy (keV) ( |
etacd | Output | real | - | auxiliary power wall plug to injector efficiency |
etacdfix | Output | real | - | secondary auxiliary power wall plug to injector efficiency |
etaech | Input | real | 0.3 | ECH wall plug to injector efficiency |
etalh | Input | real | 0.3 | lower hybrid wall plug to injector efficiency |
etanbi | Input | real | 0.3 | neutral beam wall plug to injector efficiency |
fpion | Input | real | 0.5 | fraction of beam energy to ions |
pnbitot | Output | real | - | neutral beam power entering vacuum vessel |
pscf_scene | Output | real | - | Pfirsch-Schlüter current fraction, SCENE fit |
nbshinemw | Output | real | - | neutral beam shine-through power |
feffcd | Input | real | 1.0 | current drive efficiency fudge factor ( |
forbitloss | Output | real | - | fraction of neutral beam power lost after ionisation but before thermalisation (orbit loss fraction) |
frbeam | Input | real | 1.05 | R_tangential / R_major for neutral beam injection |
ftritbm | Input | real | 1e-06 | fraction of beam that is tritium |
gamcd | Output | real | - | normalised current drive efficiency (1.0e20 A/(W m^2)) |
gamma_ecrh | Input | real | 0.35 | User input ECRH gamma (1.0e20 A/(W m^2)) |
xi_ebw | Input | real | 0.8 | User scaling input for EBW plasma heating. Default 0.43 |
iefrf | Input | integer | 5 | Switch for current drive efficiency model:
|
iefrffix | Output | integer | - | Switch for 2nd current drive efficiency model:
|
irfcd | Input | integer | 1 | Switch for current drive calculation:
|
nbshinef | Output | real | - | neutral beam shine-through fraction |
nbshield | Input | real | 0.5 | neutral beam duct shielding thickness (m) |
pheat | Output | real | - | heating power not used for current drive (MW) ( |
pheatfix | Output | real | - | secondary fixed heating power not used for current drive (MW) |
pinjalw | Input | real | 150.0 | maximum allowable value for injected power (MW) ( |
pinjemw | Output | real | - | auxiliary injected power to electrons (MW) |
pinjimw | Output | real | - | auxiliary injected power to ions (MW) |
pinjmw | Output | real | - | total auxiliary injected power (MW) |
pinjfixmw | Output | real | - | secondary total fixed auxiliary injected power (MW) |
plasma_current_internal_fraction | Output | real | - | plasma current fraction driven internally (Bootstrap + Diamagnetic + PS) |
plhybd | Output | real | - | lower hybrid injection power (MW) |
pnbeam | Output | real | - | neutral beam injection power (MW) |
porbitlossmw | Output | real | - | neutral beam power lost after ionisation but before thermalisation (orbit loss power) (MW) |
ps_current_fraction | Output | real | - | Pfirsch-Schlüter current fraction |
pwplh | Output | real | - | lower hybrid wall plug power (MW) |
pwpnb | Output | real | - | neutral beam wall plug power (MW) |
rtanbeam | Output | real | - | neutral beam centreline tangency radius (m) |
rtanmax | Output | real | - | maximum tangency radius for centreline of beam (m) |
taubeam | Output | real | - | neutral beam e-decay lengths to plasma centre |
tbeamin | Input | real | 3.0 | permitted neutral beam e-decay lengths to plasma centre |
numerics
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
ipnvars | Parameter | integer | 175 | ipnvars FIX : total number of variables available for iteration |
ipeqns | Parameter | integer | 91 | ipeqns FIX : number of constraint equations available |
ipnfoms | Parameter | integer | 19 | ipnfoms FIX : number of available figures of merit |
ipvlam | Parameter | integer | ipeqns+2*ipnvars+1 | |
iptnt | Parameter | integer | (ipeqns*(3*ipeqns+13))/2 | |
ipvp1 | Parameter | integer | ipnvars+1 | |
ioptimz | Input | integer | 1 | ioptimz /1/ : code operation switch:
minmax /7/ : switch for figure-of-merit (see lablmm for descriptions) negative => maximise, positive => minimise |
minmax | Input | integer | 7 | |
lablmm | Input | character | [b'major radius ' b'not used ' b'neutron wall load ' b'P_tf + P_pf ' b'fusion gain ' b'cost of electricity ' b'capital cost ' b'aspect ratio ' b'divertor heat load ' b'toroidal field ' b'total injected power ' b'H plant capital cost ' b'H production rate ' b'pulse length ' b'plant availability ' b'min R0, max tau_burn ' b'net electrical output ' b'Null figure of merit ' b'max Q, max t_burn '] | lablmm(ipnfoms) : labels describing figures of merit:
|
ncalls | Output | integer | - | ncalls : number of function calls during solution |
neqns | Output | integer | - | neqns /0/ : number of equality constraints to be satisfied |
nfev1 | Output | integer | - | nfev1 : number of calls to FCNHYB (HYBRD function caller) made |
nfev2 | Output | integer | - | nfev2 : number of calls to FCNVMC1 (VMCON function caller) made |
nineqns | Output | integer | - | nineqns /0/ : number of inequality constraints VMCON must satisfy (leave at zero for now) |
nvar | Input | integer | 16 | nvar /16/ : number of iteration variables to use |
nviter | Output | integer | - | nviter : number of VMCON iterations performed icc(ipeqns) /0/ : array defining which constraint equations to activate (see lablcc for descriptions) |
icc | Output | integer | - | |
active_constraints | Output | logical | - | active_constraints(ipeqns) : Logical array showing which constraints are active |
lablcc | Input | character | [b'Beta consistency ' b'Global power balance consistency ' b'Ion power balance ' b'Electron power balance ' b'Density upper limit ' b'(Epsilon x beta-pol) upper limit ' b'Beam ion density consistency ' b'Neutron wall load upper limit ' b'Fusion power upper limit ' b'Toroidal field 1/R consistency ' b'Radial build consistency ' b'Volt second lower limit ' b'Burn time lower limit ' b'NBI decay lengths consistency ' b'L-H power threshold limit ' b'Net electric power lower limit ' b'Radiation fraction upper limit ' b'Divertor heat load upper limit ' b'MVA upper limit ' b'Beam tangency radius upper limit ' b'Plasma minor radius lower limit ' b'Divertor collisionality upper lim' b'Conducting shell radius upper lim' b'Beta upper limit ' b'Peak toroidal field upper limit ' b'CS coil EOF current density limit' b'CS coil BOP current density limit' b'Fusion gain Q lower limit ' b'Inboard radial build consistency ' b'Injection power upper limit ' b'TF coil case stress upper limit ' b'TF coil conduit stress upper lim ' b'I_op / I_critical (TF coil) ' b'Dump voltage upper limit ' b'J_winding pack/J_protection limit' b'TF coil temp. margin lower limit ' b'Current drive gamma limit ' b'1st wall coolant temp rise limit ' b'First wall peak temperature limit' b'Start-up inj. power lower limit ' b'Plasma curr. ramp time lower lim ' b'Cycle time lower limit ' b'Average centrepost temperature ' b'Peak centrepost temp. upper limit' b'Edge safety factor lower limit ' b'Ip/Irod upper limit ' b'TF coil tor. thickness upper lim ' b'Poloidal beta upper limit ' b'RFP reversal parameter < 0 ' b'IFE repetition rate upper limit ' b'Startup volt-seconds consistency ' b'Tritium breeding ratio lower lim ' b'Neutron fluence on TF coil limit ' b'Peak TF coil nucl. heating limit ' b'Vessel helium concentration limit' b'Psep / R upper limit ' b'TF coil leg rad width lower limit' b'TF coil leg rad width lower limit' b'NB shine-through frac upper limit' b'CS temperature margin lower limit' b'Minimum availability value ' b'taup/taueff ' b'number of ITER-like vacuum pumps ' b'Zeff limit ' b'Dump time set by VV stress ' b'Rate of change of energy in field' b'Upper Lim. on Radiation Wall load' b'Upper Lim. on Psep * Bt / q A R ' b'pdivt < psep_kallenbach divertor ' b'Separatrix temp consistency ' b'Separatrix density consistency ' b'CS Tresca yield criterion ' b'Psep >= Plh + Paux ' b'TFC quench < tmax_croco ' b'TFC current/copper area < Max ' b'Eich critical separatrix density ' b'TFC current per turn upper limit ' b'Reinke criterion fZ lower limit ' b'Peak CS field upper limit ' b'pdivt lower limit ' b'ne0 > neped ' b'toroidalgap > tftort ' b'available_space > required_space ' b'beta > betalim_lower ' b'CP lifetime ' b'TFC turn dimension ' b'Cryogenic plant power ' b'TF coil strain absolute value ' b'CS current/copper area < Max ' b'CS stress load cycles ' b'ECRH ignitability '] | lablcc(ipeqns) : labels describing constraint equations (corresponding itvs)
|
ixc | Output | integer | - | ixc(ipnvars) /0/ : array defining which iteration variables to activate (see lablxc for descriptions) |
lablxc | Input | character | [b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' '] | lablxc(ipnvars) : labels describing iteration variables
|
name_xc | Variable | character | - | |
sqsumsq | Output | real | - | sqsumsq : sqrt of the sum of the square of the constraint residuals |
objf_name | Input | character | b' ' | Description of the objective function |
norm_objf | Output | real | - | Normalised objective function (figure of merit) |
epsfcn | Input | real | 0.001 | epsfcn /1.0e-3/ : finite difference step length for HYBRD/VMCON derivatives |
epsvmc | Input | real | 1e-06 | epsvmc /1.0e-6/ : error tolerance for VMCON |
factor | Input | real | 0.1 | factor /0.1/ : used in HYBRD for first step size |
ftol | Input | real | 0.0001 | ftol /1.0e-4/ : error tolerance for HYBRD |
boundl | Input | real | [9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99 9.e-99] | boundl(ipnvars) /../ : lower bounds used on ixc variables during VMCON optimisation runs |
boundu | Input | real | [9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99 9.e+99] | |
bondl | Output | real | - | |
bondu | Output | real | - | |
rcm | Output | real | - | |
resdl | Output | real | - | |
scafc | Output | real | - | |
scale | Output | real | - | |
xcm | Output | real | - | |
xcs | Output | real | - | |
vlam | Output | real | - |
pulse_variables
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
bctmp | Input | real | 320.0 | first wall bulk coolant temperature (C) |
dtstor | Input | real | 300.0 | maximum allowable temperature change in stainless steel thermal storage block (K) ( |
istore | Input | integer | 1 | Switch for thermal storage method:
|
itcycl | Input | integer | 1 | Switch for first wall axial stress model:
|
lpulse | Output | integer | - | Switch for reactor model:
|
impurity_radiation_module
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
nimp | Parameter | integer | 14 | nimp /14/ FIX : number of ion species in impurity radiation model |
coreradius | Input | real | 0.6 | coreradius /0.6/ : normalised radius defining the 'core' region |
coreradiationfraction | Input | real | 1.0 | coreradiationfraction /1.0/ : fraction of radiation from 'core' region that is subtracted from the loss power fimp(nimp) /1.0,0.1,0.02,0.0,0.0,0.0,0.0,0.0,0.0016,0.0,0.0,0.0,0.0,0.0/ : impurity number density fractions relative to electron density |
fimp | Input | real | [1. 0.1 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. ] | |
imp_label | Input | character | [b'H_' b'He' b'Be' b'C_' b'N_' b'O_' b'Ne' b'Si' b'Ar' b'Fe' b'Ni' b'Kr' b'Xe' b'W_'] | imp_label(nimp) : impurity ion species names:
|
all_array_hotfix_len* | Parameter | integer | 200 | |
impurity_arr_label | Input | character | [b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' ' b' '] | |
impurity_arr_z | Output | integer | - | |
impurity_arr_amass | Output | real | - | |
impurity_arr_frac | Output | real | - | |
impurity_arr_len_tab | Output | integer | - | |
impurity_arr_temp_kev | Output | real | - | |
impurity_arr_lz_wm3 | Output | real | - | |
impurity_arr_zav | Output | real | - | |
toolow | Output | logical | - | Used for reporting error in function pimpden |
ccfe_hcpb_module
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
ip | Output | integer | - | |
ofile | Output | integer | - | |
armour_density | Output | real | - | FW armour density [kg/m3] |
fw_density | Output | real | - | FW density [kg/m3] |
blanket_density | Output | real | - | Blanket density [kg/m3] |
shield_density | Output | real | - | Shield density [kg/m3] |
vv_density | Output | real | - | Vacuum vessel density [kg/m3] |
x_blanket | Output | real | - | Blanket exponent (tonne/m2) |
x_shield | Output | real | - | Shield exponent (tonne/m2) |
tfc_nuc_heating | Output | real | - | Unit nuclear heating in TF coil (W per W of fusion power) |
fw_armour_u_nuc_heating | Output | real | - | Unit heating of FW and armour in FW armour (W/kg per W of fusion power) |
shld_u_nuc_heating | Output | real | - | Unit nuclear heating in shield (W per W of fusion power) |
pnuc_tot_blk_sector | Output | real | - | Total nuclear power deposited in blanket covered sector (FW, BLKT, SHLD, TF) (MW) |
exp_blanket | Output | real | - | Exponential factors in nuclear heating calcs |
exp_shield1 | Output | real | - | Exponential factors in nuclear heating calcs |
exp_shield2 | Output | real | - | Exponential factors in nuclear heating calcs |
fson_string_m
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
block_size* | Parameter | integer | 32 |
fson_value_m
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
type_unknown | Parameter | integer | -1 | |
type_null | Parameter | integer | 0 | |
type_object | Parameter | integer | 1 | |
type_array | Parameter | integer | 2 | |
type_string | Parameter | integer | 3 | |
type_integer | Parameter | integer | 4 | |
type_real | Parameter | integer | 5 | |
type_logical | Parameter | integer | 6 |
fson_library
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
end_of_file* | Parameter | integer | -1 | |
end_of_record* | Parameter | integer | -2 | |
state_looking_for_value* | Parameter | integer | 1 | |
state_in_object* | Parameter | integer | 2 | |
state_in_pair_name* | Parameter | integer | 3 | |
state_in_pair_value* | Parameter | integer | 4 | |
pushed_index* | Variable | integer | - | |
pushed_char* | Variable | character | - |
const_and_precisions
Name | Type | Datatype | Default Value | Description |
---|---|---|---|---|
sp_ | Parameter | integer | 4 | |
dp_ | Parameter | integer | 8 | |
wp_ | Parameter | integer | dp_ | |
odep_ | Parameter | integer | dp_ | |
ypi_ | Parameter | integer | 4 | |
ypd_ | Parameter | integer | 8 | |
lfn_ | Parameter | integer | 256 | |
zero | Parameter | real | 0.0_wp_ | |
unit | Parameter | real | 1.0_wp_ | |
pi | Parameter | real | 3.141592653589793_wp_ | |
sqrt_pi | Parameter | real | 1.772453850905516_wp_ | |
sqrt_2 | Parameter | real | 1.414213562373095_wp_ | |
rad | Parameter | real | pi/180.0_wp_ | |
ex | Parameter | real | (/unit, zero, zero/) | |
ey | Parameter | real | (/zero, unit, zero/) | |
ez | Parameter | real | (/zero, zero, unit/) | |
kron | Parameter | real | reshape((/unit, zero, zero, zero, unit, zero, zero, zero, unit/), (/3, 3/)) | |
im | Parameter | complex | (0.0_wp_, 1.0_wp_) | |
czero | Parameter | complex | (0.0_wp_, 0.0_wp_) | |
cunit | Parameter | complex | (1.0_wp_, 0.0_wp_) | |
ctwo | Parameter | complex | (2.0_wp_, 0.0_wp_) | |
comp_eps | Parameter | real | EPSILON(unit) | |
comp_eps2 | Parameter | real | comp_eps**2 | |
comp_tiny | Parameter | real | TINY(unit) | |
comp_huge | Parameter | real | HUGE(unit) | |
comp_tinylog | Parameter | real | -200 | |
comp_hugelog | Parameter | real | +200 | |
output_tiny | Parameter | real | 1.0d-66 | |
output_huge | Parameter | real | 1.0d+66 | |
e_ | Parameter | real | 1.601917d-19 | |
me_ | Parameter | real | 9.109558d-31 | |
mp_ | Parameter | real | 1.672614d-27 | |
rmpe_ | Parameter | real | mp_/me_ | |
c_ | Parameter | real | 2.997925d+08 | |
eps0_ | Parameter | real | 8.854188d-12 | |
kev_ | Parameter | real | 1000*e_ | |
mc2_si | Parameter | real | me_*c_**2 | |
mc2_ | Parameter | real | mc2_SI/keV_ | |
mc_ | Parameter | real | me_*c_ | |
wce1_ | Parameter | real | e_/me_ | |
fce1_ | Parameter | real | wce1_/(2*pi) | |
wpe1_ | Parameter | real | 56.4049201 | |
fpe1_ | Parameter | real | wpe1_/(2*pi) | |
wpe12_ | Parameter | real | wpe1_**2 | |
vte1_ | Parameter | real | 1.8755328e7 | |
curr1_ | Parameter | real | e_*vte1_ | |
umax_ | Parameter | real | 7.0d0 | |
nu_ | Parameter | integer | 700 | |
npar_min | Parameter | real | 1.0d-2 |