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

[rhaas at tapir.caltech.edu]

User: rhaas
Date: 2013/08/13 09:56 AM

Modified:
 /trunk/src/
  GRHydro_Prim2ConM.F90

Log:
 Adjust loops for Prim2Con call (needs to run in the ghost zones as well after being called from ZelmaniLeak).
 
 From: Philipp Moesta <pmoesta at tapir.caltech.edu>

File Changes:

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

File [modified]: GRHydro_Prim2ConM.F90
Delta lines: +4 -4
===================================================================
--- trunk/src/GRHydro_Prim2ConM.F90	2013-08-13 14:56:43 UTC (rev 578)
+++ trunk/src/GRHydro_Prim2ConM.F90	2013-08-13 14:56:46 UTC (rev 579)
@@ -585,10 +585,10 @@
 
   if(evolve_Y_e.ne.0) then
      !$OMP PARALLEL DO PRIVATE(i, j, k)
-     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
-              Y_e_con(i,j,k) = Y_e(i,j,k) * dens(i,j,k)
+    do k = 1,cctk_lsh(3)
+      do j = 1,cctk_lsh(2)
+        do i = 1,cctk_lsh(1)
+                Y_e_con(i,j,k) = Y_e(i,j,k) * dens(i,j,k)
            enddo
         enddo
      enddo