Merge pull request #3644 from CliMA/as/fix-dyamondic

Bugfix: DYAMOND ic with EDOnly

Author: szy21

config/gpu_configs/gpu_aquaplanet_dyamond_summer.yml

@@ -12,7 +12,8 @@ rad: "allskywithclear"
 insolation: "timevarying"
 dt_rad: "1hours"
 dt_cloud_fraction: "1hours"
-vert_diff: "DecayWithHeightDiffusion"
+turbconv: "edonly_edmfx"
+edmfx_sgs_diffusive_flux: true
 implicit_diffusion: true
 approximate_linear_solve_iters: 2
 surface_setup: "DefaultMoninObukhov"

src/initial_conditions/initial_conditions.jl

@@ -487,7 +487,7 @@ function overwrite_initial_conditions!(
 
     if hasproperty(Y.c, :sgs⁰) && hasproperty(Y.c.sgs⁰, :ρatke)
         # NOTE: This is not the most consistent, but it is better than NaNs
-        fill!(Y.c.ρatke, 0)
+        fill!(Y.c.sgs⁰.ρatke, 0)
     end
 
     return nothing