diff options
| -rw-r--r-- | src/WaveToy.F90 | 10 |
1 files changed, 5 insertions, 5 deletions
diff --git a/src/WaveToy.F90 b/src/WaveToy.F90 index 780b944..915df8a 100644 --- a/src/WaveToy.F90 +++ b/src/WaveToy.F90 @@ -76,11 +76,11 @@ subroutine WaveToyFreeF90_Evolution(CCTK_ARGUMENTS) do j = jstart, jend do i = istart, iend - phi_n(i,j,k) = factor*phi(i,j,k) - & - phi_p(i,j,k) + (dt2) * & - ((phi(i+1,j,k)+phi(i-1,j,k))*dx2i & - +(phi(i,j+1,k)+phi(i,j-1,k))*dy2i & - +(phi(i,j,k+1)+phi(i,j,k-1))*dz2i) + phi(i,j,k) = factor*phi_p(i,j,k) - & + phi_p_p(i,j,k) + (dt2) * & + ((phi_p(i+1,j,k)+phi_p(i-1,j,k))*dx2i & + +(phi_p(i,j+1,k)+phi_p(i,j-1,k))*dy2i & + +(phi_p(i,j,k+1)+phi_p(i,j,k-1))*dz2i) end do end do |