Cactus::cctk_run_title = "CCZ4 Robust Stability Test"

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

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          = 0.0
ML_BSSN::dampk1              = 0.0
ML_BSSN::harmonicN           = 2     
ML_BSSN::harmonicF           = 1.0   
ML_BSSN::ShiftGammaCoeff     = 0.0
ML_BSSN::ShiftBCoeff         = 0.0 # B to zero
ML_BSSN::BetaDriver          = 0.0
ML_BSSN::LapseACoeff         = 0.0 # A to zero
ML_BSSN::LapseAdvectionCoeff = 0.0
ML_BSSN::ShiftAdvectionCoeff = 0.0

ML_BSSN::conformalMethod     = 1 # 1 for W
ML_BSSN::formulation         = 1 # 1 for CCZ4

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"

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

CoordBase::domainsize                   = minmax

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

CoordBase::boundary_size_x_upper        = 5
CoordBase::boundary_size_y_upper        = 5
CoordBase::boundary_size_z_upper        = 5
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.5
CoordBase::ymin                         = -0.5
CoordBase::zmin                         = -0.5

CoordBase::xmax                         = 0.5
CoordBase::ymax                         = 0.5
CoordBase::zmax                         = 0.5

CoordBase::dx                           = 0.02
CoordBase::dy                           = 0.02
CoordBase::dz                           = 0.02

#############################################################
# CarpetRegrid
#############################################################

CarpetRegrid::refinement_levels           = 1
CarpetRegrid::refined_regions             = "manual-coordinate-list"
CarpetRegrid::smart_outer_boundaries      = yes
CarpetRegrid::keep_same_grid_structure    = "yes"

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

Carpet::ghost_size                      = 5
Carpet::domain_from_coordbase           = "yes"
Carpet::max_refinement_levels           = 1
Carpet::init_3_timelevels               = "no"
Carpet::prolongation_order_space        = 3
Carpet::prolongation_order_time         = 2
Carpet::use_buffer_zones          = yes
Carpet::num_integrator_substeps         = 4

#############################################################
# CarpetLib
#############################################################

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

Cactus::terminate                         = "time"
Cactus::cctk_final_time                   = 1000
# Use 0.5 for RK4 and 0.25 for ICN
Time::dtfac                           = 0.5
MethodOfLines::ode_method             = "RK4"
MethodOfLines::MoL_Intermediate_Steps = 4
MethodOfLines::MoL_Num_Scratch_Levels = 1

MethodOfLines::MoL_NaN_Check          = "yes"
MethodOfLines::initial_data_is_crap   = "no"

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

ADMBase::initial_data = "exact"
ADMBase::initial_lapse = "exact"
ADMBase::initial_shift = "exact"

Exact::exact_model = "Minkowski"

Noise::apply_id_noise = yes
Noise::id_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_Theta
"
Noise::Amplitude = 1e-10

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

IO::out_dir                   = $parfile
IO::out_fileinfo              = "none"

CarpetIOScalar::outScalar_every      = 100
CarpetIOScalar::outScalar_reductions = "norm_inf norm2"
CarpetIOScalar::outScalar_vars       = "
  ADMBase::alp
  ADMBase::shift
  ADMBase::metric
  ADMBase::curv
  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_dtlapse
  ML_BSSN::ML_dtshift
  ML_BSSN::ML_shift
  ML_BSSN::ML_Theta
  ML_BSSN::ML_Ham
  ML_BSSN::ML_mom
  ML_BSSN::ML_cons_detg
  ML_BSSN::ML_cons_traceA
  ML_BSSN::ML_cons_Gamma
"

CarpetIOBasic::outInfo_every        = 100
CarpetIOBasic::outInfo_vars         = "ADMBase::alp"

IOASCII::out1D_every          = 1000
IOASCII::out1D_x              = "yes"
IOASCII::out1D_y              = "yes"
IOASCII::out1D_z              = "yes"
IOASCII::out1D_d              = "no"
IOASCII::out1D_vars           = "
  ADMBase::metric
  ADMBase::curv
  ADMBase::alp
  ADMBase::shift
  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_dtlapse
  ML_BSSN::ML_Theta
"

CarpetIOASCII::out_precision = 19