22 subroutine partdep(nc,density,fract,schmi,vset,ra,ustar,nyl,vdep)
77 real :: density(maxspec),schmi(maxspec,ni),fract(maxspec,ni)
78 real :: vset(maxspec,ni)
79 real :: vdep(maxspec),stokes,vdepj,rdp,ustar,alpha,ra,nyl
80 real,
parameter :: eps=1.e-5
85 if (density(ic).gt.0.)
then
87 if (ustar.gt.eps)
then
92 stokes=vset(ic,j)/ga*ustar*ustar/nyl
98 if (alpha.le.log10(eps))
then
99 rdp=1./(schmi(ic,j)*ustar)
101 rdp=1./((schmi(ic,j)+10.**alpha)*ustar)
103 vdepj=vset(ic,j)+1./(ra+rdp+ra*rdp*vset(ic,j))
111 vdep(ic)=vdep(ic)+vdepj*fract(ic,j)
subroutine partdep(nc, density, fract, schmi, vset, ra, ustar, nyl, vdep)