

subroutine condcalc(curth,condrate, fii, dp)
double precision a3,a2,a1,a0
if  (fii .eq. 0 ) condrate=1
if (fii .eq. 0.005 .and. dp .eq. 35) then
T=curth*(55-26)+26
a3=2.016536*0.0000001*T**3
a2=-2.722323*0.00001*T**2
a1=0.00220488*T**1
a0=1.04753176
condrate=a0+a2+a3+a1
end if
if (fii .eq. 0.01 .and. dp .eq. 35) then
T=curth*(55-26)+26
a3=2.01654*0.0000001*T**3
a2=-2.72232*0.00001*T**2
a1=0.00220488*T**1
a0=1.01753176
condrate=a0+a2+a3+a1
end if
if (fii .eq. 0.02 .and. dp .eq. 35) then
T=curth*(55-26)+26
a3=-3.62*0.000001*T**3
a2=0.000490018*T**2
a1=-0.01968784*T**1
a0=1.364428312
condrate=a0+a2+a3+a1
end if

if (fii .eq. 0.03 .and. dp .eq. 35) then
T=curth*(55-26)+26
a3=-1.282517*0.0001*T**3
a2=1.731397*0.01*T**2
a1=-7.623037*0.1*T**1
a0=1.216980*10
condrate=a0+a2+a3+a1
end if
if (fii .eq. 0.005 .and. dp .eq. 65) then
T=curth*(55-26)+26
a3=2.016536*0.0000001*T**3
a2=-2.722323*0.00001*T**2
a1=2.204880*0.001*T**1
a0=1.027532
condrate=a0+a2+a3+a1
end if
if (fii .eq. 0.01 .and. dp .eq. 65) then
T=curth*(55-26)+26
a3=1.437790*0.00001*T**3
a2=-1.641016*0.001*T**2
a1=6.190795*0.01*T**1
a0=3.470145*0.1
condrate=a0+a2+a3+a1
end if
if (fii .eq. 0.02 .and. dp .eq. 65) then
T=curth*(55-26)+26
a3=2.621496*0.000001*T**3
a2=-3.539020*0.0001*T**2
a1=1.866344*0.01*T**1
a0=8.879129*0.1
condrate=a0+a2+a3+a1
end if
if (fii .eq. 0.03 .and. dp .eq. 65) then
T=curth*(55-26)+26
a3=8.177052*0.000001*T**3
a2=-1.053902*0.001*T**2
a1=4.685788*0.01*T**1
a0=6.004129*0.1
condrate=a0+a2+a3+a1
end if
if (fii .eq. 0.005 .and. dp .eq. 90) then
T=curth*(55-26)+26
a3=1.109095*0.000001*T**3
a2=-1.497278*0.0001*T**2
a1=7.126840*0.001*T**1
a0=9.264247*0.1
condrate=a0+a2+a3+a1
end if
if (fii .eq. 0.01 .and. dp .eq. 90) then
T=curth*(55-26)+26
a3=5.666465*0.000001*T**3
a2=-6.649728*0.0001*T**2
a1=2.685713*0.01*T**1
a0=7.716425*0.1
condrate=a0+a2+a3+a1
end if
if (fii .eq. 0.02 .and. dp .eq. 90) then
T=curth*(55-26)+26
a3=-6.674733*0.000001*T**3
a2=8.010889*0.0001*T**2
a1=-2.788153*0.01*T**1
a0=1.540699
condrate=a0+a2+a3+a1
end if
if (fii .eq. 0.03 .and. dp .eq. 90) then
T=curth*(55-26)+26
a3=-3.226457*0.000001*T**3
a2=4.355717*0.0001*T**2
a1=-1.527808*0.01*T**1
a0=1.459492
condrate=a0+a2+a3+a1
end if
return
end




Subroutine viscalc(curth, fii, dp, viscorate)
double precision a3,a2,a1,a0
if (fii .eq. 0) viscorate=16.2
if (fii .eq. 0.005 .and. dp .eq. 35) then
T=curth*(55-26)+26
a3=-1.02218*0.001*T**3
a2=1.389946*0.1*T**2
a1=-6.527537*T**1
a0=1.142215*100
viscorate=a0+a2+a3+a1
end if

if (fii .eq. 0.01 .and. dp .eq. 35) then
T=curth*(55-26)+26
a3=-0.00129*T**3
a2=0.172089*T**2
a1=-7.88292*T**1
a0=132.7237
viscorate=a0+a2+a3+a1
end if
if (fii .eq. 0.02 .and. dp .eq. 35) then
T=curth*(55-26)+26
a3=-1.437992*0.001*T**3
a2=1.931289*0.1*T**2
a1=-8.871999*T**1
a0=1.493910*100
viscorate=a0+a2+a3+a1
end if
if (fii .eq. 0.03 .and. dp .eq. 35) then
T=curth*(55-26)+26
a3=-1.777576*0.001*T**3
a2=2.349728*0.1*T**2
a1=-1.058102*10*T**1
a0=1.735075*100
viscorate=a0+a2+a3+a1
end if
if (fii .eq. 0.005 .and. dp .eq. 65) then
T=curth*(55-26)+26
a3=-1.08066*0.001*T**3
a2=1.468893*0.1*T**2
a1=-6.866952*T**1
a0=1.87373*100
viscorate=a0+a2+a3+a1
end if
if (fii .eq. 0.01 .and. dp .eq. 65) then
T=curth*(55-26)+26
a3=-1.285844*0.001*T**3
a2=1.720889*0.1*T**2
a1=-7.882917*T**1
a0=1.327237*100
viscorate=a0+a2+a3+a1
end if
if (fii .eq. 0.02 .and. dp .eq. 65) then
T=curth*(55-26)+26
a3=-1.327788*0.001*T**3
a2=1.797714*0.1*T**2
a1=-8.343532*T**1
a0=1.422571*100
viscorate=a0+a2+a3+a1
end if
if (fii .eq. 0.03 .and. dp .eq. 65) then
T=curth*(55-26)+26
a3=-1.480339*0.001*T**3
a2=1.988457*0.1*T**2
a1=-9.135024*T**1
a0=1.536093*100
viscorate=a0+a2+a3+a1
end if

if (fii .eq. 0.005 .and. dp .eq. 90) then
T=curth*(55-26)+26
a3=-1.202057*0.001*T**3
a2=1.612777*0.1*T**2
a1=-7.392290*T**1
a0=1.243032*100
viscorate=a0+a2+a3+a1
end if
if (fii .eq. 0.01 .and. dp .eq. 90) then
T=curth*(55-26)+26
a3=-1.084190*0.001*T**3
a2=1.448657*0.1*T**2
a1=-6.678037*T**1
a0=1.152555*100
viscorate=a0+a2+a3+a1
end if
if (fii .eq. 0.02 .and. dp .eq. 90) then
T=curth*(55-26)+26
a3=-1.327788*0.001*T**3
a2=1.797514*0.1*T**2
a1=-8.343532*T**1
a0=1.422571*100
viscorate=a0+a2+a3+a1
end if
if (fii .eq. 0.03 .and. dp .eq. 90) then
T=curth*(55-26)+26
a3=-1.480339*0.001*T**3
a2=1.988457*0.1*T**2
a1=-9.135024*T**1
a0=1.536093*100
viscorate=a0+a2+a3+a1
end if
return
end