From 087726ee09df1a0f2418ab4238ef632c857e5bc8 Mon Sep 17 00:00:00 2001 From: goodale Date: Wed, 29 Mar 2000 07:42:10 +0000 Subject: 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 --- src/InitialData.F77 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) (limited to 'src') 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 -- cgit v1.2.3