diff options
| -rw-r--r-- | src/findboundary.F90 | 9 |
1 files changed, 3 insertions, 6 deletions
diff --git a/src/findboundary.F90 b/src/findboundary.F90 index 1e16e97..c729fe8 100644 --- a/src/findboundary.F90 +++ b/src/findboundary.F90 @@ -8,7 +8,7 @@ #include "maskvalues.h" -subroutine excision_findboundary (ierr, cctkgh, mask, ni, nj, nk) +subroutine excision_findboundary (ierr, mask, ni, nj, nk) implicit none @@ -19,10 +19,7 @@ subroutine excision_findboundary (ierr, cctkgh, mask, ni, nj, nk) ! out: zero for success, nonzero for error integer :: ierr - - ! in: pointer to Cactus GH - CCTK_POINTER :: cctkgh - + ! in: array sizes for grid functions ! (you can pass in cctk_lsh(:) for these) integer :: ni,nj,nk @@ -42,7 +39,7 @@ subroutine excision_findboundary (ierr, cctkgh, mask, ni, nj, nk) do kk=k-1,k+1 do jj=j-1,j+1 do ii=i-1,i+1 - bnd = bnd .or. mask(i,j,k)/=MASK_EXCISED + bnd = bnd .or. mask(ii,jj,kk)==MASK_ACTIVE end do end do end do |