diff --git a/phys/module_bl_keps.F b/phys/module_bl_keps.F index 8c5be78733..c804d65321 100644 --- a/phys/module_bl_keps.F +++ b/phys/module_bl_keps.F @@ -1073,16 +1073,16 @@ subroutine surface_bl_pra_ri(kms,kme,kts,kte,dz,z,rho,g,cp,z0,sflux,raten,pi1d,p zz=(zl1+z0)/zl1 zol1 = max(b_ric*psim*psim/psih,rimin) if(sfcflg)then - zol1 = min(zol1,-zfmin) + zol1 = zz*min(zol1,-zfmin) else - zol1 = max(zol1,zfmin) + zol1 = zz*max(zol1,zfmin) endif hol1 = zol1*pblh/zl1*sfcfrac if(sfcflg)then - phim = (1.-aphi16*zol1*zz)**(-1./4.) - phih = (1.-aphi16*zol1*zz)**(-1./2.) - phieps=1.-zol1*zz + phim = (1.-aphi16*zol1)**(-1./4.) + phih = (1.-aphi16*zol1)**(-1./2.) + phieps=1.-zol1 phim_sl = (1.-aphi16*hol1)**(-1./4.) phih_sl = (1.-aphi16*hol1)**(-1./2.) bfx0 = max(sflux,0.) @@ -1093,9 +1093,9 @@ subroutine surface_bl_pra_ri(kms,kme,kts,kte,dz,z,rho,g,cp,z0,sflux,raten,pi1d,p else phim_sl = (1.+aphi5*hol1) phih_sl= phim_sl - phim=(1.+aphi5*zol1*zz) - phieps=(1+2.5*(zol1*zz)**0.6)**(3./2.) - phih= phim_sl + phim=(1.+aphi5*zol1) + phieps=(1+2.5*(zol1)**0.6)**(3./2.) + phih= phim wstar3=0. wstar3_2=0. endif diff --git a/phys/module_pbl_driver.F b/phys/module_pbl_driver.F index f703071765..569c3cf0dd 100644 --- a/phys/module_pbl_driver.F +++ b/phys/module_pbl_driver.F @@ -1899,7 +1899,7 @@ SUBROUTINE pbl_driver( & ,B_E_BEP=b_e_bep & ,SF_BEP=sf_bep,VL_BEP=vl_bep & ,BR=br,ZNT=znt & - ,PSIM=psim,PSIH=psih & + ,PSIM=fm,PSIH=fhh & ,IDS=ids,IDE=ide,JDS=jds,JDE=jde,KDS=kds,KDE=kde & ,IMS=ims,IME=ime,JMS=jms,JME=jme,KMS=kms,KME=kme & ,ITS=its,ITE=ite,JTS=jts,JTE=jte,KTS=kts,KTE=kte )