Analytic data for black hole spacetimes

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

6 files changed, 166 insertions, 0 deletions

diff --git a/README b/README
new file mode 100644
index 0000000..df60551
--- /dev/null
+++ b/README
@@ -0,0 +1,20 @@
+Cactus Code Thorn IDAnalyticBH
+Authors    : ...
+Managed by : ... <...@...........>
+Version    : ...  
+CVS info   : $Header$
+--------------------------------------------------------------------------
+
+1. Purpose of the thorn
+
+This thorn does ...
+
+2. Dependencies of the thorn
+
+This thorn additionally requires thorns ...
+
+3. Thorn distribution
+
+This thorn is available to ...
+
+4. Additional information
diff --git a/interface.ccl b/interface.ccl
new file mode 100644
index 0000000..fdfbf53
--- /dev/null
+++ b/interface.ccl
@@ -0,0 +1,6 @@
+# Interface definition for thorn IDAnalyticBH
+# $Header$
+
+implements: idanalyticbh
+inherits: einstein
+
diff --git a/param.ccl b/param.ccl
new file mode 100644
index 0000000..574d28e
--- /dev/null
+++ b/param.ccl
@@ -0,0 +1,35 @@
+# Parameter definitions for thorn IDAnalyticBH
+# $Header$
+
+private:
+
+REAL mass "Mass of black hole"
+{
+: :: "Not sure if it can be negative or no"
+} 2.0
+
+REAL mu "Misner mu value"
+{
+: :: "Don't know allowed range"
+} 1.2
+
+INTEGER nmax "Number of terms to include for Misner series"
+{
+0: 
+} 30
+
+#friend:einstein
+
+#EXTENDS LOGICAL use_conformal ""
+#{ 
+#}
+
+#EXTENDS KEYWORD model ""
+#{
+#  "schwarzschild" :: "One Schwarzshild black hole"
+#} 
+
+#EXTENDS KEYWORD lapseinit ""
+#{
+#  "schwarz" :: "Set lapse to Schwarzschild"
+#}
diff --git a/schedule.ccl b/schedule.ccl
new file mode 100644
index 0000000..bc04f17
--- /dev/null
+++ b/schedule.ccl
@@ -0,0 +1,12 @@
+# Schedule definitions for thorn IDAnalyticBH
+# $Header$
+
+if (CCTK_Equals(model,"schwarzschild")) 
+{ 
+   schedule Schwarzschild at CACTUS_INITIAL
+   {
+     LANG: Fortran
+   } "Construct initial data for a single Schwarzschild black hole"
+
+}
+ 
\ No newline at end of file
diff --git a/src/Schwarzschild.F b/src/Schwarzschild.F
new file mode 100644
index 0000000..2035d5e
--- /dev/null
+++ b/src/Schwarzschild.F
@@ -0,0 +1,84 @@
+#include "cctk.h"
+#include "declare_arguments.h"
+#include "declare_parameters.h"
+
+      subroutine Schwarzschild(CCTK_FARGUMENTS)
+
+      implicit none
+
+      DECLARE_CCTK_FARGUMENTS
+      DECLARE_PARAMETERS
+
+      REAL zero,one,two,three
+
+      REAL A,B,C,t1,r0
+
+      integer CCTK_Equals
+      integer i,j,k
+      integer is,js,ks,ie,je,ke
+      REAL inval
+
+      zero = 0.d0
+      one = 1.d0
+      two = 2.d0
+      three = 3.d0
+      
+      print*, 'One Black Hole:  Throat at :',mass/two
+      
+c     conformal metric flag
+      if (use_conformal==1) then
+
+         print*,'Evolving conformal metric (psi derivs in the sources)'
+
+c        Compute conformal factor
+         psi = ( one + mass/two/r)
+         
+c        derivatives of psi / psi
+         psix = -x*mass/two/r**3 / psi
+         psiy = -y*mass/two/r**3 / psi
+         psiz = -z*mass/two/r**3 / psi
+         
+         psixy = three*x*y*mass/two/r**5 /psi
+         psixz = three*x*z*mass/two/r**5 /psi
+         psiyz = three*y*z*mass/two/r**5 /psi
+         
+         psixx  = (three*x**2 - r**2)*mass/two/r**5/psi
+         psiyy  = (three*y**2 - r**2)*mass/two/r**5/psi
+         psizz  = (three*z**2 - r**2)*mass/two/r**5/psi          
+         
+         gxx = one
+         gyy = one
+         gzz = one
+         gxy = zero 
+         gxz = zero  
+         gyz = zero   
+          
+      else
+
+         print*,'Evolving ***non-conformal*** metric'         
+
+         gxx = ( one + mass/two/r)**4
+         gyy = gxx
+         gzz = gxx
+         gxy = zero
+         gxz = zero
+         gyz = zero
+         
+      endif
+
+c     If the initial lapse is not one ...
+      if (CCTK_Equals(lapseinit,"schwarz")==1) then 
+         print *,"Initialise with schwarzschild lapse"
+         alp = (2.*r - mass)/(2.*r+mass)
+      endif
+      
+c     time symmetric initial slice
+      hxx = zero 
+      hxy = zero
+      hxz = zero
+      hyy = zero
+      hyz = zero
+      hzz = zero
+
+      return
+      end subroutine Schwarzschild
diff --git a/src/make.code.defn b/src/make.code.defn
new file mode 100644
index 0000000..d67fde6
--- /dev/null
+++ b/src/make.code.defn
@@ -0,0 +1,9 @@
+# Main make.code.defn file for thorn IDAnalyticBH
+# $Header$
+
+# Source files in this directory
+SRCS = Schwarzschild.F
+
+# Subdirectories containing source files
+SUBDIRS = 
+