################################################################################
  Cactus::cctk_run_title = "CCZ4 Binary Black Hole"
################################################################################

ActiveThorns = "
  ADMBase
  ADMCoupling
  ADMMacros
  AEILocalInterp
  AHFinderDirect
  BLAS
  Boundary
  Carpet
  CarpetInterp
  CarpetIOASCII
  CarpetIOBasic
  CarpetIOHDF5
  CarpetIOScalar
  CarpetLib
  CarpetMask
  CarpetReduce
  CarpetRegrid2
  CarpetSlab
  CarpetTracker
  CartGrid3D
  CoordBase
  CoordGauge
  Dissipation
  Formaline
  Fortran
  GenericFD
  GSL
  HDF5
  InitBase
  IOUtil
  LAPACK
  LocalInterp
  LoopControl
  ML_ADMConstraints
  ML_BSSN
  ML_BSSN_Helper
  MoL
  Multipole
  NaNChecker
  NewRad
  PunctureTracker
  QuasiLocalMeasures
  ReflectionSymmetry
  RotatingSymmetry180
  Slab
  SpaceMask
  SphericalSurface
  StaticConformal
  SummationByParts
  SymBase
  SystemStatistics
  TerminationTrigger
  Time
  TimerReport
  TmunuBase
  TwoPunctures
  Vectors
  WeylScal4
"

################################################################################
# Grid structure
################################################################################

CartGrid3D::type 			  = "coordbase"
Carpet::domain_from_coordbase = "yes"
CoordBase::domainsize 		  = "minmax"

CoordBase::xmin =    0.00
CoordBase::ymin = -120.00
CoordBase::zmin =    0.00
CoordBase::xmax = +120.00
CoordBase::ymax = +120.00
CoordBase::zmax = +120.00
CoordBase::dx   =    2.00
CoordBase::dy   =    2.00
CoordBase::dz   =    2.00

Driver::ghost_size       = 3

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_z_lower = 1

CarpetRegrid2::symmetry_rotating180  = yes

ReflectionSymmetry::reflection_z     = yes
ReflectionSymmetry::avoid_origin_z   = no

Time::dtfac = 0.25

################################################################################
# Mesh refinement
################################################################################

Carpet::max_refinement_levels           = 7
CarpetRegrid2::num_centres              = 2
CarpetRegrid2::num_levels_1             = 7
CarpetRegrid2::position_x_1             = 3
CarpetRegrid2::radius_1             [1] = 24.0
CarpetRegrid2::radius_1             [2] = 12.0
CarpetRegrid2::radius_1             [3] = 6.0
CarpetRegrid2::radius_1             [4] = 3.0
CarpetRegrid2::radius_1             [5] = 1.5
CarpetRegrid2::radius_1             [6] = 0.6
CarpetRegrid2::num_levels_2             = 7
CarpetRegrid2::position_x_2             = -3
CarpetRegrid2::radius_2             [1] = 24.0
CarpetRegrid2::radius_2             [2] = 12.0
CarpetRegrid2::radius_2             [3] = 6.0
CarpetRegrid2::radius_2             [4] = 3.0
CarpetRegrid2::radius_2             [5] = 1.5
CarpetRegrid2::radius_2             [6] = 0.6

Carpet::use_buffer_zones                = yes
Carpet::prolongation_order_space        = 5
Carpet::prolongation_order_time         = 2

CarpetRegrid2::regrid_every             = 64
CarpetRegrid2::verbose                  = yes
Carpet::grid_coordinates_filename       = "carpet-grid.asc"

Carpet::time_refinement_factors = "[1,1,2,4,8,16,32,64,128,256]"

###############################################################################
# Initial Data
###############################################################################

ADMBase::initial_data    = "twopunctures"
ADMBase::initial_lapse   = "twopunctures-averaged"
ADMBase::initial_shift   = "zero"
ADMBase::initial_dtlapse = "zero"
ADMBase::initial_dtshift = "zero"

