22 subroutine getrb(nc,ustar,nyl,diffh2o,reldiff,rb)
49 real :: ustar,diffh2o,rb(maxspec),schmidt,nyl
50 real :: reldiff(maxspec)
52 real,
parameter :: pr=0.72
55 if (reldiff(ic).gt.0.)
then
56 schmidt=nyl/diffh2o*reldiff(ic)
57 rb(ic)=2.0*(schmidt/pr)**0.67/(karman*ustar)
subroutine getrb(nc, ustar, nyl, diffh2o, reldiff, rb)