From a3ab66f6fa472da7c57dae954d8420d5c9395d2a Mon Sep 17 00:00:00 2001
From: Charles Kawczynski <kawczynski.charles@gmail.com>
Date: Wed, 20 Dec 2023 16:59:00 -0500
Subject: [PATCH] Add hyperviscosity conversion for CAM-SE model

---
 src/solver/model_getters.jl | 13 +++++++++++++
 1 file changed, 13 insertions(+)

diff --git a/src/solver/model_getters.jl b/src/solver/model_getters.jl
index be975b0126..ceb1b63d24 100644
--- a/src/solver/model_getters.jl
+++ b/src/solver/model_getters.jl
@@ -45,6 +45,19 @@ 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)^3) # ν_vort
+        κ₄_tracer = FT(0.75 * (30 / Ne * 1.1 * 10^5)^3) # ν_q
+        divergence_damping_factor = FT(5)
+        ClimaHyperdiffusion(;
+            κ₄_vorticity,
+            κ₄_tracer,
+            divergence_damping_factor,
+        )
     elseif hyperdiff_name in ("none", "false", false)
         nothing
     else