diff options
Diffstat (limited to 'src/findnormals.F90')
| -rw-r--r-- | src/findnormals.F90 | 4 |
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 |