correct computation of shift vector in Kerr metric, along with the documentation and a proof that this is right

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAnalyticBH/trunk@175 6a3ddf76-46e1-4315-99d9-bc56cac1ef84

1 files changed, 1 insertions, 1 deletions

diff --git a/src/Kerr.c b/src/Kerr.c
index ae39fd6..d7966ca 100644
--- a/src/Kerr.c
+++ b/src/Kerr.c
@@ -132,7 +132,7 @@ void KerrID(CCTK_ARGUMENTS)
     p2=a_2+r_2-a*beta_phi;
     /* drdR=sqrt_Delta/R; */
     lapse=sqrt_Delta/sqrt(p2);
-    shift_phi=-2*m*rK*a/p2;
+    shift_phi=-2*m*rK*a/(p2*Sigma);
     
     /* Kerr metric in quasi-isotropic coordinates */
     /* ds^2=Phi4*(dR^2+R^2*dth^2+R^2*Chi^2*sin(th)^2*dphi^2) */