From 15c9f06463f1db23afd383b653f26dd10255d50a Mon Sep 17 00:00:00 2001 From: Gustavo Marques Date: Wed, 2 Oct 2019 15:39:11 -0600 Subject: [PATCH] Hard-code min of SN to be 1.0e-7 --- src/parameterizations/lateral/MOM_MEKE.F90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/parameterizations/lateral/MOM_MEKE.F90 b/src/parameterizations/lateral/MOM_MEKE.F90 index bb825cdd2d..a64a73ae4c 100644 --- a/src/parameterizations/lateral/MOM_MEKE.F90 +++ b/src/parameterizations/lateral/MOM_MEKE.F90 @@ -775,7 +775,7 @@ subroutine MEKE_equilibrium(CS, MEKE, G, GV, US, SN_u, SN_v, drag_rate_visc, I_m ! Equation 1 of Jansen et al. (2015), balancing the GEOMETRIC GM coefficient against ! bottom drag (Equations 3 and 12) ! TODO: create a run time parameter for limitting SN. - MEKE%MEKE(i,j) = (CS%MEKE_GEOMETRIC_alpha * MIN(SN,1.e-5) * US%Z_to_m*G%bathyT(i,j))**2 / (CS%cdrag**2 * bottomFac2**3) + MEKE%MEKE(i,j) = (CS%MEKE_GEOMETRIC_alpha * MIN(SN,1.e-7) * US%Z_to_m*G%bathyT(i,j))**2 / (CS%cdrag**2 * bottomFac2**3) else MEKE%MEKE(i,j) = (US%Z_to_m*G%bathyT(i,j)*SN / (8*CS%cdrag))**2 endif