diff options
Diffstat (limited to 'src/GRHydro_WENOReconstruct.F90')
-rw-r--r-- | src/GRHydro_WENOReconstruct.F90 | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/src/GRHydro_WENOReconstruct.F90 b/src/GRHydro_WENOReconstruct.F90 index fb03f66..9c40d95 100644 --- a/src/GRHydro_WENOReconstruct.F90 +++ b/src/GRHydro_WENOReconstruct.F90 @@ -259,12 +259,12 @@ subroutine GRHydro_WENOReconstruct1d(order, nx, v, vminus, vplus, trivial_rp, & !!$ Calculate the reconstruction do j = 1, 5 - vplus(i) = vplus(i) + wplus1 * weno_coeffs(1,j)*v(i-3+j) & - + wplus2 * weno_coeffs(2,j)*v(i-3+j) & - + wplus3 * weno_coeffs(3,j)*v(i-3+j) - vminus(i) = vminus(i) + wminus1 * weno_coeffs(3,6-j)*v(i-3+j) & - + wminus2 * weno_coeffs(2,6-j)*v(i-3+j) & - + wminus3 * weno_coeffs(1,6-j)*v(i-3+j) + vplus(i) = vplus(i) + wplus1 * weno_coeffs(1,j)*v(i-3+j) & + + wplus2 * weno_coeffs(2,j)*v(i-3+j) & + + wplus3 * weno_coeffs(3,j)*v(i-3+j) + vminus(i) = vminus(i) + wminus1 * weno_coeffs(3,6-j)*v(i-3+j) & + + wminus2 * weno_coeffs(2,6-j)*v(i-3+j) & + + wminus3 * weno_coeffs(1,6-j)*v(i-3+j) end do !vminus(i) = v(i) !vplus(i) = v(i) |