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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/05/11 01:29 AM

Modified:
 /trunk/src/
  GRHydro_Prim2ConM.F90

Log:
 GRHydro: Fix a problem with ghost zones in Prim2ConCells calls in MHD case.
 
 From: Christian Reisswig <reisswig at tapir.caltech.edu>

File Changes:

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

File [modified]: GRHydro_Prim2ConM.F90
Delta lines: +9 -9
===================================================================
--- trunk/src/GRHydro_Prim2ConM.F90	2013-05-11 06:29:34 UTC (rev 511)
+++ trunk/src/GRHydro_Prim2ConM.F90	2013-05-11 06:29:37 UTC (rev 512)
@@ -464,9 +464,9 @@
   if(evolve_temper.ne.1) then
 
   !$OMP PARALLEL DO PRIVATE(k,j,i,det), REDUCTION(MAX:maxtau0)
-  do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1
-    do j = GRHydro_stencil,cctk_lsh(2)-GRHydro_stencil+1
-      do i = GRHydro_stencil,cctk_lsh(1)-GRHydro_stencil+1
+  do k = 1,cctk_lsh(3)
+    do j = 1,cctk_lsh(2)
+      do i = 1,cctk_lsh(1)
         
          det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k), \
               gyy(i,j,k),gyz(i,j,k),gzz(i,j,k))
@@ -505,9 +505,9 @@
 
   else 
   !$OMP PARALLEL DO PRIVATE(k,j,i,det)
-  do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1
-    do j = GRHydro_stencil,cctk_lsh(2)-GRHydro_stencil+1
-      do i = GRHydro_stencil,cctk_lsh(1)-GRHydro_stencil+1
+  do k = 1,cctk_lsh(3)
+    do j = 1,cctk_lsh(2)
+      do i = 1,cctk_lsh(1)
         
          det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k), \
               gyy(i,j,k),gyz(i,j,k),gzz(i,j,k))
@@ -747,9 +747,9 @@
   CCTK_REAL :: det
   
   !$OMP PARALLEL DO PRIVATE(k,j,i,det)
-  do k = GRHydro_stencil,cctk_lsh(3)-GRHydro_stencil+1
-    do j = GRHydro_stencil,cctk_lsh(2)-GRHydro_stencil+1
-      do i = GRHydro_stencil,cctk_lsh(1)-GRHydro_stencil+1
+  do k = 1,cctk_lsh(3)
+    do j = 1,cctk_lsh(2)
+      do i = 1,cctk_lsh(1)
         
         det = SPATIAL_DETERMINANT(gxx(i,j,k),gxy(i,j,k),gxz(i,j,k),\
              gyy(i,j,k),gyz(i,j,k),gzz(i,j,k))