Initialise module variables
subroutine init_ife_variables
!! Initialise module variables
implicit none
bldr = 1.0D0
bldrc = 1.0D0
bldzl = 4.0D0
bldzu = 4.0D0
blmatf = reshape( (/ &
0.05D0,0.05D0,0.05D0, &
0.0D0, 0.0D0, 0.0D0, &
0.45D0,0.45D0,0.45D0, &
0.0D0, 0.0D0, 0.0D0, &
0.20D0,0.20D0,0.20D0, &
0.0D0, 0.0D0, 0.0D0, &
0.30D0,0.30D0,0.30D0, &
0.0D0,0.0D0,0.0D0, &
0.0D0, 0.0D0, 0.0D0 /), shape(blmatf))
blmatm = 0.0D0
blmatv = 0.0D0
blvol = 0.0D0
cdriv0 = 154.3D0
cdriv1 = 163.2D0
cdriv2 = 244.9D0
cdriv3 = 1.463D0
chdzl = 9.0D0
chdzu = 9.0D0
chmatf = &
(/1.0D0,0.0D0,0.0D0,0.0D0,0.0D0,0.0D0,0.0D0,0.0D0,0.0D0/)
chmatm = 0.0D0
chmatv = 0.0D0
chrad = 6.5D0
chvol = 0.0D0
dcdrv0 = 111.4D0
dcdrv1 = 78.0D0
dcdrv2 = 59.9D0
drveff = 0.28D0
edrive = 5.0D6
etadrv = 0.0D0
etali = 0.4D0
etave = (/ &
0.082D0,0.079D0,0.076D0,0.073D0,0.069D0, &
0.066D0,0.062D0,0.059D0,0.055D0,0.051D0 /)
fauxbop = 0.06D0
fbreed = 0.51D0
fburn = 0.3333D0
flirad = 0.78D0
frrmax = 1.0D0
fwdr = 0.01D0
fwdzl = 0.01D0
fwdzu = 0.01D0
fwmatf = reshape( (/ &
0.05D0,0.05D0,0.05D0, &
0.0D0, 0.0D0, 0.0D0, &
0.95D0,0.95D0,0.95D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0 /), shape(fwmatf))
fwmatm = 0.0D0
fwmatv = 0.0D0
fwvol = 0.0D0
gain = 0.0D0
gainve = (/ &
60.0D0, 95.0D0,115.0D0,125.0D0,133.0D0, &
141.0D0,152.0D0,160.0D0,165.0D0,170.0D0 /)
htpmw_ife = 0.0D0
ife = 0
ifedrv = 2
ifetyp = 0
lipmw = 0.0D0
mcdriv = 1.0D0
mflibe = 0.0D0
pdrive = 23.0D6
pfusife = 1000.0D0
pifecr = 10.0D0
ptargf = 2.0D0
r1 = 0.0D0
r2 = 0.0D0
r3 = 0.0D0
r4 = 0.0D0
r5 = 0.0D0
r6 = 0.0D0
r7 = 0.0D0
reprat = 0.0D0
rrin = 6.0D0
rrmax = 20.0D0
shdr = 1.7D0
shdzl = 5.0D0
shdzu = 5.0D0
shmatf = reshape( (/ &
0.05D0,0.05D0,0.05D0, &
0.19D0,0.19D0,0.19D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.665D0,0.665D0,0.665D0, &
0.095D0,0.095D0,0.095D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0 /), shape(shmatf))
shmatm = 0.0D0
shmatv = 0.0D0
shvol = 0.0D0
sombdr = 2.7D0
somtdr = 2.7D0
taufall = 0.0D0
tdspmw = 0.01D0
tfacmw = 0.0D0
tgain = 85.0D0
uccarb = 50.0D0
ucconc = 0.1D0
ucflib = 84.0D0
uctarg = 0.3D0
v1dr = 0.0D0
v1dzl = 0.0D0
v1dzu = 0.0D0
v1matf = reshape( (/ &
1.0D0, 1.0D0, 1.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0 /), shape(v1matf))
v1matm = 0.0D0
v1matv = 0.0D0
v1vol = 0.0D0
v2dr = 2.0D0
v2dzl = 7.0D0
v2dzu = 7.0D0
v2matf = reshape( (/ &
1.0D0, 1.0D0, 1.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0 /), shape(v2matf))
v2matm = 0.0D0
v2matv = 0.0D0
v2vol = 0.0D0
v3dr = 43.3D0
v3dzl = 30.0D0
v3dzu = 20.0D0
v3matf = reshape( (/ &
1.0D0, 1.0D0, 1.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0, &
0.0D0, 0.0D0, 0.0D0 /), shape(v3matf))
v3matm = 0.0D0
v3matv = 0.0D0
v3vol = 0.0D0
zl1 = 0.0D0
zl2 = 0.0D0
zl3 = 0.0D0
zl4 = 0.0D0
zl5 = 0.0D0
zl6 = 0.0D0
zl7 = 0.0D0
zu1 = 0.0D0
zu2 = 0.0D0
zu3 = 0.0D0
zu4 = 0.0D0
zu5 = 0.0D0
zu6 = 0.0D0
zu7 = 0.0D0
end subroutine init_ife_variables