22 subroutine getrc(nc,i,j,t,gr,rh,rr,rc)
59 real :: gr,rh,rr,t,rs,rsm,corr,rluc,rclc,rgsc,rdc,rluo
67 if ((t.gt.0.).and.(t.lt.40.))
then
68 rs=ri(i,j)*(1.+(200./(gr+0.1))**2)*(400./(t*(40.-t)))
78 if ((rh.gt.0.9).or.(rr.gt.0.)) rs=rs*3.
83 rdc=100.*(1.+1000./(gr+10.))
86 corr=1000.*exp(-1.*t-4.)
88 if (reldiff(ic).gt.0.)
then
93 rsm=rs*reldiff(ic)+rm(ic)
100 rgsc=rgs(ic,i,j)+corr
106 rluo=1./(1./1000.+1./(3.*rluc))
107 rluc=1./(1./(3.*rluc)+1.e-7*henry(ic)+f0(ic)/rluo)
108 else if (rh.gt.0.9)
then
109 rluo=1./(1./3000.+1./(3.*rluc))
110 rluc=1./(1./(3.*rluc)+1.e-7*henry(ic)+f0(ic)/rluo)
116 rc(ic)=1./(1./rsm+1./rluc+1./(rdc+rclc)+1./(rac(i,j)+rgsc))
118 if (rc(ic).lt.10.) rc(ic)=10.
subroutine getrc(nc, i, j, t, gr, rh, rr, rc)