Zs/implicit tke advection

config/model_configs/prognostic_edmfx_bomex_column.yml

@@ -0,0 +1,35 @@
+job_id: "prognostic_edmfx_bomex_column"
+initial_condition: "Bomex"
+subsidence: "Bomex"
+edmf_coriolis: "Bomex"
+ls_adv: "Bomex"
+surface_setup: "Bomex"
+turbconv: "prognostic_edmfx"
+implicit_diffusion: true
+approximate_linear_solve_iters: 2
+max_newton_iters_ode: 3
+edmfx_upwinding: first_order
+edmfx_entr_model: "Generalized"
+edmfx_detr_model: "Generalized"
+edmfx_sgs_mass_flux: false
+edmfx_sgs_diffusive_flux: true
+edmfx_nh_pressure: true
+edmfx_velocity_relaxation: true
+prognostic_tke: true
+moist: "equil"
+config: "column"
+z_max: 3e3
+z_elem: 60
+z_stretch: false
+perturb_initstate: false
+dt: "60secs"
+t_end: "6hours"
+dt_save_to_disk: "10mins"
+toml: [toml/prognostic_edmfx_bomex_box.toml]
+netcdf_output_at_levels: true
+netcdf_interpolation_num_points: [2, 2, 60]
+diagnostics:
+  - short_name: [ts, ta, thetaa, ha, pfull, rhoa, ua, va, wa, hur, hus, cl, clw, cli, hussfc, evspsbl]
+    period: 10mins
+  - short_name: [arup, waup, taup, rhoaup, thetaaup, haup, husup, hurup, clwup, cliup, waen, taen, thetaaen, haen, husen, huren, clwen, clien, tke, lmix]
+    period: 10mins

post_processing/ci_plots.jl

@@ -943,6 +943,7 @@ EDMFBoxPlots = Union{
     Val{:prognostic_edmfx_gabls_box},
     Val{:prognostic_edmfx_bomex_fixtke_box},
     Val{:prognostic_edmfx_bomex_box},
+    Val{:prognostic_edmfx_bomex_column},
     Val{:prognostic_edmfx_bomex_stretched_box},
     Val{:prognostic_edmfx_dycoms_rf01_box},
     Val{:prognostic_edmfx_rico_column},
@@ -1037,6 +1038,7 @@ function make_plots(::EDMFBoxPlots, output_paths::Vector{<:AbstractString})
         "haup",
         "waup",
         "tke",
+        "lmix",
         "arup",
         "hus",
         "husup",

src/prognostic_equations/edmfx_tke.jl

@@ -55,7 +55,7 @@ function edmfx_tke_tendency!(
                 )
         end
         # pressure work
-        @. Yₜ.c.sgs⁰.ρatke[colidx] += ᶜtke_press[colidx]
+        #@. Yₜ.c.sgs⁰.ρatke[colidx] += ᶜtke_press[colidx]
     end
 
     return nothing