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))))