only limit w when buoyancy is negative

CliMA · Dec 22, 2023 · 2d55f73 · 2d55f73
1 parent c4c46ea
commit 2d55f73
Showing 1 changed file with 0 additions and 8 deletions.
diff --git a/src/cache/diagnostic_edmf_precomputed_quantities.jl b/src/cache/diagnostic_edmf_precomputed_quantities.jl
@@ -636,43 +636,35 @@ function set_diagnostic_edmf_precomputed_quantities_do_integral!(Y, p, t)
                 @. ρaʲ_level = ifelse(
                     (
                         u³ʲ_data_prev_halflevel * u³ʲ_data_prev_halflevel <
-                        abs(
                             ∇Φ³_data_prev_level *
                             (ρʲ_prev_level - ρ_prev_level) / ρʲ_prev_level,
-                        )
                     ),
                     0,
                     ρaʲ_level,
                 )
                 @. u³ʲ_halflevel = ifelse(
                     (
                         u³ʲ_data_prev_halflevel * u³ʲ_data_prev_halflevel <
-                        abs(
                             ∇Φ³_data_prev_level *
                             (ρʲ_prev_level - ρ_prev_level) / ρʲ_prev_level,
-                        )
                     ),
                     u³_halflevel,
                     u³ʲ_halflevel,
                 )
                 @. h_totʲ_level = ifelse(
                     (
                         u³ʲ_data_prev_halflevel * u³ʲ_data_prev_halflevel <
-                        abs(
                             ∇Φ³_data_prev_level *
                             (ρʲ_prev_level - ρ_prev_level) / ρʲ_prev_level,
-                        )
                     ),
                     h_tot_level,
                     h_totʲ_level,
                 )
                 @. q_totʲ_level = ifelse(
                     (
                         u³ʲ_data_prev_halflevel * u³ʲ_data_prev_halflevel <
-                        abs(
                             ∇Φ³_data_prev_level *
                             (ρʲ_prev_level - ρ_prev_level) / ρʲ_prev_level,
-                        )
                     ),
                     q_tot_level,
                     q_totʲ_level,