diff options
Diffstat (limited to 'src/GRHydro_CalcUpdate.F90')
-rw-r--r-- | src/GRHydro_CalcUpdate.F90 | 10 |
1 files changed, 2 insertions, 8 deletions
diff --git a/src/GRHydro_CalcUpdate.F90 b/src/GRHydro_CalcUpdate.F90 index 1b3b96d..f7efaa0 100644 --- a/src/GRHydro_CalcUpdate.F90 +++ b/src/GRHydro_CalcUpdate.F90 @@ -59,7 +59,7 @@ subroutine UpdateCalculation(CCTK_ARGUMENTS) do k = GRHydro_stencil + 1, cctk_lsh(3) - GRHydro_stencil do j = GRHydro_stencil + 1, cctk_lsh(2) - GRHydro_stencil do i = GRHydro_stencil + 1, cctk_lsh(1) - GRHydro_stencil - + alp_l = 0.5d0 * (alp(i,j,k) + & alp(i-xoffset,j-yoffset,k-zoffset)) alp_r = 0.5d0 * (alp(i,j,k) + & @@ -103,12 +103,6 @@ subroutine UpdateCalculation(CCTK_ARGUMENTS) divB(i,j,k) = divB(i,j,k) + ( alp_l * Bvec_l - alp_r * Bvec_r ) * idx endif -!!$ if(i.eq.75.and.j.eq.5.and.k.eq.5)write(6,*)'update:', & -!!$ Bconsrhs(i,j,k,1),Bconsxflux(i,j,k), & -!!$ Bconsxflux(i-xoffset,j-yoffset,k-zoffset), & -!!$ Bconsrhs(i,j,k,3),Bconszflux(i,j,k), & -!!$ Bconszflux(i-xoffset,j-yoffset,k-zoffset) - endif if (evolve_tracer .ne. 0) then @@ -270,7 +264,7 @@ subroutine UpdateCalculation(CCTK_ARGUMENTS) divB(i,j,k) = divB(i,j,k) + ( Bvec_l - Bvec_r ) * idx endif endif - + enddo enddo enddo |