# 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"