diff options
| author | jthorn <jthorn@e296648e-0e4f-0410-bd07-d597d9acff87> | 2002-09-03 13:49:01 +0000 |
|---|---|---|
| committer | jthorn <jthorn@e296648e-0e4f-0410-bd07-d597d9acff87> | 2002-09-03 13:49:01 +0000 |
| commit | cec6eeade81ca638c437177dcf9c70b4f678591d (patch) | |
| tree | 167910644b390f15309e681693ff29dc1a03571e /src | |
| parent | 1f7b64ac7523129f097eab1f4a1ee299c49c5101 (diff) | |
Fix a bug spotted by Peter Diener:
The param.ccl and documentation.tex define the Kerr spin parameter to be
(dimensionless) J/m^2. But (before this commit) the code used it with
(dimension = mass) J/m semantics. This meant that the spin was set
wrong whenever the BH mass was different from 1.0
This change changes the code to properly match the documentation.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/Exact/trunk@126 e296648e-0e4f-0410-bd07-d597d9acff87
Diffstat (limited to 'src')
| -rw-r--r-- | src/metrics/Kerr_KerrSchild.F | 6 |
1 files changed, 4 insertions, 2 deletions
diff --git a/src/metrics/Kerr_KerrSchild.F b/src/metrics/Kerr_KerrSchild.F index 3561612..d998fb4 100644 --- a/src/metrics/Kerr_KerrSchild.F +++ b/src/metrics/Kerr_KerrSchild.F @@ -40,13 +40,15 @@ C $Header$ cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc -C Get parameters of the exact solution. +C Get parameters of the exact solution, +C and convert from parameter file spin parameter J/m^2 +C to the J/m definition used in the code here. if (firstcall) then boostv = Kerr_KerrSchild__boost_v eps = Kerr_KerrSchild__epsilon m = Kerr_KerrSchild__mass - a = Kerr_KerrSchild__spin + a = m*Kerr_KerrSchild__spin firstcall = .false. end if |