1 files changed, 100 insertions, 0 deletions
diff --git a/src/GRHydro_Source.F90 b/src/GRHydro_Source.F90
index 6758c90..7b04b8c 100644
--- a/src/GRHydro_Source.F90
+++ b/src/GRHydro_Source.F90
@@ -93,6 +93,11 @@ subroutine SourceTerms(CCTK_ARGUMENTS)
   CCTK_REAL, DIMENSION(:,:,:), POINTER :: beta1, beta2, beta3
   CCTK_REAL, DIMENSION(:,:,:,:), POINTER :: vup
 
+  ! poison hack
+  CCTK_INT, SAVE :: reflevel = -1
+  CCTK_INT, SAVE :: last_iteration_seen = -1
+  CCTK_INT, SAVE :: mol_substep = -1
+
   if (GRHydro_UseGeneralCoordinates(cctkGH).ne.0) then
     g11 => gaa
     g12 => gab
@@ -449,6 +454,101 @@ subroutine SourceTerms(CCTK_ARGUMENTS)
 
   deallocate(force_spatial_second_order)
 
+#if(0) /* poison edges of domain */
+  if(last_iteration_seen .ne. cctk_iteration .or. reflevel .ne. grhydro_reflevel) then
+    last_iteration_seen = cctk_iteration
+    reflevel = grhydro_reflevel
+    mol_substep = 0
+  else
+    mol_substep = mol_substep + 1
+  end if
+  do k = 1, GRHydro_stencil*mol_substep
+    do j = 1, cctk_lsh(2)
+      do i = 1, cctk_lsh(1)
+        dens(i,j,k) = -1d100
+        Scon(i,j,k,1) = -1d100
+        Scon(i,j,k,2) = -1d100
+        Scon(i,j,k,3) = -1d100
+        tau(i,j,k) = -1d100
+        Bcons(i,j,k,1) = -1d100
+        Bcons(i,j,k,2) = -1d100
+        Bcons(i,j,k,3) = -1d100
+      end do
+    end do
+  end do
+  do k = cctk_lsh(3)-GRHydro_stencil*mol_substep+1, cctk_lsh(3)
+    do j = 1, cctk_lsh(2)
+      do i = 1, cctk_lsh(1)
+        dens(i,j,k) = -1d100
+        Scon(i,j,k,1) = -1d100
+        Scon(i,j,k,2) = -1d100
+        Scon(i,j,k,3) = -1d100
+        tau(i,j,k) = -1d100
+        Bcons(i,j,k,1) = -1d100
+        Bcons(i,j,k,2) = -1d100
+        Bcons(i,j,k,3) = -1d100
+      end do
+    end do
+  end do
+  do i = 1, GRHydro_stencil*mol_substep
+    do k = 1, cctk_lsh(3)
+      do j = 1, cctk_lsh(2)
+        dens(i,j,k) = -1d100
+        Scon(i,j,k,1) = -1d100
+        Scon(i,j,k,2) = -1d100
+        Scon(i,j,k,3) = -1d100
+        tau(i,j,k) = -1d100
+        Bcons(i,j,k,1) = -1d100
+        Bcons(i,j,k,2) = -1d100
+        Bcons(i,j,k,3) = -1d100
+      end do
+    end do
+  end do
+  do i = cctk_lsh(1)-GRHydro_stencil*mol_substep+1, cctk_lsh(1)
+    do k = 1, cctk_lsh(3)
+      do j = 1, cctk_lsh(2)
+        dens(i,j,k) = -1d100
+        Scon(i,j,k,1) = -1d100
+        Scon(i,j,k,2) = -1d100
+        Scon(i,j,k,3) = -1d100
+        tau(i,j,k) = -1d100
+        Bcons(i,j,k,1) = -1d100
+        Bcons(i,j,k,2) = -1d100
+        Bcons(i,j,k,3) = -1d100
+      end do
+    end do
+  end do
+  do j = 1, GRHydro_stencil*mol_substep
+    do i = 1, cctk_lsh(1)
+      do k = 1, cctk_lsh(3)
+        dens(i,j,k) = -1d100
+        Scon(i,j,k,1) = -1d100
+        Scon(i,j,k,2) = -1d100
+        Scon(i,j,k,3) = -1d100
+        tau(i,j,k) = -1d100
+        Bcons(i,j,k,1) = -1d100
+        Bcons(i,j,k,2) = -1d100
+        Bcons(i,j,k,3) = -1d100
+      end do
+    end do
+  end do
+  do j = cctk_lsh(2)-GRHydro_stencil*mol_substep+1, cctk_lsh(2)
+    do i = 1, cctk_lsh(1)
+      do k = 1, cctk_lsh(3)
+        dens(i,j,k) = -1d100
+        Scon(i,j,k,1) = -1d100
+        Scon(i,j,k,2) = -1d100
+        Scon(i,j,k,3) = -1d100
+        tau(i,j,k) = -1d100
+        Bcons(i,j,k,1) = -1d100
+        Bcons(i,j,k,2) = -1d100
+        Bcons(i,j,k,3) = -1d100
+      end do
+    end do
+  end do
+#endif
+  
+
 end subroutine SourceTerms