diff options
author | goodale <goodale@f5a6acaf-da7d-456b-b0a8-35edbc60b392> | 2000-03-29 07:42:10 +0000 |
---|---|---|
committer | goodale <goodale@f5a6acaf-da7d-456b-b0a8-35edbc60b392> | 2000-03-29 07:42:10 +0000 |
commit | 087726ee09df1a0f2418ab4238ef632c857e5bc8 (patch) | |
tree | f0016e41d9b3a0fcc7a23640b92afa8a2893b2ec | |
parent | 5037534af3cb79d0c8699acd93fd5d70ceca3a64 (diff) |
Changed to use timelevels - thanks for putting the rotation in PUGH Ed.
Note that the no-boundary condition tests assumed that phi_new had zero
at its boundaries. This is, or course, no longer true with timelevel
rotation.
Tom
git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/IDScalarWave/trunk@34 f5a6acaf-da7d-456b-b0a8-35edbc60b392
-rw-r--r-- | src/InitialData.F77 | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/src/InitialData.F77 b/src/InitialData.F77 index 2afb52d..4427e98 100644 --- a/src/InitialData.F77 +++ b/src/InitialData.F77 @@ -52,7 +52,7 @@ phi(i,j,k) = amplitude*cos(kx*x(i,j,k)+ky*y(i,j,k) & +kz*z(i,j,k)+omega*cctk_time) - phi_old(i,j,k) = amplitude*cos(kx*x(i,j,k)+ky*y(i,j,k) + phi_p(i,j,k) = amplitude*cos(kx*x(i,j,k)+ky*y(i,j,k) & +kz*z(i,j,k)+omega*(cctk_time-dt)) end do @@ -68,7 +68,7 @@ phi(i,j,k) = amplitude*exp( & -(sqrt(x(i,j,k)**2 & +y(i,j,k)**2+z(i,j,k)**2)-radius)**2/sigma**2) - phi_old(i,j,k) = amplitude*exp( + phi_p(i,j,k) = amplitude*exp( & -(sqrt(x(i,j,k)**2 & +y(i,j,k)**2+z(i,j,k)**2)-radius-dt)**2/sigma**2) @@ -89,7 +89,7 @@ c Use kx,ky,kz as number of modes in each direction. & sin(kz*(z(i,j,k)-0.5)*pi)* & cos(omega*cctk_time*pi) - phi_old(i,j,k)= amplitude*sin(kx*(x(i,j,k)-0.5)*pi)* + phi_p(i,j,k)= amplitude*sin(kx*(x(i,j,k)-0.5)*pi)* & sin(ky*(y(i,j,k)-0.5)*pi)* & sin(kz*(z(i,j,k)-0.5)*pi)* & cos(omega*(cctk_time-dt)*pi) @@ -104,7 +104,7 @@ c Use kx,ky,kz as number of modes in each direction. do i=1,cctk_lsh(1) phi(i,j,k) = 0.0d0 - phi_old(i,j,k) = 0.0d0 + phi_p(i,j,k) = 0.0d0 end do end do end do |