Merge pull request #879 from mdtoyNOAA/main_new_cdmbgwd

Added capability for cdmbgwd(1) to scale GSL blocking drag
NCAR · Mar 30, 2022 · a8e0f0c · a8e0f0c
2 parents 85f21a4 + 34d0bbe
commit a8e0f0c
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/physics/drag_suite.F90 b/physics/drag_suite.F90
@@ -531,7 +531,8 @@ subroutine drag_suite_run(                                           &
 ! non-dim sub grid mtn drag Amp (*j*)
 !     cdmb = 1.0/float(IMX/192)
 !     cdmb = 192.0/float(IMX)
-      cdmb = 4.0 * 192.0/float(IMX)
+      ! New cdmbgwd addition for GSL blocking drag
+      cdmb = 1.0
       if (cdmbgwd(1) >= 0.0) cdmb = cdmb * cdmbgwd(1)
 
 !>-# Orographic Gravity Wave Drag Section
@@ -1237,7 +1238,8 @@ subroutine drag_suite_run(                                           &
 !--------- compute flow-blocking stress
 !
                cd = max(2.0-1.0/od(i),0.0)
-               taufb(i,kts) = 0.5 * roll(i) * coefm(i) /                   &
+               ! New cdmbgwd addition for GSL blocking drag
+               taufb(i,kts) = cdmb * 0.5 * roll(i) * coefm(i) /            &
                                  max(dxmax_ls,dxy(i))**2 * cd * dxyp(i) *  &
                                  olp(i) * zblk * ulow(i)**2
                tautem = taufb(i,kts)/float(kblk-kts)