61 real :: rho1(2),tt1(2),temperature,airdens,vis_dyn,vis_kin,
viscosity
62 real :: settling,settling_old,reynolds,c_d
63 integer :: i,n,nix,njy,indzh,nsp
77 if (height(i).gt.zt)
then
88 dz=1./(height(indz+1)-height(indz))
89 dz1=(zt-height(indz))*dz
90 dz2=(height(indz+1)-zt)*dz
101 rho1(n)=rho(nix,njy,indzh,1)
102 tt1(n)=tt(nix,njy,indzh,1)
109 temperature=dz2*tt1(1)+dz1*tt1(2)
110 airdens=dz2*rho1(1)+dz1*rho1(2)
114 vis_kin=vis_dyn/airdens
116 reynolds=dquer(nsp)/1.e6*abs(vsetaver(nsp))/vis_kin
123 settling_old=vsetaver(nsp)
127 if (reynolds.lt.1.917)
then
129 else if (reynolds.lt.500.)
then
130 c_d=18.5/(reynolds**0.6)
136 sqrt(4*ga*dquer(nsp)/1.e6*density(nsp)*cunningham(nsp)/ &
139 if (abs((settling-settling_old)/settling).lt.0.01) goto 11
141 reynolds=dquer(nsp)/1.e6*abs(settling)/vis_kin
142 settling_old=settling
real function viscosity(t)
subroutine get_settling(itime, xt, yt, zt, nsp, settling)