Merge pull request #3670 from CliMA/zs/exp_sfc_flux

szy21 · web-flow · commit f2686f5806ee · 2025-03-05T07:10:15.000Z
Make surface flux fully explicit
diff --git a/NEWS.md b/NEWS.md
@@ -1,11 +1,17 @@
 ClimaAtmos.jl Release Notes
 ============================
 
-v0.28.6
+main
 -------
 
 ### Features
 
+### Add a flag for disabling surface flux tendency
+Surface flux tendency is not controlled by `vert_diff` or `edmfx_sgs_diffusive_flux` anymore.
+Instead, it is controlled by the new flag `disable_surface_flux_tendency`. When it is set to 
+true, no surface flux tendency is applied, no matter what `surface_setup` is.
+This flag is set to false by default. PR [3670](https://github.com/CliMA/ClimaAtmos.jl/pull/3670).
+
 ### Automatically determine diagnostic resolution based on model resolution
 
 If `netcdf_interpolation_num_points` is not provided, `ClimaAtmos` will
diff --git a/config/default_configs/default_config.yml b/config/default_configs/default_config.yml
@@ -175,6 +175,9 @@ divergence_damping_factor:
 rayleigh_sponge:
   help: "Rayleigh sponge [`true`, `false` (default)]"
   value: false
+disable_surface_flux_tendency:
+  help: "(Bool) Whether to disable surface flux tendencies of momentum, energy, and tracers [`true`, `false` (default)]. When this flag is true, the surface flux tendency is not applied, no matter how surface conditions are computed."
+  value: false
 surface_setup:
   help: "Surface flux scheme [`DefaultExchangeCoefficients` (default), `DefaultMoninObukhov`]"
   value: "DefaultExchangeCoefficients"
diff --git a/config/gpu_configs/baroclinic_wave_helem30.yml b/config/gpu_configs/baroclinic_wave_helem30.yml
@@ -5,5 +5,6 @@ t_end: "1days"
 dt: "90secs"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
 initial_condition: "DryBaroclinicWave"
diff --git a/config/longrun_configs/longrun_dry_baroclinic_wave.yml b/config/longrun_configs/longrun_dry_baroclinic_wave.yml
@@ -5,6 +5,7 @@ dt: "90secs"
 t_end: "120days"
 initial_condition: "DryBaroclinicWave"
 dt_save_state_to_disk: "10days"
+disable_surface_flux_tendency: true
 diagnostics:
   - short_name: [pfull, wa, va, rv, ke]
     period: 10days
diff --git a/config/longrun_configs/longrun_dry_baroclinic_wave_he60.yml b/config/longrun_configs/longrun_dry_baroclinic_wave_he60.yml
@@ -5,6 +5,7 @@ dt: "30secs"
 t_end: "30days"
 initial_condition: "DryBaroclinicWave"
 dt_save_state_to_disk: "10days"
+disable_surface_flux_tendency: true
 netcdf_interpolation_num_points: [360, 180, 2]
 diagnostics:
   - short_name: [pfull, wa, va, rv, ke]
diff --git a/config/longrun_configs/longrun_dry_held_suarez.yml b/config/longrun_configs/longrun_dry_held_suarez.yml
@@ -8,4 +8,5 @@ forcing: "held_suarez"
 rayleigh_sponge: true
 viscous_sponge: true
 dt_save_state_to_disk: "10days"
+disable_surface_flux_tendency: true
 toml: [toml/longrun_held_suarez.toml]
diff --git a/config/longrun_configs/longrun_hydrostatic_balance.yml b/config/longrun_configs/longrun_hydrostatic_balance.yml
@@ -7,3 +7,4 @@ dt: "90secs"
 t_end: "360days"
 perturb_initstate: false
 discrete_hydrostatic_balance: true
+disable_surface_flux_tendency: true
diff --git a/config/longrun_configs/longrun_moist_baroclinic_wave.yml b/config/longrun_configs/longrun_moist_baroclinic_wave.yml
@@ -7,6 +7,7 @@ initial_condition: "MoistBaroclinicWave"
 moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "10days"
+disable_surface_flux_tendency: true
 toml: [toml/longrun_baroclinic_wave.toml]
 diagnostics:
   - short_name: [pfull, wa, va, ua, ta, rhoa, rv, hus, ke, clw, cli]
diff --git a/config/longrun_configs/longrun_moist_baroclinic_wave_1M.yml b/config/longrun_configs/longrun_moist_baroclinic_wave_1M.yml
@@ -10,6 +10,7 @@ initial_condition: "MoistBaroclinicWave"
 moist: "nonequil"
 precip_model: "1M"
 dt_save_state_to_disk: "15days"
+disable_surface_flux_tendency: true
 toml: [toml/longrun_baroclinic_wave_1M.toml]
 diagnostics:
   - short_name: [pfull, wa, va, ua, ta, rhoa, rv, hus, ke, clw, cli, husra, hussn]
diff --git a/config/longrun_configs/longrun_moist_baroclinic_wave_he60.yml b/config/longrun_configs/longrun_moist_baroclinic_wave_he60.yml
@@ -7,6 +7,7 @@ initial_condition: "MoistBaroclinicWave"
 moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "10days"
+disable_surface_flux_tendency: true
 netcdf_interpolation_num_points: [360, 180, 2]
 toml: [toml/longrun_baroclinic_wave.toml]
 diagnostics:
diff --git a/config/mbw_scaling_configs/moist_baroclinic_wave_helem_120_0M_ss.yml b/config/mbw_scaling_configs/moist_baroclinic_wave_helem_120_0M_ss.yml
@@ -9,4 +9,5 @@ moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
diff --git a/config/mbw_scaling_configs/moist_baroclinic_wave_helem_120_0M_ws.yml b/config/mbw_scaling_configs/moist_baroclinic_wave_helem_120_0M_ws.yml
@@ -9,4 +9,5 @@ moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
diff --git a/config/mbw_scaling_configs/moist_baroclinic_wave_helem_170_0M_ws.yml b/config/mbw_scaling_configs/moist_baroclinic_wave_helem_170_0M_ws.yml
@@ -9,4 +9,5 @@ moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
diff --git a/config/mbw_scaling_configs/moist_baroclinic_wave_helem_240_0M_ss.yml b/config/mbw_scaling_configs/moist_baroclinic_wave_helem_240_0M_ss.yml
@@ -9,4 +9,5 @@ moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
diff --git a/config/mbw_scaling_configs/moist_baroclinic_wave_helem_240_0M_ws.yml b/config/mbw_scaling_configs/moist_baroclinic_wave_helem_240_0M_ws.yml
@@ -9,4 +9,5 @@ moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
diff --git a/config/mbw_scaling_configs/moist_baroclinic_wave_helem_30_0M_ss.yml b/config/mbw_scaling_configs/moist_baroclinic_wave_helem_30_0M_ss.yml
@@ -9,4 +9,5 @@ moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
diff --git a/config/mbw_scaling_configs/moist_baroclinic_wave_helem_30_0M_ws.yml b/config/mbw_scaling_configs/moist_baroclinic_wave_helem_30_0M_ws.yml
@@ -9,4 +9,5 @@ moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
diff --git a/config/mbw_scaling_configs/moist_baroclinic_wave_helem_340_0M_ws.yml b/config/mbw_scaling_configs/moist_baroclinic_wave_helem_340_0M_ws.yml
@@ -9,4 +9,5 @@ moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
diff --git a/config/mbw_scaling_configs/moist_baroclinic_wave_helem_42_0M_ws.yml b/config/mbw_scaling_configs/moist_baroclinic_wave_helem_42_0M_ws.yml
@@ -9,4 +9,5 @@ moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
diff --git a/config/mbw_scaling_configs/moist_baroclinic_wave_helem_480_0M_ws.yml b/config/mbw_scaling_configs/moist_baroclinic_wave_helem_480_0M_ws.yml
@@ -9,4 +9,5 @@ moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
diff --git a/config/mbw_scaling_configs/moist_baroclinic_wave_helem_60_0M_ss.yml b/config/mbw_scaling_configs/moist_baroclinic_wave_helem_60_0M_ss.yml
@@ -9,4 +9,5 @@ moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
diff --git a/config/mbw_scaling_configs/moist_baroclinic_wave_helem_60_0M_ws.yml b/config/mbw_scaling_configs/moist_baroclinic_wave_helem_60_0M_ws.yml
@@ -9,4 +9,5 @@ moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
diff --git a/config/mbw_scaling_configs/moist_baroclinic_wave_helem_84_0M_ws.yml b/config/mbw_scaling_configs/moist_baroclinic_wave_helem_84_0M_ws.yml
@@ -9,4 +9,5 @@ moist: "equil"
 precip_model: "0M"
 dt_save_state_to_disk: "Inf"
 dt_save_to_sol: "Inf"
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
diff --git a/config/model_configs/aquaplanet_equil_clearsky_tvinsol_0M_slabocean.yml b/config/model_configs/aquaplanet_equil_clearsky_tvinsol_0M_slabocean.yml
@@ -11,4 +11,4 @@ rad: "clearsky"
 insolation: "timevarying"
 dt_rad: "1hours"
 prognostic_surface: "PrognosticSurfaceTemperature"
-#check_conservation: true
+check_conservation: true
diff --git a/config/model_configs/aquaplanet_equil_clearsky_tvinsol_0M_slabocean_ft64.yml b/config/model_configs/aquaplanet_equil_clearsky_tvinsol_0M_slabocean_ft64.yml
@@ -11,4 +11,4 @@ rad: "clearsky"
 insolation: "timevarying"
 dt_rad: "1hours"
 prognostic_surface: "PrognosticSurfaceTemperature"
-#check_conservation: true
+check_conservation: true
diff --git a/config/model_configs/baroclinic_wave.yml b/config/model_configs/baroclinic_wave.yml
@@ -3,6 +3,7 @@ initial_condition: "DryBaroclinicWave"
 dt: "400secs"
 t_end: "10days"
 deep_atmosphere: false
+disable_surface_flux_tendency: true
 diagnostics:
   - short_name: [pfull, ua, wa, va, rv, ta, ke]
     period: 1days
diff --git a/config/model_configs/baroclinic_wave_conservation.yml b/config/model_configs/baroclinic_wave_conservation.yml
@@ -2,6 +2,7 @@ initial_condition: "DryBaroclinicWave"
 dt: "400secs"
 t_end: "5days"
 dt_save_state_to_disk: "5days"
+disable_surface_flux_tendency: true
 check_conservation: true
 diagnostics:
   - short_name: [massa, energya]
diff --git a/config/model_configs/baroclinic_wave_equil.yml b/config/model_configs/baroclinic_wave_equil.yml
@@ -5,6 +5,7 @@ initial_condition: "MoistBaroclinicWave"
 dt: "450secs"
 t_end: "10days"
 moist: "equil"
+disable_surface_flux_tendency: true
 diagnostics:
   - short_name: [pfull, wa, va, rv, hus, ke]
     period: 1days
diff --git a/config/model_configs/baroclinic_wave_equil_conservation.yml b/config/model_configs/baroclinic_wave_equil_conservation.yml
@@ -2,6 +2,7 @@ initial_condition: "MoistBaroclinicWave"
 dt: "400secs"
 t_end: "5days"
 moist: "equil"
+disable_surface_flux_tendency: true
 dt_save_state_to_disk: "5days"
 check_conservation: true
 diagnostics:
diff --git a/config/model_configs/baroclinic_wave_equil_conservation_ft64.yml b/config/model_configs/baroclinic_wave_equil_conservation_ft64.yml
@@ -3,6 +3,7 @@ initial_condition: "MoistBaroclinicWave"
 dt: "400secs"
 t_end: "5days"
 moist: "equil"
+disable_surface_flux_tendency: true
 dt_save_state_to_disk: "5days"
 check_conservation: true
 diagnostics:
diff --git a/config/model_configs/baroclinic_wave_hughes2023.yml b/config/model_configs/baroclinic_wave_hughes2023.yml
@@ -3,6 +3,7 @@ initial_condition: "MoistBaroclinicWave"
 topography: "Hughes2023"
 topo_smoothing: false
 moist: equil
+disable_surface_flux_tendency: true
 dt: "200secs"
 t_end: "10days"
 diagnostics:
diff --git a/config/model_configs/baroclinic_wave_topography_dcmip_rs.yml b/config/model_configs/baroclinic_wave_topography_dcmip_rs.yml
@@ -6,6 +6,7 @@ initial_condition: "DryBaroclinicWave"
 t_end: "6days"
 dt: "200secs"
 hyperdiff: "ClimaHyperdiffusion"
+disable_surface_flux_tendency: true
 netcdf_output_at_levels: false
 diagnostics:
   - short_name: [pfull, wa, va, rv]
diff --git a/config/model_configs/box_cosine_hills_float64_test.yml b/config/model_configs/box_cosine_hills_float64_test.yml
@@ -12,6 +12,7 @@ dz_bottom: 10
 dt: "0.5secs" # CFL is slightly lower than for 2D case
 t_end: "30mins"
 rayleigh_sponge: true
+disable_surface_flux_tendency: true
 toml: [toml/steady_state_test.toml]
 check_steady_state: true
 output_default_diagnostics: false
diff --git a/config/model_configs/box_density_current_test.yml b/config/model_configs/box_density_current_test.yml
@@ -14,6 +14,7 @@ z_stretch: false
 dt: "0.25secs"
 t_end: "20secs"
 dt_save_state_to_disk: "5secs"
+disable_surface_flux_tendency: true
 netcdf_interpolation_num_points: [40, 40, 80]
 diagnostics:
   - short_name: [wa, ua, va, ta, thetaa, ha]
diff --git a/config/model_configs/box_hydrostatic_balance.yml b/config/model_configs/box_hydrostatic_balance.yml
@@ -3,6 +3,7 @@ config: "box"
 hyperdiff: false
 dt: "20secs"
 perturb_initstate: false
+disable_surface_flux_tendency: true
 diagnostics:
   - short_name: [wa, ua]
     period: 1days
diff --git a/config/model_configs/extruded_plane_cosine_hills_float64_test.yml b/config/model_configs/extruded_plane_cosine_hills_float64_test.yml
@@ -12,6 +12,7 @@ dz_bottom: 10
 dt: "0.5secs" # CFL is slightly lower than for 2D case
 t_end: "9hours"
 rayleigh_sponge: true
+disable_surface_flux_tendency: true
 toml: [toml/steady_state_test.toml]
 check_steady_state: true
 output_default_diagnostics: false
diff --git a/config/model_configs/hydrostatic_balance_ft64.yml b/config/model_configs/hydrostatic_balance_ft64.yml
@@ -4,4 +4,5 @@ t_end: "8days"
 discrete_hydrostatic_balance: true
 hyperdiff: false
 perturb_initstate: false
+disable_surface_flux_tendency: true
 FLOAT_TYPE: "Float64"
diff --git a/config/model_configs/plane_agnesi_mountain_float64_test.yml b/config/model_configs/plane_agnesi_mountain_float64_test.yml
@@ -10,6 +10,7 @@ dz_bottom: 10
 dt: "0.7secs"
 t_end: "2days"
 rayleigh_sponge: true
+disable_surface_flux_tendency: true
 toml: [toml/steady_state_test.toml]
 check_steady_state: true
 output_default_diagnostics: false
diff --git a/config/model_configs/plane_cosine_hills_float64_test.yml b/config/model_configs/plane_cosine_hills_float64_test.yml
@@ -10,6 +10,7 @@ dz_bottom: 10
 dt: "0.7secs"
 t_end: "2days"
 rayleigh_sponge: true
+disable_surface_flux_tendency: true
 toml: [toml/steady_state_test.toml]
 check_steady_state: true
 output_default_diagnostics: false
diff --git a/config/model_configs/plane_density_current_test.yml b/config/model_configs/plane_density_current_test.yml
@@ -9,6 +9,7 @@ z_stretch: false
 x_elem: 80
 config: "plane"
 z_max: 6400.0
+disable_surface_flux_tendency: true
 output_default_diagnostics: false
 diagnostics:
   - short_name: thetaa
diff --git a/config/model_configs/plane_no_topography_float64_test.yml b/config/model_configs/plane_no_topography_float64_test.yml
@@ -10,6 +10,7 @@ dz_bottom: 10
 dt: "0.7secs"
 t_end: "2days"
 rayleigh_sponge: true
+disable_surface_flux_tendency: true
 toml: [toml/steady_state_test.toml]
 check_steady_state: true
 output_default_diagnostics: false
diff --git a/config/model_configs/plane_schar_mountain_float32_test.yml b/config/model_configs/plane_schar_mountain_float32_test.yml
@@ -10,6 +10,7 @@ dz_bottom: 10
 dt: "0.7secs"
 t_end: "2days"
 rayleigh_sponge: true
+disable_surface_flux_tendency: true
 toml: [toml/steady_state_test.toml]
 check_steady_state: true
 output_default_diagnostics: false
diff --git a/config/model_configs/plane_schar_mountain_float64_test.yml b/config/model_configs/plane_schar_mountain_float64_test.yml
@@ -10,6 +10,7 @@ dz_bottom: 10
 dt: "0.7secs"
 t_end: "2days"
 rayleigh_sponge: true
+disable_surface_flux_tendency: true
 toml: [toml/steady_state_test.toml]
 check_steady_state: true
 output_default_diagnostics: false
diff --git a/config/model_configs/prognostic_edmfx_simpleplume_column.yml b/config/model_configs/prognostic_edmfx_simpleplume_column.yml
@@ -1,4 +1,5 @@
 initial_condition: "SimplePlume"
+disable_surface_flux_tendency: true
 surface_setup: "SimplePlume"
 turbconv: "prognostic_edmfx"
 implicit_diffusion: false
diff --git a/config/model_configs/single_column_hydrostatic_balance_ft64.yml b/config/model_configs/single_column_hydrostatic_balance_ft64.yml
@@ -5,6 +5,7 @@ dt_save_state_to_disk: "5days"
 initial_condition: "IsothermalProfile"
 config: "column"
 hyperdiff: false
+disable_surface_flux_tendency: true
 dt: "3hours"
 FLOAT_TYPE: "Float64"
 reproducibility_test: true
diff --git a/config/model_configs/single_column_nonorographic_gravity_wave.yml b/config/model_configs/single_column_nonorographic_gravity_wave.yml
@@ -4,4 +4,5 @@ t_end: "1500secs"
 config: "column"
 hyperdiff: false
 dt: "400secs"
+disable_surface_flux_tendency: true
 non_orographic_gravity_wave: true
diff --git a/config/model_configs/single_column_radiative_equilibrium_allsky_idealized_clouds.yml b/config/model_configs/single_column_radiative_equilibrium_allsky_idealized_clouds.yml
@@ -12,6 +12,7 @@ idealized_clouds: true
 config: "column"
 insolation: "timevarying"
 dt_save_to_sol: "30hours"
+disable_surface_flux_tendency: true
 rad: "allskywithclear"
 diagnostics:
   - short_name: rhoa
diff --git a/config/model_configs/single_column_radiative_equilibrium_clearsky.yml b/config/model_configs/single_column_radiative_equilibrium_clearsky.yml
@@ -10,4 +10,5 @@ idealized_h2o: true
 t_end: "654days"
 config: "column"
 dt_save_to_sol: "30hours"
+disable_surface_flux_tendency: true
 rad: "clearsky"
diff --git a/config/model_configs/single_column_radiative_equilibrium_clearsky_prognostic_surface_temp.yml b/config/model_configs/single_column_radiative_equilibrium_clearsky_prognostic_surface_temp.yml
@@ -11,5 +11,6 @@ dt: "3hours"
 dt_rad: "3hours"
 dt_save_to_sol: "30hours"
 dt_save_state_to_disk: "100days"
+disable_surface_flux_tendency: true
 prognostic_surface: true
 toml: [toml/single_column_radiative_equilibrium_clearsky_prognostic_surface_temp.toml]
diff --git a/config/model_configs/single_column_radiative_equilibrium_gray.yml b/config/model_configs/single_column_radiative_equilibrium_gray.yml
@@ -10,6 +10,7 @@ t_end: "654days"
 dz_bottom: 30.0
 config: "column"
 dt_save_to_sol: "30hours"
+disable_surface_flux_tendency: true
 rad: "gray"
 # [2, 2, 80] instead of [1, 1, 80] because Julia ranges are inclusive of the
 # extrema. Given that our columns are 3D, we cannot map the horizontal dimension
diff --git a/config/model_configs/sphere_nonorographic_gravity_wave.yml b/config/model_configs/sphere_nonorographic_gravity_wave.yml
@@ -1,3 +1,4 @@
 t_end: "1500secs"
 dt: "400secs"
+disable_surface_flux_tendency: true
 non_orographic_gravity_wave: true
diff --git a/config/mpi_configs/mpi_make_restart.yml b/config/mpi_configs/mpi_make_restart.yml
@@ -1,2 +1,3 @@
 dt_save_state_to_disk: "5days"
 initial_condition: "DryBaroclinicWave"
+disable_surface_flux_tendency: true
diff --git a/config/mpi_configs/restart_mpi_baroclinic_wave.yml b/config/mpi_configs/restart_mpi_baroclinic_wave.yml
@@ -1,3 +1,4 @@
 t_end: "20days"
 initial_condition: "DryBaroclinicWave"
 restart_file: "mpi_make_restart/output_active/day5.0.hdf5"
+disable_surface_flux_tendency: true
diff --git a/config/perf_configs/bm_baroclinic_wave_moist.yml b/config/perf_configs/bm_baroclinic_wave_moist.yml
@@ -7,4 +7,5 @@ dt_save_to_sol: "Inf"
 log_progress: false
 moist: "equil"
 initial_condition: "MoistBaroclinicWave"
+disable_surface_flux_tendency: true
 precip_model: "0M"
diff --git a/perf/benchmark.jl b/perf/benchmark.jl
@@ -59,7 +59,7 @@ are_boundschecks_forced = Base.JLOptions().check_bounds == 1
         @test compare_mem(trials, "Wfact", 0)
         @test compare_mem(trials, "ldiv!", 0)
         @test compare_mem(trials, "T_imp!", 1000000000000000000000)
-        @test compare_mem(trials, "T_exp_T_lim!", 9920)
+        @test compare_mem(trials, "T_exp_T_lim!", 10496)
         @test compare_mem(trials, "lim!", 0)
         @test compare_mem(trials, "dss!", 0)
         @test compare_mem(trials, "cache!", 120)
diff --git a/reproducibility_tests/ref_counter.jl b/reproducibility_tests/ref_counter.jl
diff --git a/src/ClimaAtmos.jl b/src/ClimaAtmos.jl
diff --git a/src/parameters/create_parameters.jl b/src/parameters/create_parameters.jl
diff --git a/src/prognostic_equations/edmfx_sgs_flux.jl b/src/prognostic_equations/edmfx_sgs_flux.jl
diff --git a/src/prognostic_equations/implicit/implicit_solver.jl b/src/prognostic_equations/implicit/implicit_solver.jl
diff --git a/src/prognostic_equations/remaining_tendency.jl b/src/prognostic_equations/remaining_tendency.jl
diff --git a/src/prognostic_equations/surface_flux.jl b/src/prognostic_equations/surface_flux.jl
diff --git a/src/prognostic_equations/vertical_diffusion_boundary_layer.jl b/src/prognostic_equations/vertical_diffusion_boundary_layer.jl
diff --git a/src/solver/model_getters.jl b/src/solver/model_getters.jl
diff --git a/src/solver/type_getters.jl b/src/solver/type_getters.jl
diff --git a/src/solver/types.jl b/src/solver/types.jl
