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

schnetter at cct.lsu.edu schnetter at cct.lsu.edu
Mon Oct 28 11:35:35 CDT 2013


User: eschnett
Date: 2013/10/28 11:35 AM

Modified:
 /trunk/src/
  GRHydro_HLLE.F90

Log:
 Correct array indices in etaXYZ

File Changes:

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

File [modified]: GRHydro_HLLE.F90
Delta lines: +3 -3
===================================================================
--- trunk/src/GRHydro_HLLE.F90	2013-10-27 04:16:35 UTC (rev 580)
+++ trunk/src/GRHydro_HLLE.F90	2013-10-28 16:35:34 UTC (rev 581)
@@ -33,15 +33,15 @@
 
 ! eta across face in x-direction
 #define etaX(i,j,k, v, c) \
-  (0.5d0*(abs(v(i+1,j,k,0)-v(i,j,k,0)) + abs(c(i+1,j,k)-c(i,j,k))))
+  (0.5d0*(abs(v((i)+1,j,k,1)-v(i,j,k,1)) + abs(c((i)+1,j,k)-c(i,j,k))))
 
 ! eta across face in x-direction
 #define etaY(i,j,k, v, c) \
-  (0.5d0*(abs(v(i,j+1,k,1)-v(i,j,k,1)) + abs(c(i,j+1,k)-c(i,j,k))))
+  (0.5d0*(abs(v(i,(j)+1,k,2)-v(i,j,k,2)) + abs(c(i,(j)+1,k)-c(i,j,k))))
   
 ! eta across face in x-direction
 #define etaZ(i,j,k, v, c) \
-  (0.5d0*(abs(v(i,j,k+1,2)-v(i,j,k,2)) + abs(c(i,j,k+1)-c(i,j,k))))
+  (0.5d0*(abs(v(i,j,(k)+1,3)-v(i,j,k,3)) + abs(c(i,j,(k)+1)-c(i,j,k))))
 
   
 



More information about the Commits mailing list