##############################################################################
#                              GRID SETUP                                    #
##############################################################################
ActiveThorns = "CoordBase"

CoordBase::domainsize = "minmax"
CoordBase::xmin = -256
CoordBase::xmax = 256
CoordBase::ymin = 0
CoordBase::ymax = 256
CoordBase::zmin = 0
CoordBase::zmax = 256

CoordBase::spacing = "numcells"
CoordBase::ncells_x = 64
CoordBase::ncells_y = 32
CoordBase::ncells_z = 32

CoordBase::boundary_size_x_lower     = 3
CoordBase::boundary_size_y_lower     = 3
CoordBase::boundary_size_z_lower     = 3
CoordBase::boundary_size_x_upper     = 3
CoordBase::boundary_size_y_upper     = 3
CoordBase::boundary_size_z_upper     = 3

CoordBase::boundary_shiftout_x_lower = 1
CoordBase::boundary_shiftout_y_lower = 1
CoordBase::boundary_shiftout_z_lower = 1

ActiveThorns = "CartGrid3d"

CartGrid3d::type                 = "coordbase"

ActiveThorns = "Carpet CarpetLib CarpetRegrid2 CarpetTracker SphericalSurface IOUtil"

Carpet::verbose                  = 0
Carpet::domain_from_coordbase    = 1
Carpet::init_fill_timelevels     = 1
Carpet::poison_new_timelevels    = 1
Carpet::max_refinement_levels    = 10
Carpet::refinement_factor        = 2       # anything else not implemented
Carpet::refinement_centering     = "vertex"
Carpet::ghost_size               = 3
Carpet::use_buffer_zones         = 1
Carpet::prolongation_order_space = 5
Carpet::prolongation_order_time  = 2
Carpet::num_threads              = 4
CarpetLib::poison_new_memory     = 1

CarpetRegrid2::verbose                 = 1
CarpetRegrid2::regrid_every            = 256
CarpetRegrid2::num_centres             = 1
CarpetRegrid2::movement_threshold_1    = 0.1
CarpetRegrid2::position_x_1            = 0
CarpetRegrid2::num_levels_1            = 8
CarpetRegrid2::radius_1[1]             = 128
CarpetRegrid2::radius_1[2]             = 64
CarpetRegrid2::radius_1[3]             = 32
CarpetRegrid2::radius_1[4]             = 16
CarpetRegrid2::radius_1[5]             = 8
CarpetRegrid2::radius_1[6]             = 4
CarpetRegrid2::radius_1[7]             = 2
CarpetRegrid2::radius_1[8]             = 1
CarpetRegrid2::radius_1[9]             = 0.5

CarpetTracker::verbose                 = 0
CarpetTracker::surface[0]              = 1

ActiveThorns = "ReflectionSymmetry"

ReflectionSymmetry::reflection_z   = 1
ReflectionSymmetry::avoid_origin_z = 0
ReflectionSymmetry::reflection_y   = 1
ReflectionSymmetry::avoid_origin_y = 0
#ReflectionSymmetry::reflection_x   = 1
#ReflectionSymmetry::avoid_origin_x = 0


##############################################################################
#                            INITIAL DATA                                    #
##############################################################################

ActiveThorns = "ADMBase ADMCoupling ADMMacros StaticConformal Trumpet"

ADMBase::metric_timelevels        = 3
ADMBase::initial_data             = "trumpet"
ADMBase::initial_lapse            = "trumpet"
ADMBase::initial_dtlapse          = "zero"
ADMBase::initial_shift            = "trumpet"
ADMBase::initial_dtshift          = "zero"

ADMMacros::spatial_order = 4

Trumpet::boost_velocity = 0.5

##############################################################################
#                             EVOLUTION                                      #
##############################################################################

Cactus::terminate       = "time"
Cactus::cctk_final_time = 100

ActiveThorns = "MoL Time"

MoL::initial_data_is_crap   = 0
MoL::ODE_Method             = "RK3"
MoL::MoL_Intermediate_Steps = 3
MoL::MoL_Num_Scratch_Levels = 1
MoL::MoL_NaN_Check          = 1

Time::dtfac = 0.125

ActiveThorns = "ML_BSSN ML_BSSN_Helper CoordGauge GenericFD NewRad Boundary StaticConformal TmunuBase
                SymBase LoopControl"

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

