diff options
| -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 ) |