diff options
| author | cott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-01-08 01:37:53 +0000 |
|---|---|---|
| committer | cott <cott@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2011-01-08 01:37:53 +0000 |
| commit | 15ba4eb7096d838ad7601731ce6216cfa8818283 (patch) | |
| tree | 6522dedbb2b4b4c07343c24042ec856a2a77df5e | |
| parent | f2cdc7f3a147f13dd3f4861c101f995f8cbfbd05 (diff) | |
* fix a very nasty bug in PPM reconstruction that sneeked in with the very recent
PPM changes. Old version of the code were not affected.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@209 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
| -rw-r--r-- | src/GRHydro_Reconstruct.F90 | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/GRHydro_Reconstruct.F90 b/src/GRHydro_Reconstruct.F90 index 23ec0a2..b1af736 100644 --- a/src/GRHydro_Reconstruct.F90 +++ b/src/GRHydro_Reconstruct.F90 @@ -361,7 +361,7 @@ subroutine Reconstruction(CCTK_ARGUMENTS) if(clean_divergence.ne.0) then !$OMP PARALLEL DO PRIVATE(i, j) do k = GRHydro_stencil, nz - GRHydro_stencil + 1 - do j = GRHydro_stencil, ny - GRHydro_stencil + 1 + do j = GRHydro_stencil, nx - GRHydro_stencil + 1 call SimplePPM_1dM(GRHydro_eos_handle,0,ny,CCTK_DELTA_SPACE(2),& rho(j,:,k),vely(j,:,k),velz(j,:,k),velx(j,:,k),& Bvecy(j,:,k),Bvecz(j,:,k),Bvecx(j,:,k),psidc(j,:,k),eps(j,:,k),press(j,:,k),& @@ -389,7 +389,7 @@ subroutine Reconstruction(CCTK_ARGUMENTS) else !clean_divergence !$OMP PARALLEL DO PRIVATE(i, j) do k = GRHydro_stencil, nz - GRHydro_stencil + 1 - do j = GRHydro_stencil, ny - GRHydro_stencil + 1 + do j = GRHydro_stencil, nx - GRHydro_stencil + 1 call SimplePPM_1dM(GRHydro_eos_handle,0,ny,CCTK_DELTA_SPACE(2),& rho(j,:,k),vely(j,:,k),velz(j,:,k),velx(j,:,k),& Bvecy(j,:,k),Bvecz(j,:,k),Bvecx(j,:,k),dum(j,:,k),eps(j,:,k),press(j,:,k),& @@ -418,7 +418,7 @@ subroutine Reconstruction(CCTK_ARGUMENTS) else !evolve_mhd !$OMP PARALLEL DO PRIVATE(i, j) do k = GRHydro_stencil, nz - GRHydro_stencil + 1 - do j = GRHydro_stencil, ny - GRHydro_stencil + 1 + do j = GRHydro_stencil, nx - GRHydro_stencil + 1 call SimplePPM_1d(GRHydro_eos_handle,0,ny,CCTK_DELTA_SPACE(2),& rho(j,:,k),vely(j,:,k),velz(j,:,k),velx(j,:,k),eps(j,:,k),& press(j,:,k),rhominus(j,:,k),velyminus(j,:,k),velzminus(j,:,k),& |