Merge pull request #4132 from CliMA/sa/improve_stability

add Jacobain elements for du3^j_err_dq^j and drhoa_err_dq^j

Author: sajjadazimi

config/model_configs/prognostic_edmfx_rico_implicit_column.yml

@@ -7,7 +7,7 @@ tracer_upwinding: first_order
 implicit_diffusion: true
 implicit_sgs_advection: true
 implicit_sgs_entr_detr: true
-implicit_nh_pressure: true
+implicit_sgs_nh_pressure: true
 implicit_sgs_vertdiff: true
 implicit_sgs_mass_flux: true
 approximate_linear_solve_iters: 2
@@ -31,7 +31,7 @@ y_elem: 2
 z_elem: 100
 z_stretch: false
 perturb_initstate: false
-dt: "60secs"
+dt: "80secs"
 t_end: "24hours"
 dt_save_state_to_disk: "60mins"
 toml: [toml/prognostic_edmfx_1M.toml]

config/model_configs/prognostic_edmfx_trmm_implicit_column.yml

@@ -7,7 +7,7 @@ tracer_upwinding: first_order
 implicit_diffusion: true
 implicit_sgs_advection: true
 implicit_sgs_entr_detr: true
-implicit_nh_pressure: true
+implicit_sgs_nh_pressure: true
 implicit_sgs_vertdiff: true
 implicit_sgs_mass_flux: true
 approximate_linear_solve_iters: 2
@@ -31,7 +31,7 @@ x_elem: 2
 y_elem: 2
 z_elem: 82
 z_stretch: false
-dt: 50secs
+dt: 150secs
 t_end: 12hours
 dt_save_state_to_disk: 60mins
 toml: [toml/prognostic_edmfx_1M.toml]

reproducibility_tests/ref_counter.jl

@@ -1,4 +1,4 @@
-277
+278
 
 # **README**
 #
@@ -20,6 +20,9 @@
 
 
 #=
+278
+- Add ∂/∂q elements to Jacobian
+
 277
 - Update to use [email protected]. Supports Charnock-parameterization for aerodynamic roughness (default is still user-prescribed ScalarRoughness). SurfaceFluxes catch for neutrally-stable boundary layers (ζ ≈ 0) removed.