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

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
ML_BSSN::fdOrder             = 8

GenericFD::jacobian_group = "Coordinates::jacobian"
GenericFD::jacobian_derivative_group = "Coordinates::jacobian2"
GenericFD::jacobian_identity_map = 0

Boundary::radpower                     = 2

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

Carpet::domain_from_multipatch          = yes
CartGrid3D::type                        = "multipatch"
CartGrid3D::set_coordinate_ranges_on    = "all maps"
Coordinates::coordinate_system          = "Thornburg04"
Coordinates::h_cartesian                = 0.06
Coordinates::h_radial                   = 0.15

Coordinates::sphere_inner_radius        = 0.6
Coordinates::sphere_outer_radius        = 1.8
Coordinates::n_angular                  = 10

Coordinates::patch_boundary_size        = 5
Coordinates::additional_overlap_size    = 2
Coordinates::outer_boundary_size        = 5

Time::timestep_method                   = "given"
Time::timestep                          = 0.03

Interpolate::interpolator_order         = 5 # Interpatch interpolation

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

Carpet::ghost_size                      = 5
Carpet::init_3_timelevels               = "no"
#CarpetLib::max_memory_size_MB           = 500
Carpet::poison_new_timelevels           = yes
CarpetLib::poison_value                 = 0 # This is needed because the intersection of the outer and interpatch boundary is not currently initialised by anything

NaNChecker::check_every                 = 1
NaNChecker::verbose                     = "all"
NaNChecker::action_if_found             = "just warn"
NaNChecker::check_vars                  = "ADMBase::metric ADMBase::curv ADMBase::lapse ADMBase::shift 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_cons_detg ML_BSSN::ML_cons_Gamma ML_BSSN::ML_cons_traceA"
NaNChecker::report_max                  = 10

#############################################################
# 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                       = 2
# 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                   = "schwarzschild"
ADMBase::initial_lapse                  = "one"
ADMBase::initial_shift                  = "zero"
ADMBase::initial_dtlapse                = "zero"
ADMBase::initial_dtshift                = "zero"
ADMBase::metric_type                    = "physical"

IDAnalyticBH::mass = 1
IDAnalyticBH::epsilon = 1e-3

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

IO::out_dir                   = $parfile
IO::out_fileinfo              = "none"
IO::parfile_write             = "no"
CarpetIOBasic::outInfo_every        = 1
CarpetIOBasic::outInfo_vars         = "ML_BSSN::ML_log_confac"

IOASCII::out1D_every          = 1
IOASCII::out_precision        = 19
IOASCII::out1D_x              = "no"
IOASCII::out1D_y              = "no"
IOASCII::out1D_z              = "yes"
IOASCII::out1D_d              = "no"
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::At11rhs"