diff options
Diffstat (limited to 'test/teukolskyParity.par')
| -rw-r--r-- | test/teukolskyParity.par | 105 |
1 files changed, 105 insertions, 0 deletions
diff --git a/test/teukolskyParity.par b/test/teukolskyParity.par new file mode 100644 index 0000000..4a12d13 --- /dev/null +++ b/test/teukolskyParity.par @@ -0,0 +1,105 @@ +# Interpolation of psis to 2d array. Check tensorparities, especially on reflection +# Author: Mitica Vulcanov <vulcan@aei.mpg.de> +# Author: Roland Haas <roland.haas@physics.gatech.edu> +# Author: Tanja Bode <tanja.bode@physics.gatech.edu> +# $Header$ + +# Required thorns +ActiveThorns = "CartGrid3D CoordBase SymBase Boundary Time Carpet CarpetLib CarpetSlab CarpetRegrid2 CarpetReduce InitBase CarpetIOAscii IOUtil IOBasic ADMBase ADMMacros ADMConstraints GenericFD TmunuBase ADMCoupling StaticConformal CoordGauge SpaceMask IDLinearWaves LoopControl GSL WeylScal4 ReflectionSymmetry MoL InterpToArray CarpetInterp AEILocalInterp" + +# Grid +driver::ghost_size = 2 + +CartGrid3D::avoid_origin = no + +Carpet::max_refinement_levels = 3 +Carpet::prolongation_order_space = 3 # 1 seems to be buggy in Carpet/git +Carpet::prolongation_order_time = 0 +Carpet::domain_from_coordbase = yes + +Carpet::poison_new_timelevels = "yes" +Carpet::check_for_poison = "no" +CarpetLib::poison_new_memory = "yes" +CarpetLib::poison_value = 114 + +InitBase::initial_data_setup_method = init_some_levels +Carpet::init_fill_timelevels = yes + +CarpetRegrid2::min_distance = 0 +CarpetRegrid2::ensure_proper_nesting = "yes" +CarpetRegrid2::regrid_every = 0 +CarpetRegrid2::num_centres = 1 +CarpetRegrid2::num_levels_1 = 3 +CarpetRegrid2::position_x_1 = -0.0625 +CarpetRegrid2::position_y_1 = -0.0625 +CarpetRegrid2::position_z_1 = 0 +CarpetRegrid2::radius_1[1] = 2.0 +CarpetRegrid2::radius_1[2] = 1.0 + +Grid::type = "coordbase" + +CoordBase::xmin = -4.0625 +CoordBase::ymin = -4.0625 +CoordBase::zmin = 0.0 +CoordBase::xmax = 3.9375 +CoordBase::ymax = 3.9375 +CoordBase::zmax = 4.0 +CoordBase::dx = 0.5 +CoordBase::dy = 0.5 +CoordBase::dz = 0.5 + +CoordBase::boundary_size_x_lower = 2 +CoordBase::boundary_size_y_lower = 2 +CoordBase::boundary_size_z_lower = 2 +CoordBase::boundary_size_x_upper = 2 +CoordBase::boundary_size_y_upper = 2 +CoordBase::boundary_size_z_upper = 2 + +CoordBase::boundary_shiftout_z_lower = 1 +CoordBase::boundary_shiftout_z_upper = 0 + +ReflectionSymmetry::avoid_origin_z = "no" +ReflectionSymmetry::reflection_z = "yes" + +cactus::cctk_initial_time = 0 +cactus::cctk_itlast = 0 + +# Initial data +ADMBase::initial_data = "teukwaves" +ADMBase::initial_shift = "zero" +ADMBase::initial_lapse = "one" +IDLinearWaves::mvalue = 2 + +# Gauge +ADMBase::lapse_evolution_method = "static" +ADMBase::shift_evolution_method = "static" +ADMBase::dtshift_evolution_method = "static" + +# Evolution +ADMBase::evolution_method = "static" + +ADMBase::lapse_timelevels = 1 +ADMBase::shift_timelevels = 1 +ADMBase::metric_timelevels = 1 + +# Output +IO::out_dir = $parfile +IO::out_fileinfo = "none" +IO::parfile_write = "no" + +InterpToArray::narrays2d=3 +InterpToArray::array2d_vars[0] = "weylscal4::psi4r" +InterpToArray::array2d_vars[1] = "weylscal4::psi4i" +InterpToArray::array2d_vars[2] = "weylscal4::psi3r" +InterpToArray::array2d_x0 = -2.0 +InterpToArray::array2d_y0 = -2.0 +InterpToArray::array2d_z0 = -2.0 +InterpToArray::array2d_dx_i = 0.3 +InterpToArray::array2d_npoints_i = 14 +InterpToArray::array2d_dz_j = 0.3 +InterpToArray::array2d_npoints_j = 14 +IOASCII::out2d_every = 1 +IOASCII::out2d_vars = "interptoarray::arrays2d" + +WeylScal4::timelevels = 1 +WeylScal4::calc_scalars = "psis" |