diff options
| -rw-r--r-- | src/findnormals.F90 | 6 |
1 files changed, 5 insertions, 1 deletions
diff --git a/src/findnormals.F90 b/src/findnormals.F90 index d6f2fbe..2936424 100644 --- a/src/findnormals.F90 +++ b/src/findnormals.F90 @@ -209,7 +209,7 @@ subroutine excision_findnormals (ierr, mask, dirx, diry, dirz, ni, nj, nk) if (mask(i,j,k)==MASK_BOUNDARY) then ! make a snapshot - msk = mask(i-1:i+1,j-1:j+1,k-1:k+1) == MASK_ACTIVE + msk = mask(i-1:i+1,j-1:j+1,k-1:k+1) /= MASK_EXCISED ! for a template to match, all grid points have to be ! either active, or set in the template. @@ -270,6 +270,10 @@ subroutine excision_findnormals (ierr, mask, dirx, diry, dirz, ni, nj, nk) dir(:) = -dir(:) + if (mask(i+dir(1), j+dir(2), k+dir(3)) /= MASK_ACTIVE) then + call CCTK_WARN (0, "Mask boundary layer is too thick") + end if + dirx(i,j,k) = dir(1) diry(i,j,k) = dir(2) dirz(i,j,k) = dir(3) |