Skip to content

Commit

Permalink
fix
Browse files Browse the repository at this point in the history
  • Loading branch information
szy21 committed Dec 22, 2023
1 parent 2d55f73 commit d9df3b3
Showing 1 changed file with 4 additions and 8 deletions.
12 changes: 4 additions & 8 deletions src/cache/diagnostic_edmf_precomputed_quantities.jl
Original file line number Diff line number Diff line change
Expand Up @@ -636,35 +636,31 @@ function set_diagnostic_edmf_precomputed_quantities_do_integral!(Y, p, t)
@. ρaʲ_level = ifelse(
(
u³ʲ_data_prev_halflevel * u³ʲ_data_prev_halflevel <
∇Φ³_data_prev_level *
(ρʲ_prev_level - ρ_prev_level) / ρʲ_prev_level,
∇Φ³_data_prev_level * (ρʲ_prev_level - ρ_prev_level) / ρʲ_prev_level
),
0,
ρaʲ_level,
)
@. u³ʲ_halflevel = ifelse(
(
u³ʲ_data_prev_halflevel * u³ʲ_data_prev_halflevel <
∇Φ³_data_prev_level *
(ρʲ_prev_level - ρ_prev_level) / ρʲ_prev_level,
∇Φ³_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 <
∇Φ³_data_prev_level *
(ρʲ_prev_level - ρ_prev_level) / ρʲ_prev_level,
∇Φ³_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 <
∇Φ³_data_prev_level *
(ρʲ_prev_level - ρ_prev_level) / ρʲ_prev_level,
∇Φ³_data_prev_level * (ρʲ_prev_level - ρ_prev_level) / ρʲ_prev_level
),
q_tot_level,
q_totʲ_level,
Expand Down

0 comments on commit d9df3b3

Please sign in to comment.