diff options
Diffstat (limited to 'CarpetExtra/IDHydroToy/src/InitialData.F77')
-rw-r--r-- | CarpetExtra/IDHydroToy/src/InitialData.F77 | 24 |
1 files changed, 17 insertions, 7 deletions
diff --git a/CarpetExtra/IDHydroToy/src/InitialData.F77 b/CarpetExtra/IDHydroToy/src/InitialData.F77 index 7c03722ad..5bbb98032 100644 --- a/CarpetExtra/IDHydroToy/src/InitialData.F77 +++ b/CarpetExtra/IDHydroToy/src/InitialData.F77 @@ -1,5 +1,5 @@ c -*-Fortran-*- -c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDHydroToy/src/InitialData.F77,v 1.1 2001/03/18 22:37:10 eschnett Exp $ +c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDHydroToy/src/InitialData.F77,v 1.2 2001/03/19 21:30:37 eschnett Exp $ #include "cctk.h" #include "cctk_Parameters.h" @@ -17,6 +17,11 @@ c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDHydroToy/src/Initial CCTK_REAL omega integer i,j,k + CCTK_REAL vr + + external erf + real*8 erf + pi = 4*atan(1.d0) omega = sqrt(kx**2+ky**2+kz**2) @@ -44,12 +49,17 @@ c $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDHydroToy/src/Initial do j=1,cctk_lsh(2) do i=1,cctk_lsh(1) - u(i,j,k) = amplitude - $ * exp(- (spher3d_r(i,j,k) - radius + cctk_time)**2 - $ / sigma**2) - vx(i,j,k) = u(i,j,k) * cart3d_x(i,j,k)/spher3d_r(i,j,k) - vy(i,j,k) = u(i,j,k) * cart3d_y(i,j,k)/spher3d_r(i,j,k) - vz(i,j,k) = u(i,j,k) * cart3d_z(i,j,k)/spher3d_r(i,j,k) + u(i,j,k) = amplitude / spher3d_r(i,j,k) + $ * exp(- (spher3d_r(i,j,k) - radius + cctk_time)**2 / sigma**2) + +c vr = -1/2*A*(sqrt(Pi)*sigma*erf((t-r+R)/sigma)+2*exp(-(t-r+R)^2/(sigma^2))*r)/(r^2) + vr = -0.5d0 * amplitude + $ * (sqrt(pi) * sigma * erf((cctk_time - spher3d_r(i,j,k) + radius) / sigma) + $ + 2*exp(-(cctk_time - spher3d_r(i,j,k) + radius)**2 / sigma**2) * spher3d_r(i,j,k)) + $ / spher3d_r(i,j,k)**2 + vx(i,j,k) = vr * cart3d_x(i,j,k)/spher3d_r(i,j,k) + vy(i,j,k) = vr * cart3d_y(i,j,k)/spher3d_r(i,j,k) + vz(i,j,k) = vr * cart3d_z(i,j,k)/spher3d_r(i,j,k) end do end do |