diff options
author | guzman <guzman@535fb057-194f-0410-b5a5-c63992f15602> | 2002-05-03 16:04:59 +0000 |
---|---|---|
committer | guzman <guzman@535fb057-194f-0410-b5a5-c63992f15602> | 2002-05-03 16:04:59 +0000 |
commit | 89a63ca71255b010277c97bacbcd5f0d1bc65631 (patch) | |
tree | bc2fc9f25d34c794e1ad04d20a030df206c1eef8 | |
parent | c198eb23b852b1b6ae6485fe8c04812128c344b4 (diff) |
*** empty log message ***
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/RotatingDBHIVP/trunk@30 535fb057-194f-0410-b5a5-c63992f15602
-rw-r--r-- | src/ParamCheck.c | 107 | ||||
-rwxr-xr-x | test/xeon_test_rdbh.par | 54 |
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" + +################################################################## + |