fix non stochastic ePBL calculation

jiandewang · Jan 7, 2021 · 6bb9d0b · 6bb9d0b
1 parent 1d7ffa3
commit 6bb9d0b
Showing 1 changed file with 5 additions and 2 deletions.
diff --git a/src/parameterizations/vertical/MOM_energetic_PBL.F90 b/src/parameterizations/vertical/MOM_energetic_PBL.F90
@@ -971,8 +971,11 @@ subroutine ePBL_column(h, u, v, T0, S0, dSV_dT, dSV_dS, TKE_forcing, B_flux, abs
         if (Idecay_len_TKE > 0.0) exp_kh = exp(-h(k-1)*Idecay_len_TKE)
         if (CS%TKE_diagnostics) &
           eCD%dTKE_mech_decay = eCD%dTKE_mech_decay + (exp_kh-1.0) * mech_TKE * I_dtdiag
-        mech_TKE = mech_TKE * exp_kh
-        if (stoch_epbl) mech_TKE = mech_TKE * (1+(exp_kh-1) * t_rp2)
+        if (stoch_epbl) then ! perturb the TKE destruction
+           mech_TKE = mech_TKE * (1+(exp_kh-1) * t_rp2)
+        else
+           mech_TKE = mech_TKE * exp_kh
+        endif
 
         !   Accumulate any convectively released potential energy to contribute
         ! to wstar and to drive penetrating convection.