From eb1d76b0841fa77d7158ec8392d5dcfbfba8949e Mon Sep 17 00:00:00 2001 From: rhaas Date: Fri, 4 Nov 2011 18:34:36 +0000 Subject: second iteration of constraint transport * fix some indices * move poison loop to proper spot * compute conserved divergence in track_divB * add Maple worksheet to check constraint transport indices git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@298 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45 --- src/GRHydro_Source.F90 | 93 -------------------------------------------------- 1 file changed, 93 deletions(-) (limited to 'src/GRHydro_Source.F90') 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 -- cgit v1.2.3