*** empty log message ***

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/RotatingDBHIVP/trunk@30 535fb057-194f-0410-b5a5-c63992f15602

2 files changed, 161 insertions, 0 deletions

diff --git a/src/ParamCheck.c b/src/ParamCheck.c
new file mode 100644
index 0000000..1a130b8
--- /dev/null
+++ b/src/ParamCheck.c
@@ -0,0 +1,107 @@
+ /*@@
+   @file      ParamCheck.c
+   @date      Fri Apr 26 18:03:09 2002
+   @author    Tom Goodale
+   @desc 
+              Check the parameters for RotatingDBHIVP
+   @enddesc 
+   @version $Header: /arrangements/PerturbedBH2/RotatingDBHIVP/src/ParamCheck.c,v 1.1 2002/04/30 guzman
+ @@*/
+
+#include "cctk.h"
+
+#include "cctk_Arguments.h"
+#include "cctk_Parameters.h"
+
+static const char *rcsid = "$Header$";
+
+CCTK_FILEVERSION(RotatingDBHIVP_ParamCheck_c)
+
+/********************************************************************
+ *********************     Local Data Types   ***********************
+ ********************************************************************/
+
+/********************************************************************
+ ********************* Local Routine Prototypes *********************
+ ********************************************************************/
+
+/********************************************************************
+ ***************** Scheduled Routine Prototypes *********************
+ ********************************************************************/
+
+void RotatingDBHIVP_ParamCheck(CCTK_ARGUMENTS);
+void RotatingDBHIVP_ConformalCheck(CCTK_ARGUMENTS);
+
+/********************************************************************
+ ********************* Other Routine Prototypes *********************
+ ********************************************************************/
+
+/********************************************************************
+ *********************     Local Data   *****************************
+ ********************************************************************/
+
+/********************************************************************
+ *********************     External Routines   **********************
+ ********************************************************************/
+
+void RotatingDBHIVP_ParamCheck(CCTK_ARGUMENTS)
+{
+  DECLARE_CCTK_ARGUMENTS
+  DECLARE_CCTK_PARAMETERS
+
+  /* Do we know how to deal with this type of metric ? */
+  if(! CCTK_EQUALS(metric_type, "physical") &&
+     ! CCTK_EQUALS(metric_type, "static conformal"))
+  {
+    CCTK_PARAMWARN("Unknown ADMBase::metric_type - known types are \"physical\" and \"static conformal\"");
+  }
+
+  if(CCTK_EQUALS(metric_type, "static conformal") &&
+     ! CCTK_EQUALS(conformal_storage, "factor+derivs+2nd derivs"))
+  {
+    CCTK_PARAMWARN("ADMConstraints can currently only work with a physical metric or a static conf");
+  }
+   
+}
+
+/********************************************************************
+ *********************     Local Routines   *************************
+ ********************************************************************/
+
+ /*@@
+   @routine    RotatingDBHIVP_ConformalCheck
+   @date       Tu Apr 30 2002
+   @author     Goodale-Guzman
+   @desc
+   Check that the initial data has setup enough derivatives.
+   @enddesc
+   @calls
+   @calledby
+   @history
+
+   @endhistory
+
+ @@*/
+void RotatingDBHIVP_ConformalCheck(CCTK_ARGUMENTS)
+{
+  DECLARE_CCTK_ARGUMENTS;
+  DECLARE_CCTK_PARAMETERS;
+
+  if(CCTK_EQUALS(metric_type, "physical"))
+  {
+    /* Use this internally so we don't need to check the parameter again. */
+    *conformal_state = 0;
+  }
+  else if(CCTK_EQUALS(metric_type, "static conformal"))
+  {
+    if(*conformal_state < 3)
+    {
+      CCTK_WARN(0, "RotatingDBHIVP needs second derivatives of the conformal factor when running with a static conformal metric");
+    }
+  }
+
+}
+/********************************************************************
+ *********************     Local Routines   *************************
+ ********************************************************************/
+
diff --git a/test/xeon_test_rdbh.par b/test/xeon_test_rdbh.par
new file mode 100755
index 0000000..a4c25a5
--- /dev/null
+++ b/test/xeon_test_rdbh.par
@@ -0,0 +1,54 @@
+######################################################################
+#!DESC "Initial data for axisymmetric distored black hole"
+######################################################################
+
+ActiveThorns = "time boundary ellbase ADMconstraints pugh pughinterp pughreduce pughslab cartgrid3d einstein ADM RotatingDBHIVP ioascii iobasic ioutil"
+
+# GENERAL 
+
+einstein::evolution_system = "ADM"     
+
+driver::global_nx =  32
+driver::global_ny =  32
+driver::global_nz =  32
+
+grid::type = "byspacing"
+grid::dxyz = 0.4
+grid::domain = "full"
+
+time::dtfac = 0.25   
+
+cactus::cctk_initial_time = 0.   
+cactus::cctk_itlast       = 1
+
+adm::method           = "stagleap" 
+adm::bound            = "static"
+admconstraints::bound = "static"
+
+# MODEL 
+
+einstein::initial_data = "rotatingdbh"
+rotatingdbhivp::amp      = 0.5
+rotatingdbhivp::eta0     = 0.0
+rotatingdbhivp::sigma    = 1.0
+rotatingdbhivp::byJ      = 10.0
+rotatingdbhivp::mm       = 2.0
+rotatingdbhivp::etamax   = 6.0
+rotatingdbhivp::neta     = 102
+rotatingdbhivp::ntheta   = 27
+rotatingdbhivp::nphi     = 5
+rotatingdbhivp::interpolation_order = 1
+
+# GAUGE
+einstein::slicing    = "geodesic"  
+
+# OUTPUT ########################################################
+IO::outdir         = "test_rdbh"
+IOASCII::out1d_every     =  1
+IOBasic::outScalar_vars  = "einstein::metric einstein::curv admconstraints::ham"
+IOASCII::out1D_vars      = "einstein::metric einstein::curv einstein::lapse admconstraints::ham einstein::spherical_metric einstein::spherical_curv einstein::psi einstein::confac_1derivs einstein::confac_2derivs"
+IOBasic::outInfo_every   = 1
+IOBasic::outInfo_vars    = "admconstraints::ham einstein::gxx"
+
+##################################################################
+