diff options
| author | schnetter <> | 2001-12-14 16:59:00 +0000 |
|---|---|---|
| committer | schnetter <> | 2001-12-14 16:59:00 +0000 |
| commit | 24f4fe527c1d94256702a904c43c96cd86d95c9f (patch) | |
| tree | 267c58fc859249ea44cdbed9ff94114db4b269d1 /CarpetExtra/IDSpaceTimeToy | |
| parent | f78f6ff07a44eb285372f776c1d8de5db9be141b (diff) | |
Added new thorns CarpetReduce and CarpetRegrid to parameter files.
darcs-hash:20011214165954-07bb3-e68f15a15027098cd26981663394a4365aee510d.gz
Diffstat (limited to 'CarpetExtra/IDSpaceTimeToy')
| -rw-r--r-- | CarpetExtra/IDSpaceTimeToy/src/InitialData.F77 | 92 |
1 files changed, 63 insertions, 29 deletions
diff --git a/CarpetExtra/IDSpaceTimeToy/src/InitialData.F77 b/CarpetExtra/IDSpaceTimeToy/src/InitialData.F77 index 4fe16680d..69b1e710f 100644 --- a/CarpetExtra/IDSpaceTimeToy/src/InitialData.F77 +++ b/CarpetExtra/IDSpaceTimeToy/src/InitialData.F77 @@ -1,5 +1,5 @@ c -*-Fortran-*- -c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDSpaceTimeToy/src/InitialData.F77,v 1.5 2001/08/26 13:59:47 schnetter Exp $ +c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDSpaceTimeToy/src/InitialData.F77,v 1.6 2001/12/14 18:00:02 schnetter Exp $ #include "cctk.h" #include "cctk_Parameters.h" @@ -16,6 +16,7 @@ c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDSpaceTimeToy/src/Ini CCTK_REAL pi CCTK_REAL omega CCTK_REAL dt + CCTK_REAL x,y,z, r integer i,j,k pi = 4*atan(1.d0) @@ -30,20 +31,26 @@ c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDSpaceTimeToy/src/Ini do j=1,cctk_lsh(2) do i=1,cctk_lsh(1) + x = cart3d_x(i,j,k) + y = cart3d_y(i,j,k) + z = cart3d_z(i,j,k) + phi(i,j,k) = amplitude - $ * cos((kx*cart3d_x(i,j,k) + ky*cart3d_y(i,j,k) - $ + kz*cart3d_z(i,j,k) + omega*cctk_time) * pi) + $ * cos((kx*x + ky*y + kz*z + omega*cctk_time) * pi) psi(i,j,k) = - amplitude - $ * sin((kx*cart3d_x(i,j,k) + ky*cart3d_y(i,j,k) - $ + kz*cart3d_z(i,j,k) + omega*cctk_time) * pi) + $ * sin((kx*x + ky*y + kz*z + omega*cctk_time) * pi) $ * pi * omega phi_p(i,j,k) = amplitude - $ * cos((kx*cart3d_x(i,j,k) + ky*cart3d_y(i,j,k) - $ + kz*cart3d_z(i,j,k) + omega*(cctk_time-dt)) * pi) + $ * cos((kx*x + ky*y + kz*z + omega*(cctk_time - dt)) * pi) psi_p(i,j,k) = - amplitude - $ * sin((kx*cart3d_x(i,j,k) + ky*cart3d_y(i,j,k) - $ + kz*cart3d_z(i,j,k) + omega*(cctk_time-dt)) * pi) + $ * sin((kx*x + ky*y + kz*z + omega*(cctk_time - dt)) * pi) + $ * pi * omega + + phi_p_p(i,j,k) = amplitude + $ * cos((kx*x + ky*y + kz*z + omega*(cctk_time - 2*dt)) * pi) + psi_p_p(i,j,k) = - amplitude + $ * sin((kx*x + ky*y + kz*z + omega*(cctk_time - 2*dt)) * pi) $ * pi * omega end do @@ -56,15 +63,22 @@ c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDSpaceTimeToy/src/Ini do j=1,cctk_lsh(2) do i=1,cctk_lsh(1) - phi(i,j,k) = amplitude / spher3d_r(i,j,k) - $ * exp(- (spher3d_r(i,j,k) - radius - cctk_time)**2 / sigma**2) + r = spher3d_r(i,j,k) + + phi(i,j,k) = amplitude / r + $ * exp(- (r - radius - cctk_time)**2 / sigma**2) psi(i,j,k) = phi(i,j,k) - $ * 2 * (spher3d_r(i,j,k) - radius - cctk_time) / sigma**2 + $ * 2 * (r - radius - cctk_time) / sigma**2 - phi_p(i,j,k) = amplitude / spher3d_r(i,j,k) - $ * exp(- (spher3d_r(i,j,k) - radius - (cctk_time-dt))**2 / sigma**2) + phi_p(i,j,k) = amplitude / r + $ * exp(- (r - radius - (cctk_time - dt))**2 / sigma**2) psi_p(i,j,k) = phi(i,j,k) - $ * 2 * (spher3d_r(i,j,k) - radius - (cctk_time-dt)) / sigma**2 + $ * 2 * (r - radius - (cctk_time - dt)) / sigma**2 + + phi_p_p(i,j,k) = amplitude / r + $ * exp(- (r - radius - (cctk_time - 2*dt))**2 / sigma**2) + psi_p_p(i,j,k) = phi(i,j,k) + $ * 2 * (r - radius - (cctk_time - 2*dt)) / sigma**2 end do end do @@ -78,28 +92,44 @@ c Use kx,ky,kz as number of modes in each direction. do j=1,cctk_lsh(2) do i=1,cctk_lsh(1) + x = cart3d_x(i,j,k) + y = cart3d_y(i,j,k) + z = cart3d_z(i,j,k) + phi(i,j,k) = amplitude - $ * sin(kx * (cart3d_x(i,j,k) - 0.5d0) * pi) - $ * sin(ky * (cart3d_y(i,j,k) - 0.5d0) * pi) - $ * sin(kz * (cart3d_z(i,j,k) - 0.5d0) * pi) + $ * sin(kx * (x - 0.5d0) * pi) + $ * sin(ky * (y - 0.5d0) * pi) + $ * sin(kz * (z - 0.5d0) * pi) $ * cos(omega * cctk_time * pi) psi(i,j,k) = - amplitude - $ * sin(kx * (cart3d_x(i,j,k) - 0.5d0) * pi) - $ * sin(ky * (cart3d_y(i,j,k) - 0.5d0) * pi) - $ * sin(kz * (cart3d_z(i,j,k) - 0.5d0) * pi) + $ * sin(kx * (x - 0.5d0) * pi) + $ * sin(ky * (y - 0.5d0) * pi) + $ * sin(kz * (z - 0.5d0) * pi) $ * sin(omega * cctk_time * pi) $ * omega * pi phi_p(i,j,k) = amplitude - $ * sin(kx * (cart3d_x(i,j,k) - 0.5d0) * pi) - $ * sin(ky * (cart3d_y(i,j,k) - 0.5d0) * pi) - $ * sin(kz * (cart3d_z(i,j,k) - 0.5d0) * pi) - $ * cos(omega * (cctk_time-dt) * pi) + $ * sin(kx * (x - 0.5d0) * pi) + $ * sin(ky * (y - 0.5d0) * pi) + $ * sin(kz * (z - 0.5d0) * pi) + $ * cos(omega * (cctk_time - dt) * pi) psi_p(i,j,k) = - amplitude - $ * sin(kx * (cart3d_x(i,j,k) - 0.5d0) * pi) - $ * sin(ky * (cart3d_y(i,j,k) - 0.5d0) * pi) - $ * sin(kz * (cart3d_z(i,j,k) - 0.5d0) * pi) - $ * sin(omega * (cctk_time-dt) * pi) + $ * sin(kx * (x - 0.5d0) * pi) + $ * sin(ky * (y - 0.5d0) * pi) + $ * sin(kz * (z - 0.5d0) * pi) + $ * sin(omega * (cctk_time - dt) * pi) + $ * omega * pi + + phi_p_p(i,j,k) = amplitude + $ * sin(kx * (x - 0.5d0) * pi) + $ * sin(ky * (y - 0.5d0) * pi) + $ * sin(kz * (z - 0.5d0) * pi) + $ * cos(omega * (cctk_time - 2*dt) * pi) + psi_p_p(i,j,k) = - amplitude + $ * sin(kx * (x - 0.5d0) * pi) + $ * sin(ky * (y - 0.5d0) * pi) + $ * sin(kz * (z - 0.5d0) * pi) + $ * sin(omega * (cctk_time - 2*dt) * pi) $ * omega * pi end do @@ -118,6 +148,9 @@ c Use kx,ky,kz as number of modes in each direction. phi_p(i,j,k) = 0 psi_p(i,j,k) = 0 + phi_p_p(i,j,k) = 0 + psi_p_p(i,j,k) = 0 + end do end do end do @@ -132,6 +165,7 @@ c Use kx,ky,kz as number of modes in each direction. psi(i,j,k) = psi(i,j,k) - u(i,j,k) psi_p(i,j,k) = psi_p(i,j,k) - u_p(i,j,k) + psi_p_p(i,j,k) = psi_p_p(i,j,k) - u_p_p(i,j,k) end do end do |