aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorschnetter <schnetter@f5a6acaf-da7d-456b-b0a8-35edbc60b392>2009-11-05 04:32:28 +0000
committerschnetter <schnetter@f5a6acaf-da7d-456b-b0a8-35edbc60b392>2009-11-05 04:32:28 +0000
commit9503acfa1b370f7fe39b4a41fd9452dd014e75a1 (patch)
tree55afe877802ce373743f2d0338eedd2da67f868e
parent7122ac36741da4f446ac0198a300ebfbc5df3afb (diff)
Correct Gaussian initial data for R=0HEADsvnmaster
git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/IDScalarWave/trunk@75 f5a6acaf-da7d-456b-b0a8-35edbc60b392
-rw-r--r--src/InitialData.F772
1 files changed, 1 insertions, 1 deletions
diff --git a/src/InitialData.F77 b/src/InitialData.F77
index 7cbf7a0..249ee90 100644
--- a/src/InitialData.F77
+++ b/src/InitialData.F77
@@ -78,7 +78,7 @@
phi(i,j,k) = amplitude*exp(-(rp-radius)**2/sigma**2)
if (rp .eq. 0.0) then
- phi_p(i,j,k) = amplitude*(1.0 - 2.0*dt**2/sigma)*exp(-dt**2/sigma)
+ phi_p(i,j,k) = amplitude*(1.0 - 2.0*dt**2/sigma**2)*exp(-dt**2/sigma**2)
else
phi_p(i,j,k) = amplitude/2.0*(rp-dt)/rp*
& exp( - ( (rp-radius-dt)/sigma)**2 )