Add hyperviscosity conversion for CAM-SE model

CliMA · Dec 20, 2023 · 2f45d82 · 2f45d82
1 parent dd2dc95
commit 2f45d82
Showing 1 changed file with 12 additions and 0 deletions.
diff --git a/src/solver/model_getters.jl b/src/solver/model_getters.jl
@@ -45,6 +45,18 @@ function get_hyperdiffusion_model(parsed_args, ::Type{FT}) where {FT}
             κ₄_tracer,
             divergence_damping_factor,
         )
+    elseif hyperdiff_name in ("CAM_SE",)
+        # To match hyperviscosity coefficients in:
+        #    https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.1029/2017MS001257
+        #    for equation A18 and A19
+        Ne = parsed_args["h_elem"]
+        κ₄_vorticity = FT(0.15 * (30 / Ne * 1.1 * 10^5)) # ν_vort
+        κ₄_tracer = FT(0.15 * (30 / Ne * 1.1 * 10^5)) # ν_q
+        ClimaHyperdiffusion(;
+            κ₄_vorticity,
+            κ₄_tracer,
+            divergence_damping_factor,
+        )
     elseif hyperdiff_name in ("none", "false", false)
         nothing
     else