diff --git a/config/model_configs/prognostic_edmfx_bomex_column.yml b/config/model_configs/prognostic_edmfx_bomex_column.yml
index fa130cb568c..a020c69ee89 100644
--- a/config/model_configs/prognostic_edmfx_bomex_column.yml
+++ b/config/model_configs/prognostic_edmfx_bomex_column.yml
@@ -14,6 +14,7 @@ edmfx_entr_model: "Generalized"
 edmfx_detr_model: "Generalized"
 edmfx_sgs_mass_flux: false
 edmfx_sgs_diffusive_flux: true
+vert_diff: "false"
 edmfx_nh_pressure: true
 edmfx_velocity_relaxation: true
 prognostic_tke: true
@@ -34,5 +35,5 @@ diagnostics:
     period: 10mins
   - short_name: [arup, waup, taup, thetaaup, haup, husup, hurup, clwup, cliup, waen, taen, thetaaen, haen, husen, huren, clwen, clien, tke]
     period: 10mins
-  - short_name: [entr, detr, lmix, bgrad, strain, edt, evu]
+  - short_name: [entr, detr, lmix, bgrad, strain]
     period: 10mins
diff --git a/post_processing/ci_plots.jl b/post_processing/ci_plots.jl
index f072752667e..7e476ce7a6f 100644
--- a/post_processing/ci_plots.jl
+++ b/post_processing/ci_plots.jl
@@ -1050,6 +1050,8 @@ function make_plots(
         "waup",
         "tke",
         "lmix",
+        "bgrad",
+        "strain",
         "arup",
         "hus",
         "husup",
diff --git a/src/cache/precomputed_quantities.jl b/src/cache/precomputed_quantities.jl
index c10e25ed743..2f84595108b 100644
--- a/src/cache/precomputed_quantities.jl
+++ b/src/cache/precomputed_quantities.jl
@@ -38,7 +38,7 @@ function precomputed_quantities(Y, atmos)
             !(atmos.turbconv_model isa DiagnosticEDMFX)
     @assert !(atmos.moisture_model isa DryModel) ||
             !(atmos.turbconv_model isa PrognosticEDMFX)
-    @assert isnothing(atmos.turbconv_model) || isnothing(atmos.vert_diff)
+    # @assert isnothing(atmos.turbconv_model) || isnothing(atmos.vert_diff)
     TST = thermo_state_type(atmos.moisture_model, FT)
     SCT = SurfaceConditions.surface_conditions_type(atmos, FT)
     cspace = axes(Y.c)
diff --git a/src/cache/prognostic_edmf_precomputed_quantities.jl b/src/cache/prognostic_edmf_precomputed_quantities.jl
index a00c7fa4d81..f59521211d1 100644
--- a/src/cache/prognostic_edmf_precomputed_quantities.jl
+++ b/src/cache/prognostic_edmf_precomputed_quantities.jl
@@ -304,6 +304,14 @@ function set_prognostic_edmf_precomputed_quantities_closures!(Y, p, t)
     turbconv_params = CAP.turbconv_params(params)
     c_m = CAP.tke_ed_coeff(turbconv_params)
     @. ᶜK_u = c_m * ᶜmixing_length * sqrt(max(ᶜtke⁰, 0))
+    # interior_uₕ = Fields.level(Y.c.uₕ, 1)
+    # ᶜΔz_surface = Fields.Δz_field(interior_uₕ)
+    # @. ᶜK_u = eddy_diffusivity_coefficient(
+    #     params.C_E,
+    #     norm(interior_uₕ),
+    #     ᶜΔz_surface / 2,
+    #     ᶜp,
+    # )
     @. ᶜK_h = ᶜK_u / ᶜprandtl_nvec
 
     ρatke_flux_values = Fields.field_values(ρatke_flux)
diff --git a/src/initial_conditions/initial_conditions.jl b/src/initial_conditions/initial_conditions.jl
index fd1f6caa342..235ad61a2fc 100644
--- a/src/initial_conditions/initial_conditions.jl
+++ b/src/initial_conditions/initial_conditions.jl
@@ -811,7 +811,7 @@ for IC in (:Soares, :Bomex, :LifeCycleTan2018, :ARM_SGP)
                 ),
                 velocity = Geometry.UVector(u(z)),
                 turbconv_state = EDMFState(;
-                    tke = prognostic_tke ? FT(0) : tke(z),
+                    tke = prognostic_tke ? FT(0) : FT(0) #tke(z),
                 ),
             )
         end
diff --git a/src/prognostic_equations/gm_sgs_closures.jl b/src/prognostic_equations/gm_sgs_closures.jl
index 07a752e7fb2..888fd429b18 100644
--- a/src/prognostic_equations/gm_sgs_closures.jl
+++ b/src/prognostic_equations/gm_sgs_closures.jl
@@ -18,7 +18,7 @@ function smagorinsky_lilly_length(c_smag, N_eff, dz, Pr, ϵ_st)
     return N_eff > FT(0) ?
            c_smag *
            dz *
-           max(0, 1 - N_eff^2 / Pr / 2 / max(ϵ_st, eps(FT)))^(1 / 4) :
+           max(0, 1 - N_eff^2 / Pr / 2 / max(ϵ_st, FT(1e-5)))^(1 / 4) :
            c_smag * dz
 end
 
diff --git a/toml/prognostic_edmfx_bomex.toml b/toml/prognostic_edmfx_bomex.toml
index 5be8788fc4e..99e3e150722 100644
--- a/toml/prognostic_edmfx_bomex.toml
+++ b/toml/prognostic_edmfx_bomex.toml
@@ -1,3 +1,6 @@
+[C_E]
+value = 0.0176
+
 [EDMF_surface_area]
 value = 0.1