ML_BSSN::verbose                       = 1
ML_BSSN::timelevels                    = 3
ML_BSSN::fdOrder                       = 4
ML_BSSN::harmonicF                     = 2      # \dot{α} = - f α^n K
ML_BSSN::harmonicN                     = 1      # 1+log
ML_BSSN::conformalMethod               = 1      # χ method
ML_BSSN::ShiftGammaCoeff               = 0.75
ML_BSSN::HarmonicShift                 = 0
ML_BSSN::LapseACoeff                   = 0.0
ML_BSSN::ShiftBCoeff                   = 1.0
ML_BSSN::AlphaDriver                   = 0.0
ML_BSSN::BetaDriver                    = 2.0
ML_BSSN::LapseAdvectionCoeff           = 1.0
ML_BSSN::ShiftAdvectionCoeff           = 1.0
ML_BSSN::MinimumLapse                  = 1e-08
ML_BSSN::dt_lapse_shift_method         = "correct"
ML_BSSN::EpsDiss                       = 0
ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas"
ML_BSSN::my_rhs_boundary_condition     = "NewRad"
Boundary::radpower                     = 2
ML_BSSN::ML_log_confac_bound           = "none"
ML_BSSN::ML_metric_bound               = "none"
ML_BSSN::ML_Gamma_bound                = "none"
ML_BSSN::ML_trace_curv_bound           = "none"
ML_BSSN::ML_curv_bound                 = "none"
ML_BSSN::ML_lapse_bound                = "none"
ML_BSSN::ML_dtlapse_bound              = "none"
ML_BSSN::ML_shift_bound                = "none"
ML_BSSN::ML_dtshift_bound              = "none"

ActiveThorns = "Dissipation SpaceMask"

Dissipation::order = 5
Dissipation::vars  = "
        ML_BSSN::ML_log_confac
        ML_BSSN::ML_metric
        ML_BSSN::ML_trace_curv
        ML_BSSN::ML_curv
        ML_BSSN::ML_Gamma
        ML_BSSN::ML_lapse
        ML_BSSN::ML_shift
        ML_BSSN::ML_dtlapse
        ML_BSSN::ML_dtshift
"
ActiveThorns = "PunctureTracker"

PunctureTracker::track[0] = 1
PunctureTracker::initial_x[0] = 0
PunctureTracker::which_surface_to_store_info[0] = 1

###############################################################################
#                               ANALYSIS                                      #
###############################################################################

ActiveThorns = "NaNChecker CarpetReduce"

NaNChecker::check_every     = 1024
NaNChecker::action_if_found = "just warn"
NaNChecker::check_vars      = "
        ML_BSSN::ML_log_confac
        ML_BSSN::ML_metric
        ML_BSSN::ML_trace_curv
        ML_BSSN::ML_curv
        ML_BSSN::ML_Gamma
        ML_BSSN::ML_lapse
        ML_BSSN::ML_shift
        ML_BSSN::ML_dtlapse
        ML_BSSN::ML_dtshift
        ADMBase::metric
        ADMBase::curv
        ADMBase::lapse
        ADMBase::shift
        ADMBase::dtlapse
        ADMBase::dtshift
"

SphericalSurface::nsurfaces = 6
SphericalSurface::maxntheta = 39
SphericalSurface::maxnphi   = 76

SphericalSurface::ntheta      [0] = 39
SphericalSurface::nphi        [0] = 76
SphericalSurface::nghoststheta[0] = 2
SphericalSurface::nghostsphi  [0] = 2
SphericalSurface::ntheta      [1] = 39
SphericalSurface::nphi        [1] = 76
SphericalSurface::nghoststheta[1] = 2
SphericalSurface::nghostsphi  [1] = 2
SphericalSurface::ntheta      [2] = 39
SphericalSurface::nphi        [2] = 76
SphericalSurface::nghoststheta[2] = 2
SphericalSurface::nghostsphi  [2] = 2
SphericalSurface::ntheta      [3] = 39
SphericalSurface::nphi        [3] = 76
SphericalSurface::nghoststheta[3] = 2
SphericalSurface::nghostsphi  [3] = 2
SphericalSurface::ntheta      [4] = 39
SphericalSurface::nphi        [4] = 76
SphericalSurface::nghoststheta[4] = 2
SphericalSurface::nghostsphi  [4] = 2
SphericalSurface::ntheta      [5] = 39
SphericalSurface::nphi        [5] = 76
SphericalSurface::nghoststheta[5] = 2
SphericalSurface::nghostsphi  [5] = 2

ActiveThorns = "CarpetInterp AEILocalInterp"
ActiveThorns = "AHFinderDirect CarpetInterp AEILocalInterp"

AHFinderDirect::find_every   = 256
AHFinderDirect::move_origins = 1
AHFinderDirect::predict_origin_movement = 0
AHFinderDirect::N_horizons   = 5

AHFinderDirect::origin_x[1]                              = 0
AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0
AHFinderDirect::initial_guess__coord_sphere__radius[1]   = 1
AHFinderDirect::which_surface_to_store_info[1]           = 0

AHFinderDirect::surface_selection[2]                     = "areal radius"
AHFinderDirect::origin_x[2]                              = 0
AHFinderDirect::initial_guess__coord_sphere__x_center[2] = 0
AHFinderDirect::initial_guess__coord_sphere__radius[2]   = 50
AHFinderDirect::desired_value[2]                         = 50
AHFinderDirect::which_surface_to_store_info[2]           = 2

