Make cvmix compatible with maxLevelCell = 0

E3SM-Ocean-Discussion · Mar 10, 2024 · 47391d1 · 47391d1
1 parent a85a637
commit 47391d1
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/components/mpas-ocean/src/shared/mpas_ocn_vmix_cvmix.F b/components/mpas-ocean/src/shared/mpas_ocn_vmix_cvmix.F
@@ -786,6 +786,7 @@ subroutine ocn_vmix_coefs_cvmix_build(meshPool, statePool, forcingPool, err, tim
          !$omp parallel
          !$omp do schedule(runtime) private(nEdges, areaSum, edgeCount, iEdge, iNeighbor)
          do iCell=1,nCells
+            if (maxLevelCell(iCell) == 0) cycle
             nEdges = nEdgesOnCell(iCell)
             boundaryLayerDepthSmooth(iCell) = 0.0_RKIND
             areaSum = 0.0_RKIND