TwoPunctures::target_M_plus             = 0.5
TwoPunctures::target_M_minus            = 0.5

TwoPunctures::par_m_plus                = 0.476534633024028
TwoPunctures::par_m_minus               = 0.476534633024028

TwoPunctures::par_b                     = 3
TwoPunctures::center_offset[0]          = 0

TwoPunctures::par_P_plus[0]             = -0.0058677669328272
TwoPunctures::par_P_plus[1]             = 0.138357448824906
TwoPunctures::par_P_plus[2]             = 0.

TwoPunctures::par_P_minus[0]            = 0.0058677669328272
TwoPunctures::par_P_minus[1]            = -0.138357448824906
TwoPunctures::par_P_minus[2]            = 0.

TwoPunctures::par_S_plus[0]             = 0.
TwoPunctures::par_S_plus[1]             = 0.
TwoPunctures::par_S_plus[2]             = 0.

TwoPunctures::par_S_minus[0]            = 0.
TwoPunctures::par_S_minus[1]            = 0.
TwoPunctures::par_S_minus[2]            = 0.

TwoPunctures::give_bare_mass            = yes
TwoPunctures::TP_epsilon = 1.0e-6

Carpet::init_fill_timelevels = yes
InitBase::initial_data_setup_method     = "init_single_level"

###############################################################################
# Evolution
###############################################################################

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"

ML_BSSN::GammaShift          = 1.0
ML_BSSN::dampk1              = 1.0
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
ML_BSSN::formulation           = 1 # 1 for CCZ4
ML_BSSN::apply_dissipation     = never
ML_BSSN::dt_lapse_shift_method = "noLapseShiftAdvection"

################################################################################
# Boundary conditions
################################################################################

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"
ML_BSSN::ML_Theta_bound      = "none"

################################################################################
# BH tracking
################################################################################

CarpetTracker::surface                      [0] = 0
CarpetTracker::surface                      [1] = 1
PunctureTracker::track                      [0] = yes
PunctureTracker::initial_x                  [0] = 3
PunctureTracker::which_surface_to_store_info[0] = 0
PunctureTracker::track                      [1] = yes
PunctureTracker::initial_x                  [1] = -3
PunctureTracker::which_surface_to_store_info[1] = 1

###############################################################################
# Spatial finite differencing
###############################################################################

SummationByParts::order = 4
ML_BSSN::fdOrder        = 4

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
  ML_BSSN::ML_Theta
"

###############################################################################
# Time Integration
###############################################################################

MoL::ODE_Method             = "RK4"
MoL::MoL_Intermediate_Steps = 4
MoL::MoL_Num_Scratch_Levels = 1

################################################################################
# Psi4 mode decomposition by Multipole
################################################################################

WeylScal4::fd_order = "4th"

Multipole::nradii       = 4
Multipole::radius[0]    = 30
Multipole::radius[1]    = 50
Multipole::radius[2]    = 80
Multipole::radius[3]    = 100
Multipole::ntheta       = 120
Multipole::nphi         = 240
Multipole::variables    = "
  WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i' name='psi4'}
"

Multipole::out_every    = 32
Multipole::l_max        = 8
Multipole::output_hdf5  = yes
Multipole::output_ascii = no

################################################################################
# Apparent Horizons
################################################################################


AHFinderDirect::N_horizons                               = 3
AHFinderDirect::find_every                               = 128
AHFinderDirect::output_h_every                           = 0
AHFinderDirect::max_Newton_iterations__initial           = 50
AHFinderDirect::max_Newton_iterations__subsequent        = 50
AHFinderDirect::max_allowable_Theta_growth_iterations    = 10
AHFinderDirect::max_allowable_Theta_nonshrink_iterations = 10
AHFinderDirect::geometry_interpolator_name               = "Lagrange polynomial interpolation"
AHFinderDirect::geometry_interpolator_pars               = "order=4"
AHFinderDirect::surface_interpolator_name                = "Lagrange polynomial interpolation"
AHFinderDirect::surface_interpolator_pars                = "order=4"
AHFinderDirect::verbose_level                            = "physics details"
AHFinderDirect::move_origins                             = yes

