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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/01/14 08:23 AM

Modified:
 /trunk/src/
  GRHydro_Analysis.F90

Log:
 GRHydro: clear divB everywhere rather than just in interior
 
 From: Roland Haas <roland.haas at physics.gatech.edu>

File Changes:

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

File [modified]: GRHydro_Analysis.F90
Delta lines: +1 -10
===================================================================
--- trunk/src/GRHydro_Analysis.F90	2013-01-14 14:23:44 UTC (rev 460)
+++ trunk/src/GRHydro_Analysis.F90	2013-01-14 14:23:45 UTC (rev 461)
@@ -13,16 +13,7 @@
   DECLARE_CCTK_PARAMETERS
   DECLARE_CCTK_FUNCTIONS
 
-  CCTK_INT :: i,j,k
-  !$OMP PARALLEL DO PRIVATE(i,j,k)
-  do k = GRHydro_stencil + 1 - transport_constraints, cctk_lsh(3) - GRHydro_stencil ! we need to compute Evec on all faces/edges where the fluxes are defined
-     do j = GRHydro_stencil + 1 - transport_constraints, cctk_lsh(2) - GRHydro_stencil
-        do i = GRHydro_stencil + 1 - transport_constraints, cctk_lsh(1) - GRHydro_stencil
-           divB(i,j,k) = 0.0d0
-        end do
-     end do
-  end do
-  !$OMP END PARALLEL DO
+  divB = 0.0d0
 
 end subroutine GRHydro_Analysis_Init