From 9503acfa1b370f7fe39b4a41fd9452dd014e75a1 Mon Sep 17 00:00:00 2001 From: schnetter Date: Thu, 5 Nov 2009 04:32:28 +0000 Subject: Correct Gaussian initial data for R=0 git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/IDScalarWave/trunk@75 f5a6acaf-da7d-456b-b0a8-35edbc60b392 --- src/InitialData.F77 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) 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 ) -- cgit v1.2.3