From a291b4379f0e35ee9616ebb895944bc22b569545 Mon Sep 17 00:00:00 2001 From: miguel Date: Fri, 30 Mar 2001 21:37:16 +0000 Subject: Debugging. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@13 f75ba9e5-694f-0410-ac2c-87ea7ce7132b --- src/findnormals.F90 | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) 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) -- cgit v1.2.3