diff options
author | allen <allen@f5a6acaf-da7d-456b-b0a8-35edbc60b392> | 2002-08-21 15:44:37 +0000 |
---|---|---|
committer | allen <allen@f5a6acaf-da7d-456b-b0a8-35edbc60b392> | 2002-08-21 15:44:37 +0000 |
commit | 5f565bb6944f8896ca8ea80931be2e2b78344589 (patch) | |
tree | a3c13f8afe14ad00fa92916172f544ac640bd2ea | |
parent | edcb8c19f8a53cc647da3d5f876cc41257eb6999 (diff) |
Fixed error setting solution at special case of r=0
Closes Cactus/1085
Also made the gaussian in IDScalarWaveCXX consistent with gaussians in
the other ID thorns, hence new testsuite results.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/IDScalarWave/trunk@59 f5a6acaf-da7d-456b-b0a8-35edbc60b392
-rw-r--r-- | src/InitialData.F77 | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/InitialData.F77 b/src/InitialData.F77 index 8a8da75..08de276 100644 --- a/src/InitialData.F77 +++ b/src/InitialData.F77 @@ -65,7 +65,7 @@ do k=1, cctk_lsh(3) do j=1, cctk_lsh(2) do i=1, cctk_lsh(1) - + xp=x(i,j,k) yp=y(i,j,k) zp=z(i,j,k) @@ -75,7 +75,7 @@ phi(i,j,k) = amplitude*exp(-(rp-radius)**2/sigma**2) if (rp .eq. 0.0) then - phi_p(i,j,k) = (1.0 - 2.0*dt**2/sigma)*exp(-dt**2/sigma) + phi_p(i,j,k) = amplitude*(1.0 - 2.0*dt**2/sigma)*exp(-dt**2/sigma) else phi_p(i,j,k) = amplitude/2.0*(rp-dt)/rp* & exp( - ( (rp-radius-dt)/sigma)**2 ) |