58 real :: ps,tsurf,tdsurf,tlev,ustar,hf,e,
ew,tv,rhoa,plev
59 real :: theta,thetastar
63 tv=tsurf*(1.+0.378*e/ps)
65 theta=tlev*(100000./plev)**(r_air/cpa)
66 if (ustar.le.0.) ustar=1.e-8
67 thetastar=hf/(rhoa*cpa*ustar)
68 if(abs(thetastar).gt.1.e-10)
then
real function obukhov_gfs(ps, tsurf, tdsurf, tlev, ustar, hf, plev)