diff options
Diffstat (limited to 'src/GRHydro_Source.F90')
-rw-r--r-- | src/GRHydro_Source.F90 | 93 |
1 files changed, 0 insertions, 93 deletions
diff --git a/src/GRHydro_Source.F90 b/src/GRHydro_Source.F90 index 7b04b8c..c72ea5e 100644 --- a/src/GRHydro_Source.F90 +++ b/src/GRHydro_Source.F90 @@ -454,99 +454,6 @@ 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 |