Fixing small bug.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@12 f75ba9e5-694f-0410-ac2c-87ea7ce7132b

1 files changed, 3 insertions, 1 deletions

diff --git a/src/findnormals.F90 b/src/findnormals.F90
index d34f25a..d6f2fbe 100644
--- a/src/findnormals.F90
+++ b/src/findnormals.F90
@@ -205,10 +205,11 @@ subroutine excision_findnormals (ierr, mask, dirx, diry, dirz, ni, nj, nk)
   do k=2,nk-1
      do j=2,nj-1
         do i=2,ni-1
+
            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) /= 0
+              msk = mask(i-1:i+1,j-1:j+1,k-1:k+1) == MASK_ACTIVE
               
               ! for a template to match, all grid points have to be
               ! either active, or set in the template.
@@ -280,6 +281,7 @@ subroutine excision_findnormals (ierr, mask, dirx, diry, dirz, ni, nj, nk)
   end do
   
   ! Success
+
   ierr = 0
   
 end subroutine excision_findnormals