User: rhaas Date: 2013/04/10 12:54 PM

Modified: /trunk/src/ GRHydro_Analysis.F90

Log: GRHydro: remove superfluous "divb(i,j,k) +" from divb computation

From: Roland Haas rhaas@tapir.caltech.edu

File Changes:

Directory: /trunk/src/ ======================

File [modified]: GRHydro_Analysis.F90 Delta lines: +2 -2 =================================================================== --- trunk/src/GRHydro_Analysis.F90 2013-04-10 17:54:08 UTC (rev 504) +++ trunk/src/GRHydro_Analysis.F90 2013-04-10 17:54:10 UTC (rev 505) @@ -41,7 +41,7 @@ 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 @@ 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