sustar(i) = sustar(i) + aitr*ustar(i)
_RL ustar
ustar = sqrt(ustmp/atmrho)
ustar = 0. _d 0
if ( ustar .eq. 0. _d 0 ) then
else if ( ustar .lt. ustofu11 ) then
tmp2 = sqrt(tmp1*tmp1 + ustar*ustar/cquadrag_1)
tmp4 = ustar*ustar/clindrag_1/2 - tmp3**3
tmp5 = sqrt(ustar*ustar/clindrag_1*
& (ustar*ustar/clindrag_1/4 - tmp3**3))
rdn = ustar/us
ustar = rdn*sh
& ustar**2
ustar = rd*sh
tau = atmrho*ustar**2
hs(i,j,bi,bj) = atmcp*tau*tstar/ustar
hl(i,j,bi,bj) = flamb*tau*qstar/ustar
evap(i,j,bi,bj) = -recip_rhonil*tau*qstar/ustar