diff --git a/physics/drag_suite.F90 b/physics/drag_suite.F90 index 7fea98b13..b4bd4e4d9 100644 --- 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)