diff options
author | eschnett <eschnett@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-10-28 16:35:34 +0000 |
---|---|---|
committer | eschnett <eschnett@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2013-10-28 16:35:34 +0000 |
commit | 786972179a0774a422afa8bb7031bcdcf404b20c (patch) | |
tree | 26a8c206fc87e5de96a437c8239061998145fc85 | |
parent | e937d71b0b60c2e1d59332d2dfa650b840f06ebc (diff) |
Correct array indices in etaXYZ
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@581 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
-rw-r--r-- | src/GRHydro_HLLE.F90 | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/GRHydro_HLLE.F90 b/src/GRHydro_HLLE.F90 index 5d38062..34455cf 100644 --- a/src/GRHydro_HLLE.F90 +++ b/src/GRHydro_HLLE.F90 @@ -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)))) |