[Commits] [svn:einsteintoolkit] GRHydro/ (Rev. 514)

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/05/13 12:30 AM

Modified:
 /branches/ET_2013_05/
 /branches/ET_2013_05/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>

Directory Changes:

Directory: /branches/ET_2013_05/
================================

Property changes on: branches/ET_2013_05
___________________________________________________________________
Directory: /svn:mergeinfo/
==========================

   - /branches/divclean:231-243
/branches/hot_and_MHD_temp_dev:166-185
/trunk:510
   + /branches/divclean:231-243
/branches/hot_and_MHD_temp_dev:166-185
/trunk:510,512

File Changes:

Directory: /branches/ET_2013_05/src/
====================================

File [modified]: GRHydro_Prim2ConM.F90
Delta lines: +9 -9
===================================================================
--- branches/ET_2013_05/src/GRHydro_Prim2ConM.F90	2013-05-13 05:27:57 UTC (rev 513)
+++ branches/ET_2013_05/src/GRHydro_Prim2ConM.F90	2013-05-13 05:29:59 UTC (rev 514)
@@ -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))