diff options
| -rw-r--r-- | doc/documentation.tex | 7 | ||||
| -rw-r--r-- | src/InitialData.F77 | 18 |
2 files changed, 14 insertions, 11 deletions
diff --git a/doc/documentation.tex b/doc/documentation.tex index 3a63740..5cf9d16 100644 --- a/doc/documentation.tex +++ b/doc/documentation.tex @@ -77,7 +77,12 @@ The initial profile is \begin{equation} \phi(r) = A \exp (- r^2/\sigma) \end{equation} -and is set with the parameters +with the solution at the origin being +\begin{equation} +\Psi(r=0,t) = \left(1-2\frac{t^2}{\sigma}\right)\exp(-t^2/\sigma) +\end{equation} + +The Gaussian solution is set with the parameters \begin{itemize} \item {\tt amplitude} = $A$ diff --git a/src/InitialData.F77 b/src/InitialData.F77 index d3af842..9b56815 100644 --- a/src/InitialData.F77 +++ b/src/InitialData.F77 @@ -73,17 +73,15 @@ rp = sqrt(xp*xp+yp*yp+zp*zp) phi(i,j,k) = amplitude*exp(-(rp-radius)**2/sigma**2) - phi_p(i,j,k) = amplitude/2.0*(rp-dt)/rp* - & exp( - ( (rp-radius-dt)/sigma)**2 ) - & + amplitude/2.0*(rp+dt)/rp* - & exp( - ( (rp-radius+dt)/sigma)**2 ) - -c phi_p[index] = amplitude/2.0*(R-dt)/R* -c exp( - sqr( (R - radius - dt)/ sigma ) ) -c + amplitude/2.0*(R+dt)/R* -c exp( - sqr( (R - radius + dt)/ sigma ) ) -c ; + if (rp .eq. 0.0) then + phi_p(i,j,k) = (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 ) + & + amplitude/2.0*(rp+dt)/rp* + & exp( - ( (rp-radius+dt)/sigma)**2 ) + endif end do end do |