min mixing length 10x

CliMA · Aug 17, 2023 · 60088d3 · 60088d3
1 parent b8345ea
commit 60088d3
Showing 1 changed file with 1 addition and 2 deletions.
diff --git 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)