22 subroutine part0(dquer,dsigma,density,fract,schmi,cun,vsh)
78 real,
parameter :: tr=293.15
81 real :: dquer,dsigma,density,xdummy,d01,d02,delta,x01,x02,fract(ni)
82 real :: dmean,alpha,cun,dc,schmidt,schmi(ni),vsh(ni),kn,
erf
83 real,
parameter :: myl=1.81e-5,nyl=0.15e-4
84 real,
parameter :: lam=6.53e-8,kb=1.38e-23,eps=1.2e-38
90 xdummy=sqrt(2.)*alog(dsigma)
99 d01=dquer*dsigma**(-3)
102 d01=dquer*dsigma**(-3.+delta*
real(i))
103 x01=alog(d01/dquer)/xdummy
104 x02=alog(d02/dquer)/xdummy
110 fract(i)=0.5*(
erf(x01)-
erf(x02))
116 dmean=1.e-6*exp(0.5*alog(d01*d02))
123 if ((-1.1/kn).le.log10(eps)*log(10.))
then
126 alpha=1.257+0.4*exp(-1.1/kn)
129 dc=kb*tr*cun/(3.*pi*myl*dmean)
131 schmi(i)=schmidt**(-2./3.)
132 vsh(i)=ga*density*dmean*dmean*cun/(18.*myl)
subroutine part0(dquer, dsigma, density, fract, schmi, cun, vsh)