Initialize normalVelocity=0 at the top of the domain

- Initialize normalVelocity=0 at the top of the domain to fix the restart issue
E3SM-Project · Oct 26, 2023 · 54a3220 · 54a3220
1 parent 755247f
commit 54a3220
Showing 1 changed file with 8 additions and 3 deletions.
diff --git a/components/mpas-ocean/src/mode_forward/mpas_ocn_time_integration_split_ab2.F b/components/mpas-ocean/src/mode_forward/mpas_ocn_time_integration_split_ab2.F
@@ -3238,10 +3238,15 @@ subroutine ocn_time_integration_split_ab2_init(domain)!{{{
          call mpas_pool_get_array(meshPool, 'minLevelEdgeBot', minLevelEdgeBot)
          call mpas_pool_get_array(meshPool, 'maxLevelEdgeTop', maxLevelEdgeTop)
 
+         ! normalVelocity=0 on inactive cells
          do iEdge = 1, nEdges
-            kmax  = maxLevelEdgeTop(iEdge)
-            ! normalVelocity=0 on land cells
-            do k = kmax+1, nVertLevels
+            ! at top
+            do k = 1, minLevelEdgeBot(iEdge)-1
+               normalVelocity(k,iEdge) = 0.0_RKIND
+            enddo
+
+            ! at bottom
+            do k = maxLevelEdgeTop(iEdge)+1, nVertLevels
                normalVelocity(k,iEdge) = 0.0_RKIND
             enddo
          enddo ! edge loop