diff options
Diffstat (limited to 'src/InitialData.F77')
| -rw-r--r-- | src/InitialData.F77 | 18 |
1 files changed, 8 insertions, 10 deletions
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 |