From 60088d3b8985c7b4f2ef23edae4bd4e665c6cfde Mon Sep 17 00:00:00 2001 From: Zhaoyi Shen <11598433+szy21@users.noreply.github.com> Date: Wed, 16 Aug 2023 22:49:12 -0700 Subject: [PATCH] min mixing length 10x --- src/prognostic_equations/edmfx_tke.jl | 3 +-- 1 file changed, 1 insertion(+), 2 deletions(-) diff --git a/src/prognostic_equations/edmfx_tke.jl b/src/prognostic_equations/edmfx_tke.jl index ef8f26b9b2..238b614bbf 100644 --- a/src/prognostic_equations/edmfx_tke.jl +++ b/src/prognostic_equations/edmfx_tke.jl @@ -53,8 +53,7 @@ function edmfx_tke_tendency!( end # pressure work # dissipation - #min_mixing_length = @. c_d * dt * max(ᶜtke⁰[colidx], eps(FT))^(FT(1) / 2) - min_mixing_length = FT(100) + min_mixing_length = @. FT(10) * c_d * dt * max(ᶜtke⁰[colidx], eps(FT))^(FT(1) / 2) @. Yₜ.c.sgs⁰.ρatke[colidx] -= Y.c.ρ[colidx] * c_d * max(ᶜtke⁰[colidx], 0)^(FT(3) / 2) / max(ᶜmixing_length[colidx], min_mixing_length)