Move interpolation to own test case, adding two more variables.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/WeylScal4/trunk@117 4f5cb9a8-4dd8-4c2d-9bbd-173fa4467843

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"