diff options
| author | eschnett <> | 2001-03-17 15:04:00 +0000 |
|---|---|---|
| committer | eschnett <> | 2001-03-17 15:04:00 +0000 |
| commit | b0432e3310d34454b72e6b674e249e64561c7e1c (patch) | |
| tree | ae7932c903957d805919c8bfea7877f0116573b9 /CarpetExtra/IDScalarWave | |
| parent | 91a1b9a72e06ca12338fa2b50065052ac84eccec (diff) | |
Removed unnecessary parameter checking. Small changes.
darcs-hash:20010317150459-f6438-65b8ca470951c76852d25e999003fe363358e54e.gz
Diffstat (limited to 'CarpetExtra/IDScalarWave')
| -rw-r--r-- | CarpetExtra/IDScalarWave/src/InitialData.F77 | 8 |
1 files changed, 3 insertions, 5 deletions
diff --git a/CarpetExtra/IDScalarWave/src/InitialData.F77 b/CarpetExtra/IDScalarWave/src/InitialData.F77 index 63fadc2b7..0bb3d8082 100644 --- a/CarpetExtra/IDScalarWave/src/InitialData.F77 +++ b/CarpetExtra/IDScalarWave/src/InitialData.F77 @@ -72,11 +72,9 @@ c call CCTK_INFO ("IDScalarWave_InitialData") do i=1, cctk_lsh(1) phi(i,j,k) = amplitude - $ * exp(- (sqrt(cart3d_x(i,j,k)**2 + cart3d_y(i,j,k)**2 - $ + cart3d_z(i,j,k)**2) - radius)**2 / sigma**2) + $ * exp(- (spher3d_r(i,j,k) - radius)**2 / sigma**2) phi_p(i,j,k) = amplitude - $ * exp(- (sqrt(cart3d_x(i,j,k)**2 + cart3d_y(i,j,k)**2 - $ + cart3d_z(i,j,k)**2) - radius - dt)**2 / sigma**2) + $ * exp(- (spher3d_r(i,j,k)**2 - radius - dt)**2 / sigma**2) end do end do @@ -96,7 +94,7 @@ c Use kx,ky,kz as number of modes in each direction. $ * sin(kz * (cart3d_z(i,j,k) - 0.5d0) * pi) $ * cos(omega * cctk_time * pi) - phi_p(i,j,k)= amplitude + 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) |