AHFinderDirect::origin_x                             [1] = 3
AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 3
AHFinderDirect::initial_guess__coord_sphere__radius  [1] = 0.25
AHFinderDirect::which_surface_to_store_info          [1] = 2 
AHFinderDirect::set_mask_for_individual_horizon      [1] = no
AHFinderDirect::reset_horizon_after_not_finding      [1] = no
AHFinderDirect::track_origin_from_grid_scalar        [1] = yes
AHFinderDirect::track_origin_source_x                [1] = "PunctureTracker::pt_loc_x[0]"
AHFinderDirect::track_origin_source_y                [1] = "PunctureTracker::pt_loc_y[0]"
AHFinderDirect::track_origin_source_z                [1] = "PunctureTracker::pt_loc_z[0]"
AHFinderDirect::max_allowable_horizon_radius         [1] = 3

AHFinderDirect::origin_x                             [2] = -3
AHFinderDirect::initial_guess__coord_sphere__x_center[2] = -3
AHFinderDirect::initial_guess__coord_sphere__radius  [2] = 0.25
AHFinderDirect::which_surface_to_store_info          [2] = 3 
AHFinderDirect::set_mask_for_individual_horizon      [2] = no
AHFinderDirect::reset_horizon_after_not_finding      [2] = no
AHFinderDirect::track_origin_from_grid_scalar        [2] = yes
AHFinderDirect::track_origin_source_x                [2] = "PunctureTracker::pt_loc_x[1]"
AHFinderDirect::track_origin_source_y                [2] = "PunctureTracker::pt_loc_y[1]"
AHFinderDirect::track_origin_source_z                [2] = "PunctureTracker::pt_loc_z[1]"
AHFinderDirect::max_allowable_horizon_radius         [2] = 3

AHFinderDirect::origin_x                             [3] = 0
AHFinderDirect::find_after_individual                [3] = 0
AHFinderDirect::initial_guess__coord_sphere__x_center[3] = 0
AHFinderDirect::initial_guess__coord_sphere__radius  [3] = 1.0
AHFinderDirect::which_surface_to_store_info          [3] = 4 
AHFinderDirect::set_mask_for_individual_horizon      [3] = no
AHFinderDirect::max_allowable_horizon_radius         [3] = 6

################################################################################
# Spherical surfaces
################################################################################

SphericalSurface::nsurfaces = 5
SphericalSurface::maxntheta = 66
SphericalSurface::maxnphi   = 124
SphericalSurface::verbose               = no

# Surfaces 0 and 1 are used by PunctureTracker

# Horizon 1
SphericalSurface::ntheta            [2] = 41
SphericalSurface::nphi              [2] = 80
SphericalSurface::nghoststheta      [2] = 2
SphericalSurface::nghostsphi        [2] = 2

# Horizon 2
SphericalSurface::ntheta            [3] = 41
SphericalSurface::nphi              [3] = 80
SphericalSurface::nghoststheta      [3] = 2
SphericalSurface::nghostsphi        [3] = 2

# Common horizon
SphericalSurface::ntheta            [4] = 41
SphericalSurface::nphi              [4] = 80
SphericalSurface::nghoststheta      [4] = 2
SphericalSurface::nghostsphi        [4] = 2

################################################################################
# Isolated Horizons
################################################################################

QuasiLocalMeasures::verbose                = yes
QuasiLocalMeasures::veryverbose            = no
QuasiLocalMeasures::interpolator           = "Lagrange polynomial interpolation"
QuasiLocalMeasures::interpolator_options   = "order=4"
QuasiLocalMeasures::spatial_order          = 4
QuasiLocalMeasures::num_surfaces           = 3
QuasiLocalMeasures::surface_index      [0] = 2
QuasiLocalMeasures::surface_index      [1] = 3
QuasiLocalMeasures::surface_index      [2] = 4

