[Commits] [svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 209)

[cott at tapir.caltech.edu]

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),&