56 real :: akm(nwzmax),bkm(nwzmax)
57 real :: ps,tsurf,tdsurf,tlev,ustar,hf,e,
ew,tv,rhoa,plev
58 real :: ak1,bk1,theta,thetastar
62 tv=tsurf*(1.+0.378*e/ps)
64 ak1=(akm(1)+akm(2))/2.
65 bk1=(bkm(1)+bkm(2))/2.
67 theta=tlev*(100000./plev)**(r_air/cpa)
68 if (ustar.le.0.) ustar=1.e-8
69 thetastar=hf/(rhoa*cpa*ustar)
70 if(abs(thetastar).gt.1.e-10)
then
real function obukhov_ecmwf(ps, tsurf, tdsurf, tlev, ustar, hf, akm, bkm)