################################################################################
# Correctness Checking
################################################################################

Carpet::poison_new_timelevels = yes
Carpet::poison_value          = 100
CarpetLib::poison_new_memory  = yes
CarpetLib::poison_value       = 200
Carpet::check_for_poison      = no

NaNChecker::check_every     = 128
NanChecker::check_after     = 0
NaNChecker::report_max      = 10
NaNChecker::verbose         = "all"
NaNChecker::action_if_found = "terminate"
NaNChecker::out_NaNmask     = yes
NaNChecker::check_vars      = "
  ML_BSSN::ML_log_confac
"

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

Cactus::cctk_timer_output               = full
TimerReport::out_every                  = 256
TimerReport::n_top_timers               = 40
TimerReport::output_all_timers_together = yes
TimerReport::output_schedule_timers     = no

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

IO::out_dir                  = $parfile
IOScalar::one_file_per_group = yes
IOASCII::one_file_per_group  = yes

IOBasic::outInfo_every      = 1
IOBasic::outInfo_reductions = "minimum maximum"
IOBasic::outInfo_vars       = "
  Carpet::physical_time_per_hour
  ML_BSSN::ML_log_confac
  SystemStatistics::maxrss_mb
"

IOScalar::outScalar_every               = 128
IOScalar::outScalar_reductions          = "minimum maximum average"
IOScalar::outScalar_vars                = "SystemStatistics::process_memory_mb"

IOASCII::out0D_every = 32
IOASCII::out0D_vars  = "
  Carpet::timing
  PunctureTracker::pt_loc
  QuasiLocalMeasures::qlm_scalars{out_every = 128}
"

IOHDF5::out1D_every                    = 256
IOHDF5::out1D_d                        = no
IOHDF5::out1D_x                        = yes
IOHDF5::out1D_y                        = no
IOHDF5::out1D_z                        = no
IOHDF5::out1D_vars                     = "
  ML_BSSN::ML_log_confac
"

IOHDF5::out2D_every                    = 1024
IOHDF5::out2D_xy                       = yes
IOHDF5::out2D_xz                       = no
IOHDF5::out2D_yz                       = no
IOHDF5::out2D_vars                     = "
  ML_BSSN::ML_log_confac
"

IOHDF5::out_every                      = 8192
IOHDF5::out_vars                       = "
  ML_BSSN::ML_log_confac
"

Carpet::verbose           = no
Carpet::veryverbose       = no
Carpet::schedule_barriers = no
Carpet::storage_verbose   = no
Carpet::timers_verbose    = no
CarpetLib::output_bboxes  = no

Cactus::cctk_full_warnings         = yes
Cactus::highlight_warning_messages = no

###############################################################################
# Checkpointing and recovery
###############################################################################

CarpetIOHDF5::checkpoint            		= yes
IO::checkpoint_ID                   		= no
IO::recover     							= "autoprobe"
IO::checkpoint_every_walltime_hours 		= 3
IO::out_proc_every                          = 2
IO::checkpoint_keep                         = 1
IO::checkpoint_on_terminate         		= yes
IO::checkpoint_dir                  		= "../checkpoints"
IO::recover_dir                     		= "../checkpoints"
IO::abort_on_io_errors                      = yes
CarpetIOHDF5::open_one_input_file_at_a_time = yes
CarpetIOHDF5::compression_level             = 9

################################################################################
# Run termination
################################################################################

TerminationTrigger::max_walltime                 = @WALLTIME_HOURS@   # hours
TerminationTrigger::on_remaining_walltime        = 30 # minutes
TerminationTrigger::output_remtime_every_minutes = 30
Cactus::terminate       = "time"
Cactus::cctk_final_time = 300.0