diff --git a/src/core/MOM_continuity_PPM.F90 b/src/core/MOM_continuity_PPM.F90 index 86d05a15b4..3ee272122e 100644 --- a/src/core/MOM_continuity_PPM.F90 +++ b/src/core/MOM_continuity_PPM.F90 @@ -237,14 +237,14 @@ subroutine continuity_PPM(u, v, hin, h, uh, vh, dt, G, GV, CS, uhbt, vhbt, OBC, h(i,j,k) = h_input(i+1,j,k) enddo enddo - do J=LB%jsh-1,LB%jeh - do i=LB%ish-1,LB%ieh+1 - if (OBC%OBC_kind_v(i,J) == OBC_FLATHER .and. (OBC%OBC_direction_v(i,J) == OBC_DIRECTION_E)) & - v(i,J,k) = v(i-1,J,k) - if (OBC%OBC_kind_v(i,J) == OBC_FLATHER .and. (OBC%OBC_direction_v(i,J) == OBC_DIRECTION_W)) & - v(i,J,k) = v(i+1,J,k) +! do J=LB%jsh-1,LB%jeh +! do i=LB%ish-1,LB%ieh+1 +! if (OBC%OBC_kind_v(i,J) == OBC_FLATHER .and. (OBC%OBC_direction_v(i,J) == OBC_DIRECTION_E)) & +! v(i,J,k) = v(i-1,J,k) +! if (OBC%OBC_kind_v(i,J) == OBC_FLATHER .and. (OBC%OBC_direction_v(i,J) == OBC_DIRECTION_W)) & +! v(i,J,k) = v(i+1,J,k) enddo - enddo ; enddo +! enddo ; enddo endif LB%ish = G%isc ; LB%ieh = G%iec ; LB%jsh = G%jsc ; LB%jeh = G%jec @@ -272,12 +272,12 @@ subroutine continuity_PPM(u, v, hin, h, uh, vh, dt, G, GV, CS, uhbt, vhbt, OBC, if (OBC%OBC_kind_v(i,J) == OBC_FLATHER .and. (OBC%OBC_direction_v(i,J) == OBC_DIRECTION_S)) & h(i,j,k) = h_input(i,j+1,k) enddo ; enddo - do j=LB%jsh-1,LB%jeh+1 ; do I=LB%ish-1,LB%ieh - if (OBC%OBC_kind_u(I,j) == OBC_FLATHER .and. (OBC%OBC_direction_u(I,j) == OBC_DIRECTION_N)) & - u(I,j,k) = u(I,j-1,k) - if (OBC%OBC_kind_u(I,j) == OBC_FLATHER .and. (OBC%OBC_direction_u(I,j) == OBC_DIRECTION_S)) & - u(I,j,k) = u(I,j+1,k) - enddo ; enddo +! do j=LB%jsh-1,LB%jeh+1 ; do I=LB%ish-1,LB%ieh +! if (OBC%OBC_kind_u(I,j) == OBC_FLATHER .and. (OBC%OBC_direction_u(I,j) == OBC_DIRECTION_N)) & +! u(I,j,k) = u(I,j-1,k) +! if (OBC%OBC_kind_u(I,j) == OBC_FLATHER .and. (OBC%OBC_direction_u(I,j) == OBC_DIRECTION_S)) & +! u(I,j,k) = u(I,j+1,k) +! enddo ; enddo enddo endif else ! .not. x_first @@ -305,12 +305,12 @@ subroutine continuity_PPM(u, v, hin, h, uh, vh, dt, G, GV, CS, uhbt, vhbt, OBC, if (OBC%OBC_kind_v(i,J) == OBC_FLATHER .and. (OBC%OBC_direction_v(i,J) == OBC_DIRECTION_S)) & h(i,j,k) = h_input(i,j+1,k) enddo ; enddo - do j=LB%jsh-1,LB%jeh+1 ; do I=LB%ish-1,LB%ieh - if (OBC%OBC_kind_u(I,j) == OBC_FLATHER .and. (OBC%OBC_direction_u(I,j) == OBC_DIRECTION_N)) & - u(I,j,k) = u(I,j-1,k) - if (OBC%OBC_kind_u(I,j) == OBC_FLATHER .and. (OBC%OBC_direction_u(I,j) == OBC_DIRECTION_S)) & - u(I,j,k) = u(I,j+1,k) - enddo ; enddo +! do j=LB%jsh-1,LB%jeh+1 ; do I=LB%ish-1,LB%ieh +! if (OBC%OBC_kind_u(I,j) == OBC_FLATHER .and. (OBC%OBC_direction_u(I,j) == OBC_DIRECTION_N)) & +! u(I,j,k) = u(I,j-1,k) +! if (OBC%OBC_kind_u(I,j) == OBC_FLATHER .and. (OBC%OBC_direction_u(I,j) == OBC_DIRECTION_S)) & +! u(I,j,k) = u(I,j+1,k) +! enddo ; enddo enddo endif @@ -340,14 +340,14 @@ subroutine continuity_PPM(u, v, hin, h, uh, vh, dt, G, GV, CS, uhbt, vhbt, OBC, h(i,j,k) = h_input(i+1,j,k) enddo enddo - do J=LB%jsh-1,LB%jeh - do i=LB%ish-1,LB%ieh+1 - if (OBC%OBC_kind_v(i,J) == OBC_FLATHER .and. (OBC%OBC_direction_v(i,J) == OBC_DIRECTION_E)) & - v(i,J,k) = v(i-1,J,k) - if (OBC%OBC_kind_v(i,J) == OBC_FLATHER .and. (OBC%OBC_direction_v(i,J) == OBC_DIRECTION_W)) & - v(i,J,k) = v(i+1,J,k) +! do J=LB%jsh-1,LB%jeh +! do i=LB%ish-1,LB%ieh+1 +! if (OBC%OBC_kind_v(i,J) == OBC_FLATHER .and. (OBC%OBC_direction_v(i,J) == OBC_DIRECTION_E)) & +! v(i,J,k) = v(i-1,J,k) +! if (OBC%OBC_kind_v(i,J) == OBC_FLATHER .and. (OBC%OBC_direction_v(i,J) == OBC_DIRECTION_W)) & +! v(i,J,k) = v(i+1,J,k) enddo - enddo ; enddo +! enddo ; enddo endif endif