From 2d55f7336ff204c8cc56d3f5862c73733e3470be Mon Sep 17 00:00:00 2001
From: Zhaoyi Shen <11598433+szy21@users.noreply.github.com>
Date: Fri, 22 Dec 2023 13:07:18 -0800
Subject: [PATCH] only limit w when buoyancy is negative

---
 src/cache/diagnostic_edmf_precomputed_quantities.jl | 8 --------
 1 file changed, 8 deletions(-)

diff --git a/src/cache/diagnostic_edmf_precomputed_quantities.jl b/src/cache/diagnostic_edmf_precomputed_quantities.jl
index 050fb0e90db..667e4e2e0b1 100644
--- a/src/cache/diagnostic_edmf_precomputed_quantities.jl
+++ b/src/cache/diagnostic_edmf_precomputed_quantities.jl
@@ -636,10 +636,8 @@ 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,
@@ -647,10 +645,8 @@ function set_diagnostic_edmf_precomputed_quantities_do_integral!(Y, p, t)
                 @. 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,
@@ -658,10 +654,8 @@ function set_diagnostic_edmf_precomputed_quantities_do_integral!(Y, p, t)
                 @. 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,
@@ -669,10 +663,8 @@ function set_diagnostic_edmf_precomputed_quantities_do_integral!(Y, p, t)
                 @. 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,