diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:45 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-01-14 14:23:45 +0000 |
commit | 5e3fd99e900f957b320e9dd428fa062fc1713980 (patch) | |
tree | 5210d699194b7cb01f74db97e24aa0855c5d0b47 /src | |
parent | 9a0acd9e6addbcdbbe2905b5a986664ec646a861 (diff) |
GRHydro: clear divB everywhere rather than just in interior
From: Roland Haas <roland.haas@physics.gatech.edu>
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@461 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src')
-rw-r--r-- | src/GRHydro_Analysis.F90 | 11 |
1 files changed, 1 insertions, 10 deletions
diff --git a/src/GRHydro_Analysis.F90 b/src/GRHydro_Analysis.F90 index 3bed392..8417fb7 100644 --- a/src/GRHydro_Analysis.F90 +++ b/src/GRHydro_Analysis.F90 @@ -13,16 +13,7 @@ subroutine GRHydro_Analysis_Init(CCTK_ARGUMENTS) 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 |