diff options
| author | eschnett <> | 2001-03-19 20:30:00 +0000 |
|---|---|---|
| committer | eschnett <> | 2001-03-19 20:30:00 +0000 |
| commit | fc0e7af205874d55377d18ffdb68555df2843d32 (patch) | |
| tree | ae2b89c76f3a3af299a8699833e391f42dddc05f /CarpetExtra/IDHydroToy | |
| parent | 61b00f25a9820d71401ee6ec3188f6d7b0fa3412 (diff) | |
small
darcs-hash:20010319203034-f6438-ab6a4b4cd4ff91f151ade70f9b3e2a813ee0503b.gz
Diffstat (limited to 'CarpetExtra/IDHydroToy')
| -rw-r--r-- | CarpetExtra/IDHydroToy/src/InitialData.F77 | 24 | ||||
| -rw-r--r-- | CarpetExtra/IDHydroToy/src/make.code.defn | 4 |
2 files changed, 19 insertions, 9 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 diff --git a/CarpetExtra/IDHydroToy/src/make.code.defn b/CarpetExtra/IDHydroToy/src/make.code.defn index d766da22c..3b82fb82f 100644 --- a/CarpetExtra/IDHydroToy/src/make.code.defn +++ b/CarpetExtra/IDHydroToy/src/make.code.defn @@ -1,8 +1,8 @@ # Main make.code.defn file for thorn IDHydroToy -# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDHydroToy/src/make.code.defn,v 1.1 2001/03/18 22:37:10 eschnett Exp $ +# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/IDHydroToy/src/make.code.defn,v 1.2 2001/03/19 21:30:37 eschnett Exp $ # Source files in this directory -SRCS = InitialData.F77 Startup.F77 +SRCS = InitialData.F77 Startup.F77 erf.f gammln.f gammp.f gcf.f gser.f # Subdirectories containing source files SUBDIRS = |