diff options
| author | miguel <miguel@f75ba9e5-694f-0410-ac2c-87ea7ce7132b> | 2001-03-30 21:37:16 +0000 |
|---|---|---|
| committer | miguel <miguel@f75ba9e5-694f-0410-ac2c-87ea7ce7132b> | 2001-03-30 21:37:16 +0000 |
| commit | a291b4379f0e35ee9616ebb895944bc22b569545 (patch) | |
| tree | bf83c1505617499f6def40c884c39c9c77656854 | |
| parent | ced01f2c056ba69da517666210a3663c4860e7ab (diff) | |
Debugging.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@13 f75ba9e5-694f-0410-ac2c-87ea7ce7132b
| -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) |