diff options
Diffstat (limited to 'CarpetExtra/WaveToyF77')
-rw-r--r-- | CarpetExtra/WaveToyF77/src/WaveToy.F77 | 28 |
1 files changed, 14 insertions, 14 deletions
diff --git a/CarpetExtra/WaveToyF77/src/WaveToy.F77 b/CarpetExtra/WaveToyF77/src/WaveToy.F77 index adcdea62a..bbdf70457 100644 --- a/CarpetExtra/WaveToyF77/src/WaveToy.F77 +++ b/CarpetExtra/WaveToyF77/src/WaveToy.F77 @@ -65,13 +65,13 @@ c ----------------- dy2i = 1/dy2 dz2i = 1/dz2 - istart = 2 - jstart = 2 - kstart = 2 + istart = 1+cctk_nghostzones(1) + jstart = 1+cctk_nghostzones(2) + kstart = 1+cctk_nghostzones(3) - iend = cctk_lsh(1)-1 - jend = cctk_lsh(2)-1 - kend = cctk_lsh(3)-1 + iend = cctk_lsh(1)-cctk_nghostzones(1) + jend = cctk_lsh(2)-cctk_nghostzones(2) + kend = cctk_lsh(3)-cctk_nghostzones(3) factor = 2 * (1 - dt2 * (dx2i + dy2i + dz2i)) @@ -81,11 +81,11 @@ c ---------------- 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 @@ -132,9 +132,9 @@ c call CCTK_INFO ("WaveToyF77_Boundaries") c Set the stencil width c --------------------- - sw(1) = 1 - sw(2) = 1 - sw(3) = 1 + sw(1) = cctk_nghostzones(1) + sw(2) = cctk_nghostzones(2) + sw(3) = cctk_nghostzones(3) c Apply the outer boundary conditions c ----------------------------------- |