diff options
Diffstat (limited to 'src')
-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 |