Initialise fwbs variables
subroutine init_fwbs_variables
!! Initialise fwbs variables
implicit none
life_blkt_fpy = 0.0D0
life_blkt = 0.0D0
m_fw_blkt_div_coolant_total = 0.0D0
m_vv = 0.0D0
denstl = 7800.0D0
denwc = 15630.0D0
dewmkg = 0.0D0
emult = 1.269D0
emultmw = 0.0D0
fblss = 0.09705D0
f_ster_div_single = 0.115D0
f_a_fw_hcd = 0.0D0
fhole = 0.0D0
i_fw_blkt_vv_shape = 2
life_fw_fpy = 0.0D0
m_fw_total = 0.0D0
fw_armour_mass = 0.0D0
fw_armour_thickness = 0.005D0
fw_armour_vol = 0.0D0
i_blanket_type = 1
i_blkt_inboard = 1
inuclear = 0
qnuc = 0.0D0
f_blkt_li6_enrichment = 30.0D0
p_blkt_nuclear_heat_total_mw = 0.0D0
p_div_nuclear_heat_total_mw = 0.0D0
p_fw_nuclear_heat_total_mw = 0.0D0
p_fw_hcd_nuclear_heat_mw = 0.0D0
pnucloss = 0.0D0
pnucvvplus = 0.0D0
pnucshld = 0.0D0
m_blkt_total = 0.0D0
m_blkt_steel_total = 0.0D0
armour_fw_bl_mass = 0.0D0
breeder_f = 0.5D0
breeder_multiplier = 0.75D0
vfcblkt = 0.05295D0
vfpblkt = 0.1D0
m_blkt_li4sio4 = 0.0D0
m_blkt_tibe12 = 0.0D0
f_neut_shield = -1.0D0
f_a_fw_coolant_inboard = 0.0D0
f_a_fw_coolant_outboard = 0.0D0
psurffwi = 0.0D0
psurffwo = 0.0D0
vol_fw_total = 0.0D0
f_vol_blkt_steel = 0.0D0
f_vol_blkt_li4sio4 = 0.0D0
f_vol_blkt_tibe12 = 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
m_blkt_beryllium = 0.0D0
i_coolant_pumping = 2
i_shield_mat = 0
i_thermal_electric_conversion = 0
secondary_cycle_liq = 4
i_blkt_coolant_type = 1
i_fw_coolant_type = 'helium'
dr_fw_wall = 0.003D0
radius_fw_channel = 0.006D0
dx_fw_module = 0.02D0
temp_fw_coolant_in = 573.0D0
temp_fw_coolant_out = 823.0D0
pres_fw_coolant = 15.5D6
temp_fw_peak = 873.0D0
roughness = 1.0D-6
len_fw_channel = 4.0D0
f_fw_peak = 1.0D0
pres_blkt_coolant = 15.50D6
temp_blkt_coolant_in = 573.0D0
temp_blkt_coolant_out = 823.0D0
coolp = 15.5D6
n_blkt_outboard_modules_poloidal = 8
n_blkt_inboard_modules_poloidal = 7
n_blkt_outboard_modules_toroidal = 48
n_blkt_inboard_modules_toroidal = 32
temp_fw_max = 823.0D0
fw_th_conductivity = 28.34D0
fvoldw = 1.74D0
fvolsi = 1.0D0
fvolso = 0.64D0
fwclfr = 0.15D0
p_div_rad_total_mw = 0.0D0
p_fw_rad_total_mw = 0.0D0
p_fw_hcd_rad_total_mw = 0.0D0
pradloss = 0.0D0
ptfnuc = 0.0D0
ptfnucpm3 = 0.0D0
r_cryostat_inboard = 0.0D0
z_cryostat_half_inside = 0.0D0
dr_pf_cryostat = 0.5D0
vol_cryostat = 0.0D0
vol_cryostat_internal = 0.0D0
vol_vv = 0.0D0
vfshld = 0.25D0
vol_blkt_total = 0.0D0
vol_blkt_inboard = 0.0D0
vol_blkt_outboard = 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
m_blkt_li2o = 0.0D0
wtbllipb = 0.0D0
m_blkt_vanadium = 0.0D0
m_blkt_lithium = 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
i_fw_blkt_shared_coolant = 0
i_blkt_liquid_breeder_type = 0
i_blkt_dual_coolant = 0
i_blkt_liquid_breeder_channel_type = 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
den_fw_coolant = 0.0D0
visc_fw_coolant = 0.0D0
den_blkt_coolant = 0.0D0
visc_blkt_coolant = 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