diff options
| -rw-r--r-- | src/InitialData.F77 | 2 |
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 ) |