diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-04-10 17:54:10 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-04-10 17:54:10 +0000 |
commit | 62681dd7b5e2fbeb6fc55ce05c446ffec617de08 (patch) | |
tree | d35952f16805da4b6b9e13a9bcef8de965add0c3 | |
parent | 8d7fce6a38570b9f5ef2ea637d46ff8effb4f135 (diff) |
GRHydro: remove superfluous "divb(i,j,k) +" from divb computation
From: Roland Haas <rhaas@tapir.caltech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@505 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
-rw-r--r-- | src/GRHydro_Analysis.F90 | 4 |
1 files changed, 2 insertions, 2 deletions
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 |