User: cott Date: 2011/01/07 07:37 PM

Modified: /trunk/src/ GRHydro_Reconstruct.F90

Log: * 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.

File Changes:

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

File [modified]: GRHydro_Reconstruct.F90 Delta lines: +3 -3 =================================================================== --- trunk/src/GRHydro_Reconstruct.F90 2011-01-06 01:13:02 UTC (rev 208) +++ trunk/src/GRHydro_Reconstruct.F90 2011-01-08 01:37:53 UTC (rev 209) @@ -361,7 +361,7 @@ 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 @@ 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 @@ 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),&