User: rhaas Date: 2012/07/17 12:08 PM

Modified: /trunk/src/ GRHydro_UpdateMask.F90

Log: GRHydro: only loop over points where rhs is defined for UpdateMask

From: Roland Haas roland.haas@physics.gatech.edu

File Changes:

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

File [modified]: GRHydro_UpdateMask.F90 Delta lines: +6 -6 =================================================================== --- trunk/src/GRHydro_UpdateMask.F90 2012-07-17 17:08:35 UTC (rev 402) +++ trunk/src/GRHydro_UpdateMask.F90 2012-07-17 17:08:37 UTC (rev 403) @@ -38,9 +38,9 @@

if(evolve_temper.ne.1.and.evolve_Y_e.ne.1) then !$OMP PARALLEL DO PRIVATE(k,j,i) - do k = 1, cctk_lsh(3) - do j = 1, cctk_lsh(2) - do i = 1, cctk_lsh(1) + do k = 1+GRHydro_stencil, cctk_lsh(3)-GRHydro_stencil + do j = 1+GRHydro_stencil, cctk_lsh(2)-GRHydro_stencil + do i = 1+GRHydro_stencil, cctk_lsh(1)-GRHydro_stencil if ( GRHydro_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0) .or. & (atmosphere_mask(i,j,k) .ne. 0) .or. & (tau(i,j,k) + frac * taurhs(i,j,k) .le. 0.d0) .or. & @@ -59,9 +59,9 @@ else ! allow negative total energy !$OMP PARALLEL DO PRIVATE(k,j,i) - do k = 1, cctk_lsh(3) - do j = 1, cctk_lsh(2) - do i = 1, cctk_lsh(1) + do k = 1+GRHydro_stencil, cctk_lsh(3)-GRHydro_stencil + do j = 1+GRHydro_stencil, cctk_lsh(2)-GRHydro_stencil + do i = 1+GRHydro_stencil, cctk_lsh(1)-GRHydro_stencil if ( GRHydro_enable_internal_excision /= 0 .and. (hydro_excision_mask(i,j,k) .ne. 0) .or. & (atmosphere_mask(i,j,k) .ne. 0) .or. & (dens(i,j,k) + frac * densrhs(i,j,k) .le. 0.d0) ) then