diff options
Diffstat (limited to 'src/GRHydro_MP5Reconstruct.F90')
| -rw-r--r-- | src/GRHydro_MP5Reconstruct.F90 | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/GRHydro_MP5Reconstruct.F90 b/src/GRHydro_MP5Reconstruct.F90 index be42731..72b2f73 100644 --- a/src/GRHydro_MP5Reconstruct.F90 +++ b/src/GRHydro_MP5Reconstruct.F90 @@ -40,7 +40,7 @@ subroutine GRHydro_MP5Reconstruct1d(nx, v, vminus, vplus, trivial_rp, & DECLARE_CCTK_PARAMETERS - CCTK_INT :: nx, i, j, k, r + CCTK_INT :: nx, i CCTK_REAL, dimension(nx) :: v, vplus, vminus CCTK_INT, dimension(nx) :: hydro_excision_mask @@ -95,7 +95,7 @@ subroutine GRHydro_MP5Reconstruct1d(nx, v, vminus, vplus, trivial_rp, & #define MINMOD4(w,x,y,z) \ 0.125d0*( sign(one,w)+sign(one,x) )*abs( (sign(one,w)+sign(one,y)) * (sign(one,w)+sign(one,z)) )*min(abs(w), abs(x), abs(y), abs(z)) -#define MP5(am2, am1, a, ap1, ap2, arecon) \ +#define MP5(am2, am1, a, ap1, ap2, arecon) \ vl = (2.0d0*am2 - 13.0d0*am1 + 47.0d0*a + 27.0d0*ap1 - 3.0d0*ap2)/60.0d0 &&\ vmp = a + MINMOD( ap1-a, mp5_alpha*(a-am1) ) &&\ if ((vl-a)*(vl-vmp) .le. mp5_eps) then &&\ |