Skip to content

Commit

Permalink
wip
Browse files Browse the repository at this point in the history
  • Loading branch information
szy21 committed Aug 30, 2023
1 parent 29cfde9 commit 40f9b7f
Show file tree
Hide file tree
Showing 2 changed files with 3 additions and 3 deletions.
2 changes: 1 addition & 1 deletion src/cache/diagnostic_edmf_precomputed_quantities.jl
Original file line number Diff line number Diff line change
Expand Up @@ -205,7 +205,7 @@ function set_diagnostic_edmf_precomputed_quantities!(Y, p, t)
ρʲ_int_level = Fields.field_values(Fields.level(ᶜρʲ, 1))

@. u³ʲ_int_halflevel = CT3(
Geometry.WVector($(FT(0.11)), local_geometry_int_halflevel),
Geometry.WVector($(FT(0.1)), local_geometry_int_halflevel),
local_geometry_int_halflevel,
)
@. h_totʲ_int_level = sgs_scalar_first_interior_bc(
Expand Down
4 changes: 2 additions & 2 deletions src/prognostic_equations/edmfx_closures.jl
Original file line number Diff line number Diff line change
Expand Up @@ -184,9 +184,9 @@ function pi_groups_entr_detr(
# TODO - Temporary: Switch to Π groups after simple tests are done
# (kinematic, bubble, Bomex)
# and/or we can calibrate things in ClimaAtmos
entr = max(0, min(FT(0.12) * (ᶜbuoyʲ - ᶜbuoy⁰) / ᶜwʲ, 1 / dt))
#entr = max(0, min(FT(0.12) * (ᶜbuoyʲ - ᶜbuoy⁰) / ᶜwʲ, 1 / dt))
#detr = max(0, min(detr_coeff * ᶜwʲ, 1 / dt))
detr = max(0, min(FT(0.12) * (ᶜbuoyʲ - ᶜbuoy⁰) / ᶜwʲ, 1 / dt))
#detr = max(0, min(FT(0.12) * (ᶜbuoyʲ - ᶜbuoy⁰) / ᶜwʲ, 1 / dt))

return (; entr, detr)
end
Expand Down

0 comments on commit 40f9b7f

Please sign in to comment.