diff options
Diffstat (limited to 'src/extrapolate.F90')
-rw-r--r-- | src/extrapolate.F90 | 11 |
1 files changed, 7 insertions, 4 deletions
diff --git a/src/extrapolate.F90 b/src/extrapolate.F90 index 19ca1b6..c09e35c 100644 --- a/src/extrapolate.F90 +++ b/src/extrapolate.F90 @@ -5,9 +5,8 @@ ! diry, and dirz for extrapolation. #include "cctk.h" -#include "cctk_Parameters.h" - #include "cctk_Functions.h" +#include "cctk_Parameters.h" #include "maskvalues.h" @@ -46,15 +45,19 @@ subroutine excision_extrapolate (ierr, var, oldvar, & integer i,j,k integer ii,jj,kk + if (CCTK_IsThornActive(CCTK_THORNSTRING) == 0) then + call CCTK_WARN (0, "The routine excision_extrapolate was called, but thorn " // CCTK_THORNSTRING // " is not active") + end if + do k=2,nk-1 do j=2,nj-1 do i=2,ni-1 - if (mask(i,j,k)==MASK_BOUNDARY) then + if (abs(mask(i,j,k)-MASK_BOUNDARY)<0.01) then ii = i + int(dirx(i,j,k)) jj = j + int(diry(i,j,k)) kk = k + int(dirz(i,j,k)) var(i,j,k) = oldvar(i,j,k) + var(ii,jj,kk) - oldvar(ii,jj,kk) - else if (mask(i,j,k)==MASK_EXCISED) then + else if (abs(mask(i,j,k)-MASK_EXCISED)<0.01) then var(i,j,k) = var0 end if end do |