correction of stability functions for the tke-epsilon-tpe for consist…

…ency with other pbl schemes
wrf-model · Oct 16, 2024 · e5a930c · e5a930c
1 parent 0a11865
commit e5a930c
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 9 deletions.
diff --git 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
@@ -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  )