diff options
author | schnetter <schnetter@e296648e-0e4f-0410-bd07-d597d9acff87> | 2004-11-01 21:50:43 +0000 |
---|---|---|
committer | schnetter <schnetter@e296648e-0e4f-0410-bd07-d597d9acff87> | 2004-11-01 21:50:43 +0000 |
commit | f5817d0d77263575f2e3bd59dcc8f7f40c94ace2 (patch) | |
tree | f7bb2939ba530db021a8d00bb5ed4ea8fc5576c8 /src/metrics | |
parent | b76eaee65aee5a6309a84f91d3a8a3424f8fbbf4 (diff) |
Add "Schwarzschild/BL" initial data, which are a Brill-Lindquist black
hole.
Allow "Schwarzschild/KS" initial data to be offset from the origin.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/Exact/trunk@212 e296648e-0e4f-0410-bd07-d597d9acff87
Diffstat (limited to 'src/metrics')
-rw-r--r-- | src/metrics/Kerr_KerrSchild.F77 | 8 | ||||
-rw-r--r-- | src/metrics/Schwarzschild_BL.F77 | 71 | ||||
-rw-r--r-- | src/metrics/make.code.defn | 1 |
3 files changed, 76 insertions, 4 deletions
diff --git a/src/metrics/Kerr_KerrSchild.F77 b/src/metrics/Kerr_KerrSchild.F77 index 79619cd..f2fb45a 100644 --- a/src/metrics/Kerr_KerrSchild.F77 +++ b/src/metrics/Kerr_KerrSchild.F77 @@ -66,10 +66,10 @@ C t0 is never used, but is here for illustration, and we introduce C x0 and y0 also only for clarity. C Note that z0 = 0 means z = vt for the BH. - t0 = gamma * (t - boostv * z) - z0 = gamma * (z - boostv * t) - x0 = x - y0 = y + t0 = gamma * ((t - Kerr_KerrSchild__t) - boostv * (z - Kerr_KerrSchild__z)) + z0 = gamma * ((z - Kerr_KerrSchild__z) - boostv * (t - Kerr_KerrSchild__t)) + x0 = x - Kerr_KerrSchild__x + y0 = y - Kerr_KerrSchild__y C Coordinate distance to center of black hole. Note it moves! diff --git a/src/metrics/Schwarzschild_BL.F77 b/src/metrics/Schwarzschild_BL.F77 new file mode 100644 index 0000000..f9d63f7 --- /dev/null +++ b/src/metrics/Schwarzschild_BL.F77 @@ -0,0 +1,71 @@ +c Schwarzschild spacetime in Brill-Lindquist coordinates. +C +c $Header$ + +#include "cctk.h" +#include "cctk_Parameters.h" + + subroutine Exact__Schwarzschild_BL( + $ x, y, z, t, + $ gdtt, gdtx, gdty, gdtz, + $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, + $ gutt, gutx, guty, gutz, + $ guxx, guyy, guzz, guxy, guyz, guzx, + $ psi, Tmunu_flag) + + implicit none + + DECLARE_CCTK_PARAMETERS + +c input arguments + CCTK_REAL x, y, z, t + +c output arguments + CCTK_REAL gdtt, gdtx, gdty, gdtz, + $ gdxx, gdyy, gdzz, gdxy, gdyz, gdzx, + $ gutt, gutx, guty, gutz, + $ guxx, guyy, guzz, guxy, guyz, guzx + CCTK_REAL psi + LOGICAL Tmunu_flag + +c local variables + CCTK_REAL eps, m + +c local variables + CCTK_REAL r, psi4 + +C This is a vacuum spacetime with no cosmological constant + Tmunu_flag = .false. + +C Get parameters of the exact solution. + + eps = Schwarzschild_BL__epsilon + m = Schwarzschild_BL__mass + + r = ((x**2 + y**2 + z**2)**2 + eps**4) ** 0.25d0 + psi4 = (1 + m / (2 * r)) ** 4 + + gdtt = -1 + gdtx = 0 + gdty = 0 + gdtz = 0 + gdxx = psi4 + gdyy = psi4 + gdzz = psi4 + gdxy = 0 + gdyz = 0 + gdzx = 0 + + gutt = -1 + gutx = 0 + guty = 0 + gutz = 0 + guxx = 1 / psi4 + guyy = 1 / psi4 + guzz = 1 / psi4 + guxy = 0 + guyz = 0 + guzx = 0 + + return + end diff --git a/src/metrics/make.code.defn b/src/metrics/make.code.defn index b50ab0f..ba33ccd 100644 --- a/src/metrics/make.code.defn +++ b/src/metrics/make.code.defn @@ -14,6 +14,7 @@ SRCS = Minkowski.F77 \ \ Schwarzschild_EF.F77 \ Schwarzschild_PG.F77 \ + Schwarzschild_BL.F77 \ Schwarzschild_Novikov.F77 \ Kerr_BoyerLindquist.F77 \ Kerr_KerrSchild.F77 \ |