diff options
-rw-r--r-- | src/InitialData.F77 | 19 |
1 files changed, 10 insertions, 9 deletions
diff --git a/src/InitialData.F77 b/src/InitialData.F77 index 1483002..fc90ac3 100644 --- a/src/InitialData.F77 +++ b/src/InitialData.F77 @@ -103,8 +103,9 @@ c Use kx,ky,kz as number of modes in each direction. do j=1,cctk_lsh(2) do i=1,cctk_lsh(1) - phi(i,j,k) = 0.0d0 + phi(i,j,k) = 0.0d0 phi_p(i,j,k) = 0.0d0 + end do end do end do @@ -112,14 +113,14 @@ c Use kx,ky,kz as number of modes in each direction. end if c zero out corners of phi next */ - phi_n(1, 1, 1) = 0.0 - phi_n(1, 1, cctk_lsh(3)) = 0.0 - phi_n(1, cctk_lsh(2), 1) = 0.0 - phi_n(1, cctk_lsh(2), cctk_lsh(3)) = 0.0 - phi_n(cctk_lsh(1), 1, 1) = 0.0 - phi_n(cctk_lsh(1), 1, cctk_lsh(3)) = 0.0 - phi_n(cctk_lsh(1), cctk_lsh(2), 1) = 0.0 - phi_n(cctk_lsh(1), cctk_lsh(2), cctk_lsh(3)) = 0.0 +c phi_n(1, 1, 1) = 0.0 +c phi_n(1, 1, cctk_lsh(3)) = 0.0 +c phi_n(1, cctk_lsh(2), 1) = 0.0 +c phi_n(1, cctk_lsh(2), cctk_lsh(3)) = 0.0 +c phi_n(cctk_lsh(1), 1, 1) = 0.0 +c phi_n(cctk_lsh(1), 1, cctk_lsh(3)) = 0.0 +c phi_n(cctk_lsh(1), cctk_lsh(2), 1) = 0.0 +c phi_n(cctk_lsh(1), cctk_lsh(2), cctk_lsh(3)) = 0.0 return end |