Initialise module variables
subroutine init_physics_variables
!! Initialise module variables
implicit none
m_beam_amu = 0.0D0
m_fuel_amu = 0.0D0
m_ions_total_amu = 0.0D0
alphaj = 1.0D0
alphan = 0.25D0
alphap = 0.0D0
alpha_rate_density_total = 0.0D0
alpha_rate_density_plasma = 0.0D0
alphat = 0.5D0
aspect = 2.907D0
beamfus0 = 1.0D0
beta = 0.042D0
beta_fast_alpha = 0.0D0
beta_max = 0.0D0
beta_min = 0.0D0
beta_beam = 0.0D0
beta_poloidal = 0.0D0
beta_poloidal_eps = 0.0D0
beta_toroidal = 0.0D0
beta_thermal = 0.0D0
beta_thermal_poloidal = 0.0D0
beta_thermal_toroidal = 0.0D0
beta_norm_total = 0.0D0
beta_norm_thermal = 0.0D0
beta_norm_poloidal = 0.0D0
e_plasma_beta_thermal = 0.0D0
beta_norm_toroidal = 0.0D0
betbm0 = 1.5D0
bp = 0.0D0
bt = 5.68D0
btot = 0.0D0
burnup = 0.0D0
burnup_in = 0.0D0
bvert = 0.0D0
c_beta = 0.5D0
csawth = 1.0D0
f_vol_plasma = 1.0D0
f_r_conducting_wall = 1.35D0
dene = 9.8D19
nd_fuel_ions = 0.0D0
dlamee = 0.0D0
dlamie = 0.0D0
dlimit = 0.0D0
nd_alphas = 0.0D0
nd_beam_ions = 0.0D0
beam_density_out = 0.0D0
beta_norm_max = 3.5D0
dnelimt = 0.0D0
nd_ions_total = 0.0D0
dnla = 0.0D0
nd_protons = 0.0D0
dntau = 0.0D0
nd_impurities = 0.0D0
beta_poloidal_eps_max = 1.38D0
eps = 0.34399724802D0
aux_current_fraction = 0.0D0
inductive_current_fraction = 0.0D0
f_alpha_electron = 0.0D0
f_alpha_plasma = 0.95D0
f_alpha_ion = 0.0D0
f_deuterium = 0.5D0
ftar = 1.0D0
ffwal = 0.92D0
fgwped = 0.85D0
fgwsep = 0.50D0
f_helium3 = 0.0D0
figmer = 0.0D0
fkzohm = 1.0D0
fplhsep = 1.0D0
fpdivlim = 1.0D0
fne0 = 1.0D0
f_tritium = 0.5D0
fusion_rate_density_total = 0.0D0
fusion_rate_density_plasma = 0.0D0
fvsbrnni = 1.0D0
gamma = 0.4D0
f_beta_alpha_beam_thermal = 0.0D0
hfac = 0.0D0
hfact = 1.0D0
taumax = 10.0D0
i_bootstrap_current = 3
i_beta_component = 0
i_plasma_current = 4
i_diamagnetic_current = 0
i_density_limit = 8
idivrt = 2
i_beta_fast_alpha = 1
ignite = 0
iinvqd = 1
ipedestal = 1
i_pfirsch_schluter_current = 0
neped = 4.0D19
nesep = 3.0D19
alpha_crit = 0.0D0
nesep_crit = 0.0D0
plasma_res_factor = 1.0D0
rhopedn = 1.0D0
rhopedt = 1.0D0
tbeta = 2.0D0
teped = 1.0D0
tesep = 0.1D0
iprofile = 1
iradloss = 1
isc = 34
i_plasma_wall_gap = 1
i_plasma_geometry = 0
i_plasma_shape = 0
itart = 0
itartpf = 0
iwalld = 1
plasma_square = 0.0D0
kappa = 1.792D0
kappa95 = 1.6D0
kappaa = 0.0D0
kappaa_IPB = 0.d0
ne0 = 0.0D0
ni0 = 0.0D0
m_s_limit = 0.3D0
p0 = 0.0D0
j_plasma_0 = 0.0D0
f_dd_branching_trit = 0.0D0
alpha_power_density_plasma = 0.0D0
alpha_power_density_total = 0.0D0
alpha_power_electron_density = 0.0D0
palpfwmw = 0.0D0
alpha_power_ions_density = 0.0D0
alpha_power_total = 0.0D0
alpha_power_plasma = 0.0D0
alpha_power_beams = 0.0D0
non_alpha_charged_power = 0.0D0
charged_power_density = 0.0D0
pcoef = 0.0D0
p_plasma_inner_rad_mw = 0.0D0
pcoreradpv = 0.0D0
dd_power = 0.0D0
dhe3_power = 0.0D0
pdivt = 0.0D0
pdivl = 0.0D0
pdivu = 0.0D0
pdivmax = 0.0D0
dt_power_total = 0.0D0
dt_power_plasma = 0.0D0
p_plasma_outer_rad_mw = 0.0D0
pedgeradpv = 0.0D0
charged_particle_power = 0.0D0
phiint = 0.0D0
pflux_fw_rad_mw = 0.0D0
piepv = 0.0D0
plasma_current = 0.0D0
neutron_power_plasma = 0.0D0
neutron_power_total = 0.0D0
neutron_power_density_total = 0.0D0
neutron_power_density_plasma = 0.0D0
p_plasma_ohmic_mw = 0.0D0
pden_plasma_ohmic_mw = 0.0D0
powerht = 0.0D0
fusion_power = 0.0D0
len_plasma_poloidal = 0.0D0
p_plasma_rad_mw = 0.0D0
pden_plasma_rad_mw = 0.0D0
pradsolmw = 0.0D0
proton_rate_density = 0.0D0
psolradmw = 0.0D0
pden_plasma_sync_mw = 0.0D0
ilhthresh = 19
plhthresh = 0.0D0
pthrmw = 0.0D0
ptremw = 0.0D0
ptrepv = 0.0D0
ptrimw = 0.0D0
pscalingmw = 0.0D0
ptripv = 0.0D0
q = 3.0D0
q0 = 1.0D0
q95 = 0.0D0
qfuel = 0.0D0
tauratio = 1.0D0
q95_min = 0.0D0
qstar = 0.0D0
rad_fraction_sol = 0.8D0
rad_fraction_total = 0.0D0
f_nd_alpha_electron = 0.10D0
f_nd_protium_electrons = 0.0D0
rli = 0.9D0
rlp = 0.0D0
rmajor = 8.14D0
rminor = 0.0D0
f_nd_beam_electron = 0.005D0
rncne = 0.0D0
rndfuel = 0.0D0
rnfene = 0.0D0
rnone = 0.0D0
rpfac = 0.0D0
res_plasma = 0.0D0
res_time = 0.0D0
a_plasma_surface = 0.0D0
a_plasma_surface_outboard = 0.0D0
i_single_null = 1
f_sync_reflect = 0.6D0
tauee = 0.0D0
tauee_in = 0.0D0
taueff = 0.0D0
tauei = 0.0D0
taup = 0.0D0
te = 12.9D0
te0 = 0.0D0
ten = 0.0D0
ti = 12.9D0
ti0 = 0.0D0
tin = 0.0D0
tratio = 1.0D0
triang = 0.36D0
triang95 = 0.24D0
vol_plasma = 0.0D0
vsbrn = 0.0D0
vshift = 0.0D0
vsind = 0.0D0
vsres = 0.0D0
vsstt = 0.0D0
wallmw = 0.0D0
wtgpd = 0.0D0
a_plasma_poloidal = 0.0D0
zeff = 0.0D0
zeffai = 0.0D0
end subroutine init_physics_variables