From a31d08bdd46b7d88efe52d3523c85db1983b2021 Mon Sep 17 00:00:00 2001 From: Michael Toy Date: Fri, 11 Mar 2022 20:42:20 +0000 Subject: [PATCH] Added capability for cdmbgwd(1) to scale GSL blocking drag --- physics/drag_suite.F90 | 6 ++++-- 1 file changed, 4 insertions(+), 2 deletions(-) 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)