Modifying samfdeepcnv.f, samfdeepcnv.meta, samfshalcnv.f, samfshalcnv.meta so that convective adjustment time can be tunable, introducing two parameters cat_adj_deep and cat_adj_shal for this

Author: JunghoonShin-NOAA

physics/CONV/SAMF/samfdeepcnv.f

@@ -85,7 +85,8 @@ subroutine samfdeepcnv_run (im,km,first_time_step,restart,        &
      &    clam,c0s,c1,betal,betas,evef,pgcon,asolfac,cscale,            &
      &    do_ca, ca_closure, ca_entr, ca_trigger, nthresh,ca_deep,      &
      &    rainevap,sigmain,sigmaout,omegain,omegaout,betadcu,betamcu,   &
-     &    betascu,maxMF,do_mynnedmf,sigmab_coldstart,errmsg,errflg)
+     &    betascu,maxMF,do_mynnedmf,sigmab_coldstart,cat_adj_deep,      &
+     &    errmsg,errflg)
 
 !
       use machine , only : kind_phys
@@ -139,6 +140,9 @@ subroutine samfdeepcnv_run (im,km,first_time_step,restart,        &
      &                     evef,  pgcon
       character(len=*), intent(out) :: errmsg
       integer,          intent(out) :: errflg
+
+      ! for HAFS
+      real(kind_phys),           intent(in) :: cat_adj_deep
 !
 !------local variables
       integer              i, indx, jmn, k, kk, km1, n
@@ -2952,7 +2956,8 @@ subroutine samfdeepcnv_run (im,km,first_time_step,restart,        &
            umean(i) = max(umean(i), 1.)
            tauadv = gdx(i) / umean(i)
            advfac(i) = tauadv / dtconv(i)
-           advfac(i) = min(advfac(i), 1.)
+           advfac(i) = min(cat_adj_deep*advfac(i), 1.)
+           !advfac(i) = min(0.85*advfac(i), 1.)  !shin
         endif
       enddo
 

physics/CONV/SAMF/samfdeepcnv.meta

@@ -762,6 +762,14 @@
   type = real
   kind = kind_phys
   intent = out
+[cat_adj_deep]
+  standard_name = Adjustment_for_convective_advection_time_for_deep
+  long_name = Adjustment for convective advection time for deep
+  units = none
+  dimensions = ()
+  type = real
+  kind = kind_phys
+  intent = in
 [errmsg]
   standard_name = ccpp_error_message
   long_name = error message for error handling in CCPP

physics/CONV/SAMF/samfshalcnv.f

@@ -59,7 +59,7 @@ subroutine samfshalcnv_run(im,km,itc,ntc,cliq,cp,cvap,            &
      &     dot,ncloud,hpbl,ud_mf,dt_mf,cnvw,cnvc,                       &
      &     clam,c0s,c1,evef,pgcon,asolfac,hwrf_samfshal,                & 
      &     sigmain,sigmaout,omegain,omegaout,betadcu,betamcu,betascu,   &
-     &     errmsg,errflg)
+     &     cat_adj_shal,errmsg,errflg)
 !
       use machine , only : kind_phys
       use funcphys , only : fpvs
@@ -100,6 +100,9 @@ subroutine samfshalcnv_run(im,km,itc,ntc,cliq,cp,cvap,            &
      &     restart,progsigma,progomega
       character(len=*), intent(out) :: errmsg
       integer,          intent(out) :: errflg
+
+      ! for HAFS
+      real(kind_phys),           intent(in) :: cat_adj_shal
 !
 !  local variables
       integer              i,j,indx, k, kk, km1, n
@@ -1992,7 +1995,8 @@ subroutine samfshalcnv_run(im,km,itc,ntc,cliq,cp,cvap,            &
            umean(i) = max(umean(i), 1.)
            tauadv = gdx(i) / umean(i)
            advfac(i) = tauadv / dtconv(i)
-           advfac(i) = min(advfac(i), 1.)
+           advfac(i) = min(cat_adj_shal*advfac(i), 1.)
+           !advfac(i) = min(0.85*advfac(i), 1.)  !shin
         endif
       enddo
 c

physics/CONV/SAMF/samfshalcnv.meta

@@ -550,6 +550,14 @@
   dimensions = ()
   type = real
   intent = in
+[cat_adj_shal]
+  standard_name = Adjustment_for_convective_advection_time_for_shallow
+  long_name = Adjustment for convective advection time for shallow
+  units = none
+  dimensions = ()
+  type = real
+  kind = kind_phys
+  intent = in
 [errmsg]
   standard_name = ccpp_error_message
   long_name = error message for error handling in CCPP