re = ren/( exf_one + ren*(ztln-psixh)/karman )
qstar = re*delq
dEvdT = (tau*re)*ssq*qsat_exp/Ts2
_RL re ! = Ce / sqrt(Cd) [-]
re = rh
clha = rhoa*lath*us*re*rd
re = ren/(1. _d 0 + ren*(ztln-psixh)/xkar)
clha = lath*tau*re
dEvdT = tau*re*ssq*dEvdT/Ts2
clha = rhoa*lath*usm*re*rd
_RL rd,re,rh,tta,toa
re = ren/(exf_one + ren*(ztln-psixh)/karman)
qstar(i,j) = re*delq(i,j)