Bug fixes from Hyun's code review

Co-authored-by: Hyun (Hyun-Gyu) Kang <[email protected]>
E3SM-Ocean-Discussion · Dec 3, 2024 · 41827e8 · 41827e8
1 parent 948241d
commit 41827e8
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/components/mpas-ocean/src/shared/mpas_ocn_manufactured_solution.F b/components/mpas-ocean/src/shared/mpas_ocn_manufactured_solution.F
@@ -187,13 +187,13 @@ subroutine ocn_manufactured_solution_tend_vel(tend, err)!{{{
             if (.not. config_disable_vel_hmix) then
                if (config_use_mom_del2) then 
                   viscDel2 = config_mom_del2
-                  u = u + viscDel2 * eta0 * kx**2 * cos(phase)
-                  v = v + viscDel2 * eta0 * ky**2 * cos(phase)
+                  u = u + viscDel2 * eta0 * (kx**2 + ky**2) * cos(phase)
+                  v = v + viscDel2 * eta0 * (kx**2 + ky**2) * cos(phase)
                endif
                if (config_use_mom_del4) then 
                   viscDel4 = config_mom_del4
-                  u = u - viscDel4 * eta0 * (kx**4 * cos(phase) + kx**2 * ky**2 * cos(phase))
-                  v = v - viscDel4 * eta0 * (ky**4 * cos(phase) + kx**2 * ky**2 * cos(phase))
+                  u = u - viscDel4 * eta0 * ((kx**4 + ky**4 + kx**2 * ky**2) * cos(phase))
+                  v = v - viscDel4 * eta0 * ((kx**4 + ky**4 + kx**2 * ky**2) * cos(phase))
                endif
             endif