From 62681dd7b5e2fbeb6fc55ce05c446ffec617de08 Mon Sep 17 00:00:00 2001 From: rhaas Date: Wed, 10 Apr 2013 17:54:10 +0000 Subject: GRHydro: remove superfluous "divb(i,j,k) +" from divb computation From: Roland Haas git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@505 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45 --- src/GRHydro_Analysis.F90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'src/GRHydro_Analysis.F90') diff --git a/src/GRHydro_Analysis.F90 b/src/GRHydro_Analysis.F90 index 8417fb7..cbb6f6a 100644 --- a/src/GRHydro_Analysis.F90 +++ b/src/GRHydro_Analysis.F90 @@ -41,7 +41,7 @@ subroutine GRHydro_CalcDivB(CCTK_ARGUMENTS) if(track_divB.ne.0) then if(transport_constraints.ne.0) then ! edge based divergence (see WhiskyMHD & Bruno's thesis, Eq. 7.27) - divB(i,j,k) = divB(i,j,k) + & + divB(i,j,k) = & 0.25d0*(Bcons(i+1,j,k,1)-Bcons(i,j,k,1)+ & Bcons(i+1,j+1,k,1)-Bcons(i,j+1,k,1)+ & Bcons(i+1,j,k+1,1)-Bcons(i,j,k+1,1)+ & @@ -68,7 +68,7 @@ subroutine GRHydro_CalcDivB(CCTK_ARGUMENTS) Bcons_r3 = 0.5d0 * (Bcons(i,j,k,3) + & Bcons(i,j,k+1,3)) - divB(i,j,k) = divB(i,j,k) + (Bcons_l1 - Bcons_r1) * idx + (Bcons_l2 - Bcons_r2) * idy + (Bcons_l3 - Bcons_r3) * idz + divB(i,j,k) = (Bcons_l1 - Bcons_r1) * idx + (Bcons_l2 - Bcons_r2) * idy + (Bcons_l3 - Bcons_r3) * idz endif endif -- cgit v1.2.3