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