wip

CliMA · Oct 10, 2023 · 503265d · 503265d
1 parent 2599c62
commit 503265d
Show file tree

Hide file tree

Showing 2 changed files with 24 additions and 0 deletions.
diff --git a/config/model_configs/advective_edmfx_adv_test_box.yml b/config/model_configs/advective_edmfx_adv_test_box.yml
@@ -0,0 +1,23 @@
+job_id: "advective_edmfx_adv_test_box"
+initial_condition: "MoistAdiabaticProfileEDMFX"
+advection_test: true
+turbconv: "advective_edmfx"
+ode_algo: "SSP33ShuOsher"
+edmfx_upwinding: "first_order"
+config: "box"
+moist: "equil"
+hyperdiff: "true"
+kappa_4: 1e8
+x_max: 1e4
+y_max: 1e4
+z_max: 5.5e4
+x_elem: 2
+y_elem: 2
+z_elem: 63
+dz_bottom: 30.0
+dz_top: 3000.0
+dt: "10secs"
+t_end: "3600secs"
+dt_save_to_disk: "100secs"
+FLOAT_TYPE: "Float64"
+toml: [toml/edmfx_box_advection.toml]
diff --git a/src/prognostic_equations/edmfx_entr_detr.jl b/src/prognostic_equations/edmfx_entr_detr.jl
@@ -303,6 +303,7 @@ function detrainment(
     detr_tau = CAP.detr_tau(params)
     detr = min(1 / detr_tau, 1 / dt)
     return detr
+end
 
 edmfx_entr_detr_tendency!(Yₜ, Y, p, t, colidx, turbconv_model) = nothing
 function edmfx_entr_detr_tendency!(Yₜ, Y, p, t, colidx, turbconv_model::EDMFX)