Initialise fwbs variables
subroutine init_fwbs_variables
!! Initialise fwbs variables
implicit none
bktlife = 0.0D0
bktlife_cal = 0.0D0
coolmass = 0.0D0
vvmass = 0.0D0
denstl = 7800.0D0
denw = 19250.0D0
denwc = 15630.0D0
dewmkg = 0.0D0
emult = 1.269D0
emultmw = 0.0D0
fblss = 0.09705D0
fdiv = 0.115D0
fhcd = 0.0D0
fhole = 0.0D0
fwbsshape = 2
fwlife = 0.0D0
fwmass = 0.0D0
fw_armour_mass = 0.0D0
fw_armour_thickness = 0.005D0
fw_armour_vol = 0.0D0
iblanket = 1
iblnkith = 1
inuclear = 0
qnuc = 0.0D0
li6enrich = 30.0D0
pnucblkt = 0.0D0
pnucdiv = 0.0D0
pnucfw = 0.0D0
pnuchcd = 0.0D0
pnucloss = 0.0D0
pnucvvplus = 0.0D0
pnucshld = 0.0D0
whtblkt = 0.0D0
whtblss = 0.0D0
armour_fw_bl_mass = 0.0D0
breeder_f = 0.5D0
breeder_multiplier = 0.75D0
vfcblkt = 0.05295D0
vfpblkt = 0.1D0
whtblli4sio4 = 0.0D0
whtbltibe12 = 0.0D0
f_neut_shield = -1.0D0
vffwi = 0.0D0
vffwo = 0.0D0
psurffwi = 0.0D0
psurffwo = 0.0D0
volfw = 0.0D0
fblss_ccfe = 0.0D0
fblli2sio4 = 0.0D0
fbltibe12 = 0.0D0
breedmat = 1
densbreed = 0.0D0
fblbe = 0.6D0
fblbreed = 0.154D0
fblhebmi = 0.4D0
fblhebmo = 0.4D0
fblhebpi = 0.6595D0
fblhebpo = 0.6713D0
hcdportsize = 1
nflutf = 0.0D0
npdiv = 2
nphcdin = 2
nphcdout = 2
tbr = 0.0D0
tritprate = 0.0D0
wallpf = 1.21D0
whtblbreed = 0.0D0
whtblbe = 0.0D0
iblanket_thickness = 2
primary_pumping = 2
i_shield_mat = 0
secondary_cycle = 0
secondary_cycle_liq = 4
coolwh = 1
afwi = 0.008D0
afwo = 0.008D0
fwcoolant = 'helium'
fw_wall = 0.003D0
afw = 0.006D0
pitch = 0.02D0
fwinlet = 573.0D0
fwoutlet = 823.0D0
fwpressure = 15.5D6
tpeak = 873.0D0
roughness = 1.0D-6
fw_channel_length = 4.0D0
peaking_factor = 1.0D0
blpressure = 15.50D6
inlet_temp = 573.0D0
outlet_temp = 823.0D0
coolp = 15.5D6
nblktmodpo = 8
nblktmodpi = 7
nblktmodto = 48
nblktmodti = 32
tfwmatmax = 823.0D0
fw_th_conductivity = 28.34D0
fvoldw = 1.74D0
fvolsi = 1.0D0
fvolso = 0.64D0
fwclfr = 0.15D0
praddiv = 0.0D0
pradfw = 0.0D0
pradhcd = 0.0D0
pradloss = 0.0D0
ptfnuc = 0.0D0
ptfnucpm3 = 0.0D0
rdewex = 0.0D0
zdewex = 0.0D0
rpf2dewar = 0.5D0
vdewex = 0.0D0
vdewin = 0.0D0
vfshld = 0.25D0
volblkt = 0.0D0
volblkti = 0.0D0
volblkto = 0.0D0
volshld = 0.0D0
whtshld = 0.0D0
wpenshld = 0.0D0
wtshldi = 0.0D0
wtshldo = 0.0D0
irefprop = 1
fblli = 0.0D0
fblli2o = 0.08D0
fbllipb = 0.68D0
fblvd = 0.0D0
wtblli2o = 0.0D0
wtbllipb = 0.0D0
whtblvd = 0.0D0
whtblli = 0.0D0
vfblkt = 0.25D0
blktmodel = 0
declblkt = 0.075D0
declfw = 0.075D0
declshld = 0.075D0
blkttype = 3
etaiso = 0.85D0
etahtp = 0.95D0
pnuc_cp = 0.0D0
pnuc_cp_sh = 0.0D0
pnuc_cp_tf = 0.0D0
neut_flux_cp = 0.0D0
ipump = 0
i_bb_liq = 0
icooldual = 0
ifci = 0
ims = 0
n_liq_recirc = 10
r_f_liq_ib=0.5
r_f_liq_ob=0.5
w_f_liq_ib=0.5
w_f_liq_ob=0.5
den_ceramic = 3.21D3
th_wall_secondary = 1.25D-2
bz_channel_conduct_liq = 8.33D5
a_bz_liq = 0.2D0
b_bz_liq = 0.2D0
nopol = 2
nopipes = 4
den_liq = 9.5D3
specific_heat_liq = 1.9D2
thermal_conductivity_liq = 30.0
wht_liq = 0.0D0
wht_liq_ib = 0.0D0
wht_liq_ob = 0.0D0
dynamic_viscosity_liq = 0.0D0
electrical_conductivity_liq = 0.0D0
hartmann_liq = (0.0D0, 0.0D0)
b_mag_blkt = (5.0D0, 5.0D0)
etaiso_liq = 0.85D0
blpressure_liq = 1.7D6
inlet_temp_liq = 570.0D0
outlet_temp_liq = 720.0D0
rhof_fw = 0.0D0
visc_fw = 0.0D0
rhof_bl = 0.0D0
visc_bl = 0.0D0
cp_fw = 0.0D0
cv_fw = 0.0D0
cp_bl = 0.0D0
cv_bl = 0.0D0
f_nuc_pow_bz_struct = 0.34D0
f_nuc_pow_bz_liq = 0.66D0
pnuc_fw_ratio_dcll = 0.14D0
pnuc_blkt_ratio_dcll = 0.86D0
bzfllengi_n_rad = 4
bzfllengi_n_pol = 2
bzfllengo_n_rad = 4
bzfllengo_n_pol = 2
bzfllengi_n_rad_liq = 2
bzfllengi_n_pol_liq = 2
bzfllengo_n_rad_liq = 2
bzfllengo_n_pol_liq = 2
end subroutine init_fwbs_variables