Initialise module variables
subroutine init_physics_variables
!! Initialise module variables
implicit none
abeam = 0.0D0
afuel = 0.0D0
aion = 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
betaft = 0.0D0
betalim = 0.0D0
betalim_lower = 0.0D0
beta_beam = 0.0D0
betap = 0.0D0
normalised_total_beta = 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
cvol = 1.0D0
cwrmax = 1.35D0
dene = 9.8D19
deni = 0.0D0
dlamee = 0.0D0
dlamie = 0.0D0
dlimit = 0.0D0
dnalp = 0.0D0
dnbeam = 0.0D0
beam_density_out = 0.0D0
dnbeta = 3.5D0
dnelimt = 0.0D0
dnitot = 0.0D0
dnla = 0.0D0
dnprot = 0.0D0
dntau = 0.0D0
dnz = 0.0D0
epbetmax = 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
gammaft = 0.0D0
hfac = 0.0D0
hfact = 1.0D0
taumax = 10.0D0
i_bootstrap_current = 3
iculbl = 0
i_plasma_current = 4
i_diamagnetic_current = 0
i_density_limit = 8
idivrt = 2
ifalphap = 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
iscrp = 1
ishape = 0
itart = 0
itartpf = 0
iwalld = 1
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
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
pinnerzoneradmw = 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
pouterzoneradmw = 0.0D0
pedgeradpv = 0.0D0
charged_particle_power = 0.0D0
phiint = 0.0D0
photon_wall = 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
pohmmw = 0.0D0
pohmpv = 0.0D0
powerht = 0.0D0
fusion_power = 0.0D0
pperim = 0.0D0
pradmw = 0.0D0
pradpv = 0.0D0
pradsolmw = 0.0D0
proton_rate_density = 0.0D0
psolradmw = 0.0D0
psyncpv = 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
qlim = 0.0D0
qstar = 0.0D0
rad_fraction_sol = 0.8D0
rad_fraction_total = 0.0D0
ralpne = 0.10D0
protium = 0.0D0
rli = 0.9D0
rlp = 0.0D0
rmajor = 8.14D0
rminor = 0.0D0
rnbeam = 0.005D0
rncne = 0.0D0
rndfuel = 0.0D0
rnfene = 0.0D0
rnone = 0.0D0
rpfac = 0.0D0
rplas = 0.0D0
res_time = 0.0D0
sarea = 0.0D0
sareao = 0.0D0
sf = 0.0D0
i_single_null = 1
ssync = 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
plasma_volume = 0.0D0
vsbrn = 0.0D0
vshift = 0.0D0
vsind = 0.0D0
vsres = 0.0D0
vsstt = 0.0D0
wallmw = 0.0D0
wtgpd = 0.0D0
xarea = 0.0D0
zeff = 0.0D0
zeffai = 0.0D0
end subroutine init_physics_variables