wip

CliMA · Dec 12, 2023 · 988ff6a · 988ff6a
1 parent ddb106a
commit 988ff6a
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/cache/diagnostic_edmf_precomputed_quantities.jl b/src/cache/diagnostic_edmf_precomputed_quantities.jl
@@ -369,7 +369,11 @@ function set_diagnostic_edmf_precomputed_quantities_do_integral!(Y, p, t)
                 p.atmos.edmfx_entr_model,
             )
 
-            if i > 2
+            # We don't have an upper limit to entrainment for the first level 
+            # (calculated at i=2), as the vertical at the first level is zero
+            if i == 2
+                @. entrʲ_prev_level = max(entrʲ_prev_level, 0)
+            else
                 @. entrʲ_prev_level = limit_entrainment(
                     entrʲ_prev_level,
                     draft_area(ρaʲ_prev_level, ρʲ_prev_level),