23 l,gr,rh,rr,snow,vdepo,lnest)
60 integer :: yyyymmdd,hhmmss,yyyy,mmdd,n,lseason,i,j,ix,jy,lnest
61 real :: vdepo(maxspec),vd,rb(maxspec),rc(maxspec),raquer,ylat
62 real ::
raerod,ra,ust,temp,tc,pa,l,gr,rh,rr,myl,nyl,rhoa,diffh2o,snow
63 real :: slanduse(numclass)
64 real,
parameter :: eps=1.e-5
70 jul=bdate+
real(wftime(n),kind=dp)/86400._dp
78 call
caldate(jul,yyyymmdd,hhmmss)
80 mmdd=yyyymmdd-10000*yyyy
82 if ((ylat.gt.-20).and.(ylat.lt.20))
then
86 if ((mmdd.ge.1201).or.(mmdd.le.301))
then
88 else if ((mmdd.ge.1101).or.(mmdd.le.331))
then
90 else if ((mmdd.ge.401).and.(mmdd.le.515))
then
92 else if ((mmdd.ge.516).and.(mmdd.le.915))
then
100 diffh2o=2.11e-5*(temp/273.15)**1.94*(101325/pa)
111 myl=(1.718+0.0049*tc-1.2e-05*tc**2)*1.e-05
113 myl=(1.718+0.0049*tc)*1.e-05
134 call
getrb(nspec,ust,nyl,diffh2o,reldiff,rb)
138 if (snow.gt.0.001)
then
145 slanduse(j)=xlandusen(ix,jy,j,lnest)
152 if (slanduse(j).gt.eps)
then
158 raquer=raquer+ra*slanduse(j)
163 call
getrc(nspec,lseason,j,tc,gr,rh,rr,rc)
170 if (reldiff(i).gt.0.)
then
171 if ((ra+rb(i)+rc(i)).gt.0.)
then
172 vd=1./(ra+rb(i)+rc(i))
179 vdepo(i)=vdepo(i)+vd*slanduse(j)
189 call
partdep(nspec,density,fract,schmi,vset,raquer,ust,nyl,vdepo)
197 if ((reldiff(i).lt.0.).and.(density(i).lt.0.).and. &
198 (dryvel(i).gt.0.))
then
real function raerod(l, ust, z0)
subroutine caldate(juldate, yyyymmdd, hhmiss)
subroutine getrc(nc, i, j, t, gr, rh, rr, rc)
subroutine getrb(nc, ustar, nyl, diffh2o, reldiff, rb)
subroutine getvdep_nests(n, ix, jy, ust, temp, pa, L, gr, rh, rr, snow, vdepo, lnest)
subroutine partdep(nc, density, fract, schmi, vset, ra, ustar, nyl, vdep)