Skip to content

Commit f881969

Browse files
committed
Set vertical velocity top to 0 from k=1:minLevelCell-1
1 parent c85f36c commit f881969

File tree

1 file changed

+1
-0
lines changed

1 file changed

+1
-0
lines changed

src/core_ocean/shared/mpas_ocn_diagnostics.F

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1194,6 +1194,7 @@ subroutine ocn_diagnostic_solve_vortVel(relativeVorticity, &
11941194
end do
11951195
end do
11961196
! Vertical velocity at bottom (maxLevelCell(iCell)+1) is zero, initialized above.
1197+
vertVelocityTop(1:minLevelCell(iCell)-1, iCell) = 0.0_RKIND
11971198
vertVelocityTop(maxLevelCell(iCell)+1, iCell) = 0.0_RKIND
11981199
do k = maxLevelCell(iCell), 1, -1
11991200
vertVelocityTop(k,iCell) = vertVelocityTop(k+1,iCell) - div_hu(k)

0 commit comments

Comments
 (0)