[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 505)

[rhaas at tapir.caltech.edu]

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 at 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