Merge pull request #303 from hafs-community/feature/3dtke_update

Update scale-aware 3DTKE EDMF and introduce TTE-EDMF in GFS PBL scheme

Author: rhaesung

physics/PBL/SATMEDMF/mfscuq.f

@@ -15,9 +15,7 @@ subroutine mfscuq(im,ix,km,kmscu,ntcw,ntrac1,delt,
      &   cnvflg,zl,zm,q1,t1,u1,v1,plyr,pix,
      &   thlx,thvx,thlvx,gdx,thetae,
      &   krad,mrad,radmin,buo,wush,tkemean,vez0fun,xmfd,
-     &   tcdo,qcdo,ucdo,vcdo,xlamdeq,a1,
-!The following flag is for SA-3D-TKE (kyf)
-     &   sa3dtke)
+     &   tcdo,qcdo,ucdo,vcdo,xlamdeq,a1)
 !
       use machine , only : kind_phys
       use funcphys , only : fpvs
@@ -33,7 +31,6 @@ subroutine mfscuq(im,ix,km,kmscu,ntcw,ntrac1,delt,
       integer   krad(im), mrad(im)
 !
       logical cnvflg(im)
-      logical sa3dtke !flag for SA-3D-TKE scheme (kyf)
       real(kind=kind_phys) delt
       real(kind=kind_phys) q1(ix,km,ntrac1),t1(ix,km),
      &                     u1(ix,km),      v1(ix,km),
@@ -430,16 +427,6 @@ subroutine mfscuq(im,ix,km,kmscu,ntcw,ntrac1,delt,
         endif
       enddo
 !
-!> - Set updraft fraction to 0 when using SA-3D-TKE scheme (kyf)
-!! Scale-aware capability is done with pfnl in satmedmfvdifq.F
-!! Zhu et al. (2025)
-!
-      if (sa3dtke) then
-        do i = 1, im
-          sigma(i) = 0.
-        enddo
-      endif
-!
 !> - Compute scale-aware function based on 
 !! Arakawa and Wu (2013) \cite arakawa_and_wu_2013
 !
@@ -460,6 +447,9 @@ subroutine mfscuq(im,ix,km,kmscu,ntcw,ntrac1,delt,
         do i = 1, im
           if(cnvflg(i) .and.
      &       (k >= mrad(i) .and. k < krad(i))) then
+             if (sigma(i) > ra1(i)) then
+               xmfd(i,k) = sigma(i) * xmfd(i,k) / ra1(i)
+             endif
              xmfd(i,k) = scaldfunc(i) * xmfd(i,k)
              dz   = zl(i,k+1) - zl(i,k)
              xmmx = dz / dt2