ActiveThorns = "CoordBase SymBase Boundary NanChecker CartGrid3d Time MoL CarpetIOBasic CarpetIOScalar IOUtil Carpet CarpetLib CarpetReduce CarpetInterp CarpetSlab CarpetIOASCII ADMBase StaticConformal SpaceMask Slab Periodic ShiftedGaugeWave GenericFD CoordGauge ADMCoupling LoopControl ML_BSSN ML_BSSN_Helper GenericFD TMuNuBase SphericalSurface ADMMacros TimerReport"

ADMBase::evolution_method         = "ML_BSSN"
ADMBase::lapse_evolution_method   = "ML_BSSN"
ADMBase::shift_evolution_method   = "ML_BSSN"
ADMBase::dtlapse_evolution_method = "ML_BSSN"
ADMBase::dtshift_evolution_method = "ML_BSSN"

# We use the BBH gauge to ensure it is regression-tested.  Also, there
# is no harmonic shift in McLachlan
ML_BSSN::harmonicN           = 1      # 1+log
ML_BSSN::harmonicF           = 2.0    # 1+log
ML_BSSN::ShiftGammaCoeff     = 0.75
ML_BSSN::BetaDriver          = 1.0
ML_BSSN::LapseAdvectionCoeff = 1.0
ML_BSSN::ShiftAdvectionCoeff = 1.0

ML_BSSN::MinimumLapse        = 1.0e-8
ML_BSSN::conformalMethod = 1 # 1 for W, 0 for phi
#ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas"
ML_BSSN::my_rhs_boundary_condition     = "static"

ML_BSSN::apply_dissipation   = "always"
ML_BSSN::EpsDiss             = 0.1

Boundary::radpower                     = 2

#############################################################
# Grid
#############################################################

CoordBase::domainsize                   = minmax

CoordBase::boundary_size_x_lower        = 3
CoordBase::boundary_size_y_lower        = 3
CoordBase::boundary_size_z_lower        = 3
CoordBase::boundary_shiftout_x_lower    = 1
CoordBase::boundary_shiftout_y_lower    = 1
CoordBase::boundary_shiftout_z_lower    = 1

CoordBase::boundary_size_x_upper        = 3
CoordBase::boundary_size_y_upper        = 3
CoordBase::boundary_size_z_upper        = 3
CoordBase::boundary_shiftout_x_upper    = 0
CoordBase::boundary_shiftout_y_upper    = 0
CoordBase::boundary_shiftout_z_upper    = 0

CartGrid3D::type                        = "coordbase"
CartGrid3D::domain                      = "full"
CartGrid3D::avoid_origin                = "no"

Periodic::periodic = "yes"

CoordBase::xmin                         = 0
CoordBase::ymin                         = 0
CoordBase::zmin                         = 0

CoordBase::xmax                         = 1
CoordBase::ymax                         = 1
CoordBase::zmax                         = 1

CoordBase::dx                           = 0.1
CoordBase::dy                           = 0.1
CoordBase::dz                           = 0.1

#############################################################
# Carpet
#############################################################

Carpet::ghost_size                      = 3
Carpet::domain_from_coordbase           = "yes"
Carpet::init_3_timelevels               = "no"

#############################################################
# Timers
#############################################################

# Do not commit these as testsuite output as they will change on
# different machines.  They are useful for performance monitoring in
# automated build and test systems.
TimerReport::output_all_timers_readable = yes
TimerReport::n_top_timers               = 40
TimerReport::output_schedule_timers     = no

#############################################################
# Time integration
#############################################################

Cactus::terminate                         = "iteration"
Cactus::cctk_itlast                       = 0
# Use 0.5 for RK4 and 0.25 for ICN
Time::dtfac                           = 0.5

MethodOfLines::ode_method             = "RK4"
MethodOfLines::MoL_Intermediate_Steps = 4

#MethodOfLines::ode_method             = "generic"
#MethodOfLines::generic_type             = "RK"
#MethodOfLines::MoL_Intermediate_Steps = 1

MethodOfLines::MoL_Num_Scratch_Levels = 1
MethodOfLines::MoL_NaN_Check          = "no"
MethodOfLines::initial_data_is_crap   = "no"

#############################################################
# Initial data
#############################################################

ADMBase::initial_data = "ShiftedGaugeWave"
ADMBase::initial_lapse = "ShiftedGaugeWave"
ADMBase::initial_shift = "ShiftedGaugeWave"
ADMBase::initial_dtlapse = "ShiftedGaugeWave"
ADMBase::initial_dtshift = "ShiftedGaugeWave"

ShiftedGaugeWave::amp = 0.1
ShiftedGaugeWave::period = 0.577350269189626

ShiftedGaugeWave::psi = -1.9216757376671543544
ShiftedGaugeWave::theta = 0.66214523564555227398
ShiftedGaugeWave::phi = 1.2199169159226388841

#############################################################
# Output
#############################################################

IO::out_dir                   = $parfile
IO::out_fileinfo              = "none"
IO::parfile_write             = "no"
CarpetIOBasic::outInfo_every        = 1
CarpetIOBasic::outInfo_vars         = "ADMBase::alp"

IOASCII::out1D_every                    = 1
IOASCII::out_precision        = 19
IOASCII::out1D_d                        = no
IOASCII::out1D_x                        = yes
IOASCII::out1D_y                        = yes
IOASCII::out1D_z                        = yes
IOASCII::out1D_vars                     = "ADMBase::gxx ADMBase::kxx ML_BSSN::ML_cons_detg ML_BSSN::ML_cons_Gamma ML_BSSN::ML_cons_traceA ML_BSSN::ML_curv ML_BSSN::ML_dtlapse ML_BSSN::ML_dtshift ML_BSSN::ML_Gamma ML_BSSN::ML_lapse ML_BSSN::ML_log_confac ML_BSSN::ML_metric ML_BSSN::ML_shift ML_BSSN::ML_trace_curv ML_BSSN::ML_curvrhs ML_BSSN::ML_dtlapserhs ML_BSSN::ML_dtshiftrhs ML_BSSN::ML_Gammarhs ML_BSSN::ML_lapserhs ML_BSSN::ML_log_confacrhs ML_BSSN::ML_metricrhs ML_BSSN::ML_shiftrhs ML_BSSN::ML_trace_curvrhs"

Carpet::poison_new_timelevels           = yes