AHFinderDirect::surface_selection[3]                     = "areal radius"
AHFinderDirect::origin_x[3]                              = 0
AHFinderDirect::initial_guess__coord_sphere__x_center[3] = 0
AHFinderDirect::initial_guess__coord_sphere__radius[3]   = 100
AHFinderDirect::desired_value[3]                         = 100
AHFinderDirect::which_surface_to_store_info[3]           = 3

AHFinderDirect::surface_selection[4]                     = "areal radius"
AHFinderDirect::origin_x[4]                              = 0
AHFinderDirect::initial_guess__coord_sphere__x_center[4] = 0
AHFinderDirect::initial_guess__coord_sphere__radius[4]   = 150
AHFinderDirect::desired_value[4]                         = 150
AHFinderDirect::which_surface_to_store_info[4]           = 4

AHFinderDirect::surface_selection[5]                     = "areal radius"
AHFinderDirect::origin_x[5]                              = 0
AHFinderDirect::initial_guess__coord_sphere__x_center[5] = 0
AHFinderDirect::initial_guess__coord_sphere__radius[5]   = 250
AHFinderDirect::desired_value[5]                         = 250
AHFinderDirect::which_surface_to_store_info[5]           = 5

ActiveThorns = "Multipole WeylScal4"

WeylScal4::verbose = 1
Multipole::integration_method = "simpson"
Multipole::nradii = 4
Multipole::radius[0] = 50
Multipole::radius[1] = 100
Multipole::radius[2] = 150
Multipole::radius[3] = 250
Multipole::variables = "WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i'}"
Multipole::l_max = 4
Multipole::interpolator_name = "Lagrange polynomial interpolation"

ActiveThorns = "ML_Kretschmann ML_Kretschmann_Helper"
ML_Kretschmann::conformalMethod = 1
ML_Kretschmann::ML_Kretschmann_kretschmann_calc_every = 128

ActiveThorns = "QuasiLocalMeasures"
QuasiLocalMeasures::num_surfaces = 5
QuasiLocalMeasures::surface_index[0] = 0
QuasiLocalMeasures::surface_index[1] = 2
QuasiLocalMeasures::surface_index[2] = 3
QuasiLocalMeasures::surface_index[3] = 4
QuasiLocalMeasures::surface_index[4] = 5

###############################################################################
#                             OUTPUT                                          #
###############################################################################

ActiveThorns = "CarpetIOBasic CarpetIOHDF5"

IO::out_dir       = $parfile
IO::out_criterion = "time"
IO::out_dt        = 0.25
IO::out_mode      = "onefile"
IO::out_unchunked = 1
IOHDF5::checkpoint = 1
IO::checkpoint_every = 4096
IO::checkpoint_on_terminate = 1
IO::checkpoint_dir = $parfile
IO::recover = "autoprobe"
IO::recover_dir = $parfile

IOBasic::outInfo_criterion  = "iteration"
IOBasic::outInfo_every      = 128
IOBasic::outInfo_reductions = "minimum maximum average norm1 norm2"
IOBasic::outInfo_vars       = "
        ML_BSSN::H
        PunctureTracker::pt_loc_x[0]
        PunctureTracker::pt_vel_x[0]
        QuasiLocalMeasures::qlm_adm_energy[0]
        QuasiLocalMeasures::qlm_adm_energy[1]
        QuasiLocalMeasures::qlm_adm_energy[2]
        QuasiLocalMeasures::qlm_adm_energy[3]
        QuasiLocalMeasures::qlm_adm_energy[4]
"

IOHDF5::compression_level      = 5
IOHDF5::output_symmetry_points = 0
IOHDF5::out3D_ghosts        = 0
IOHDF5::out3D_outer_ghosts  = 0
IOHDF5::out_vars           = "
    ML_BSSN::ML_metric
    ML_BSSN::ML_curv
    ML_BSSN::ML_trace_curv
    ML_BSSN::ML_log_confac
    ML_BSSN::ML_lapse
    ML_BSSN::ML_shift
    ML_BSSN::ML_dtshift
    ML_BSSN::ML_Ham
    ML_BSSN::ML_mom
    ML_BSSN::ML_cons_detg
    ML_BSSN::ML_cons_traceA
    ML_BSSN::ML_Gamma
    AHFinderDirect::ah_radius
    AHFinderDirect::ah_origin
    ML_Kretschmann::Kretsch
    PunctureTracker::pt_loc_x[0]
    PunctureTracker::pt_vel_x[0]
    SphericalSurface::sf_info
    SphericalSurface::sf_radius
    QuasiLocalMeasures::qlm_scalars
    "