Try dss on Covariant123Vector

CliMA · Sep 8, 2023 · d805c04 · d805c04
1 parent 4205c01
commit d805c04
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/prognostic_equations/hyperdiffusion.jl b/src/prognostic_equations/hyperdiffusion.jl
@@ -14,6 +14,7 @@ function hyperdiffusion_cache(Y, atmos, do_dss)
     # Grid scale quantities
     ᶜ∇²u = similar(Y.c, C123{FT})
     gs_quantities = (;
+        ᶜ∇²u = similar(Y.c, C123{FT}),
         ᶜ∇²uₕ = similar(Y.c, C12{FT}),
         ᶜ∇²uᵥ = similar(Y.c, C3{FT}),
         ᶜ∇²specific_energy = similar(Y.c, FT),
@@ -133,11 +134,14 @@ function hyperdiffusion_tendency!(Yₜ, Y, p, t)
                 # DSS on Grid scale quantities
                 # Need to split the DSS computation here, because our DSS 
                 # operations do not accept Covariant123Vector types
+                #=
                 @. ᶜ∇²uₕ = C12(ᶜ∇²u)
                 @. ᶜ∇²uᵥ = C3(ᶜ∇²u)
                 dss_op!(ᶜ∇²uₕ, buffer.ᶜ∇²uₕ)
                 dss_op!(ᶜ∇²uᵥ, buffer.ᶜ∇²uᵥ)
                 @. ᶜ∇²u = C123(ᶜ∇²uₕ) + C123(ᶜ∇²uᵥ)
+                =#
+                dss_op!(ᶜ∇²u, buffer.ᶜ∇²u)
                 dss_op!(ᶜ∇²specific_energy, buffer.ᶜ∇²specific_energy)
                 if turbconv_model isa EDMFX && diffuse_tke
                     dss_op!(ᶜ∇²tke⁰, buffer.ᶜ∇²tke⁰)