Update example parameter files:

- remove outdated examples
- ensure the vacuum example works
- add matter example

20 files changed, 1875 insertions, 224 deletions

diff --git a/par/ks-mclachlan.par b/par/ks-mclachlan.par
new file mode 100644
index 0000000..7fc652f
--- /dev/null
+++ b/par/ks-mclachlan.par
@@ -0,0 +1,400 @@
+Cactus::cctk_run_title = "Kerr-Schild"
+
+Cactus::cctk_full_warnings         = yes
+Cactus::highlight_warning_messages = no
+
+Cactus::terminate       = "time"
+Cactus::cctk_final_time = 100.0
+
+
+
+ActiveThorns = "IOUtil"
+
+IO::out_dir = $parfile
+
+
+
+ActiveThorns = "AEILocalInterp"
+
+ActiveThorns = "Fortran"
+
+ActiveThorns = "GSL"
+
+ActiveThorns = "GenericFD"
+
+ActiveThorns = "HDF5"
+
+ActiveThorns = "LocalInterp"
+
+ActiveThorns = "LoopControl"
+
+ActiveThorns = "Slab"
+
+
+
+ActiveThorns = "InitBase"
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab"
+
+Carpet::verbose           = no
+Carpet::veryverbose       = no
+Carpet::schedule_barriers = no
+Carpet::storage_verbose   = no
+Carpet::timers_verbose    = no
+CarpetLib::output_bboxes  = no
+
+Carpet::domain_from_coordbase = yes
+Carpet::max_refinement_levels = 10
+
+driver::ghost_size       = 3
+Carpet::use_buffer_zones = yes
+
+Carpet::prolongation_order_space = 5
+Carpet::prolongation_order_time  = 2
+
+Carpet::convergence_level = 0
+
+Carpet::init_fill_timelevels = yes
+
+Carpet::poison_new_timelevels = yes
+CarpetLib::poison_new_memory  = yes
+
+Carpet::output_timers_every      = 5120
+CarpetLib::print_timestats_every = 5120
+CarpetLib::print_memstats_every  = 5120
+
+
+
+ActiveThorns = "NaNChecker"
+
+NaNChecker::check_every     = 512
+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
+"
+
+
+
+ActiveThorns = "Boundary CartGrid3D CoordBase ReflectionSymmetry RotatingSymmetry90 SymBase"
+
+CoordBase::domainsize = "minmax"
+
+CoordBase::xmin =    0.00
+CoordBase::ymin =    0.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
+
+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
+
+CartGrid3D::type = "coordbase"
+
+ReflectionSymmetry::reflection_z   = yes
+ReflectionSymmetry::avoid_origin_z = no
+
+
+
+ActiveThorns = "SphericalSurface"
+
+SphericalSurface::nsurfaces = 1
+SphericalSurface::maxntheta = 39
+SphericalSurface::maxnphi   = 76
+
+SphericalSurface::ntheta      [0] = 39
+SphericalSurface::nphi        [0] = 76
+SphericalSurface::nghoststheta[0] = 2
+SphericalSurface::nghostsphi  [0] = 2
+
+
+
+ActiveThorns = "CarpetMask"
+
+CarpetMask::excluded_surface       [0] = 0
+CarpetMask::excluded_surface_factor[0] = 1.0
+
+
+
+ActiveThorns = "CarpetRegrid2 CarpetTracker"
+
+CarpetTracker::surface[0] = 0
+
+CarpetRegrid2::regrid_every            = 0
+CarpetRegrid2::freeze_unaligned_levels = yes
+CarpetRegrid2::symmetry_rotating90     = yes
+CarpetRegrid2::verbose                 = yes
+
+CarpetRegrid2::num_centres = 1
+
+CarpetRegrid2::num_levels_1 = 6
+CarpetRegrid2::radius_1[1]  = 60.0
+CarpetRegrid2::radius_1[2]  = 30.0
+CarpetRegrid2::radius_1[3]  = 15.0
+CarpetRegrid2::radius_1[4]  =  7.5
+CarpetRegrid2::radius_1[5]  =  3.75
+
+
+
+ActiveThorns = "MoL Time"
+
+MoL::ODE_Method             = "RK4"
+MoL::MoL_Intermediate_Steps = 4
+MoL::MoL_Num_Scratch_Levels = 1
+
+Time::dtfac = 0.40
+
+
+
+ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase"
+
+ADMMacros::spatial_order = 4
+
+
+
+ActiveThorns = "Exact"
+
+ADMBase::initial_data    = "exact"
+ADMBase::initial_lapse   = "exact"
+ADMBase::initial_shift   = "exact"
+ADMBase::initial_dtlapse = "exact"
+ADMBase::initial_dtshift = "exact"
+
+Exact::exact_model = "Kerr/Kerr-Schild"
+Exact::Kerr_KerrSchild__mass      = 1.0
+Exact::Kerr_KerrSchild__spin      = 0.8
+Exact::Kerr_KerrSchild__parabolic = yes
+Exact::Kerr_KerrSchild__power     = 8
+Exact::Kerr_KerrSchild__epsilon   = 1.0
+
+
+
+ActiveThorns = "ML_BSSN ML_BSSN_Helper 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"
+
+ML_BSSN::timelevels = 3
+
+ML_BSSN::harmonicN           = 1      # 1+log
+ML_BSSN::harmonicF           = 2.0    # 1+log
+ML_BSSN::ShiftGammaCoeff     = 0.75
+ML_BSSN::AlphaDriver         = 1.0
+ML_BSSN::BetaDriver          = 1.0
+ML_BSSN::LapseAdvectionCoeff = 0.0
+ML_BSSN::ShiftAdvectionCoeff = 0.0
+
+ML_BSSN::MinimumLapse        = 1.0e-8
+
+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"
+
+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 = "ADMConstraints"
+
+ADMConstraints::constraints_persist    = yes
+ADMConstraints::constraints_timelevels = 3
+
+
+
+ActiveThorns = "AHFinderDirect"
+
+AHFinderDirect::find_every = 1280
+
+AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no
+
+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::output_h_every = 0
+
+AHFinderDirect::N_horizons = 1
+
+AHFinderDirect::initial_guess__coord_sphere__radius[1] = 2.0
+AHFinderDirect::which_surface_to_store_info        [1] = 0
+AHFinderDirect::reset_horizon_after_not_finding    [1] = no
+
+
+
+ActiveThorns = "CarpetIOBasic"
+
+IOBasic::outInfo_every      = 1280
+IOBasic::outInfo_reductions = "norm2"
+IOBasic::outInfo_vars       = "
+        Carpet::physical_time_per_hour
+        ADMConstraints::ham
+"
+
+
+
+ActiveThorns = "CarpetIOScalar"
+
+IOScalar::one_file_per_group = yes
+
+IOScalar::outScalar_every = 1280
+IOScalar::outScalar_vars  = "
+        CarpetReduce::weight
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+        ADMConstraints::hamiltonian
+        ADMConstraints::momentum
+        SphericalSurface::sf_radius
+"
+
+
+
+ActiveThorns = "CarpetIOASCII"
+
+IOASCII::one_file_per_group = yes
+
+IOASCII::output_symmetry_points = no
+IOASCII::out3D_ghosts           = no
+
+IOASCII::out0D_every = 1280
+IOASCII::out0D_vars  = "
+        Carpet::timing
+        CarpetReduce::weight
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+        ADMConstraints::hamiltonian
+        ADMConstraints::momentum
+        SphericalSurface::sf_active
+        SphericalSurface::sf_valid
+        SphericalSurface::sf_info
+        SphericalSurface::sf_radius
+        SphericalSurface::sf_origin
+        SphericalSurface::sf_coordinate_descriptors
+"
+
+IOASCII::out1D_every = 1280
+IOASCII::out1D_vars  = "
+        CarpetReduce::weight
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+        ADMConstraints::hamiltonian
+        ADMConstraints::momentum
+        SphericalSurface::sf_radius
+"
+
+IOASCII::out2D_every = 1280
+IOASCII::out2D_vars  = "
+        SphericalSurface::sf_radius
+"
+
+
+
+Activethorns = "CarpetIOHDF5"
+
+IOHDF5::out_every              = 5120
+IOHDF5::one_file_per_group     = yes
+IOHDF5::output_symmetry_points = no
+IOHDF5::out3D_ghosts           = no
+IOHDF5::compression_level      = 1
+IOHDF5::use_checksums          = yes
+IOHDF5::out_vars               = "
+        CarpetReduce::weight
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+        ADMConstraints::hamiltonian
+        ADMConstraints::momentum
+"
+
+#IOHDF5::checkpoint                  = yes
+#IO::checkpoint_dir                  = $parfile
+#IO::checkpoint_ID                   = yes
+#IO::checkpoint_every_walltime_hours = 6.0
+#IO::checkpoint_on_terminate         = yes
+#
+#IO::recover     = "autoprobe"
+#IO::recover_dir = $parfile
+
+
+
+ActiveThorns = "Formaline"
+
+
+
+ActiveThorns = "TimerReport"
+
+TimerReport::out_every                  = 5120
+TimerReport::out_filename               = "TimerReport"
+TimerReport::output_all_timers_together = yes
+TimerReport::output_all_timers_readable = yes
+TimerReport::n_top_timers               = 20
diff --git a/par/old/bl-mclachlan-ab4.par b/par/old/bl-mclachlan-ab4.par
new file mode 100644
index 0000000..56c7d42
--- /dev/null
+++ b/par/old/bl-mclachlan-ab4.par
@@ -0,0 +1,380 @@
+Cactus::cctk_run_title = "Brill-Lindquist"
+
+Cactus::cctk_full_warnings         = yes
+Cactus::highlight_warning_messages = no
+
+Cactus::terminate       = "time"
+Cactus::cctk_final_time = 100.0
+
+
+
+ActiveThorns = "IOUtil"
+
+IO::out_dir = $parfile
+
+
+
+ActiveThorns = "AEILocalInterp"
+
+ActiveThorns = "Fortran"
+
+ActiveThorns = "LocalInterp"
+
+ActiveThorns = "LoopControl"
+
+ActiveThorns = "GSL"
+
+ActiveThorns = "GenericFD"
+
+ActiveThorns = "HDF5"
+
+ActiveThorns = "Slab"
+
+
+
+ActiveThorns = "InitBase"
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab"
+
+Carpet::verbose           = no
+Carpet::schedule_barriers = no
+Carpet::veryverbose       = no
+Carpet::storage_verbose   = no
+Carpet::timers_verbose    = no
+CarpetLib::output_bboxes  = no
+
+Carpet::domain_from_coordbase = yes
+Carpet::max_refinement_levels = 10
+
+driver::ghost_size       = 3
+Carpet::use_buffer_zones = yes
+
+Carpet::prolongation_order_space = 5
+Carpet::prolongation_order_time  = 2
+
+Carpet::convergence_level = 0
+
+Carpet::init_fill_timelevels = yes
+
+Carpet::poison_new_timelevels = yes
+CarpetLib::poison_new_memory  = yes
+
+Carpet::output_timers_every      = 5120
+CarpetLib::print_timestats_every = 5120
+CarpetLib::print_memstats_every  = 5120
+
+
+
+ActiveThorns = "NaNChecker"
+
+NaNChecker::check_every     = 1 # 512
+NaNChecker::verbose         = "all"
+NaNChecker::action_if_found = "just warn"
+NaNChecker::check_vars      = "
+        ML_BSSN::ML_Gamma
+        ML_BSSN::ML_lapse
+        ML_BSSN::ML_shift
+        ML_BSSN::ML_log_confac
+        ML_BSSN::ML_metric
+        ML_BSSN::ML_trace_curv
+        ML_BSSN::ML_curv
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+"
+
+
+
+ActiveThorns = "Boundary CartGrid3D CoordBase ReflectionSymmetry RotatingSymmetry180 SymBase"
+
+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
+
+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
+
+CartGrid3D::type = "coordbase"
+
+ReflectionSymmetry::reflection_z   = yes
+ReflectionSymmetry::avoid_origin_z = no
+
+
+
+ActiveThorns = "SphericalSurface"
+
+SphericalSurface::nsurfaces = 1
+SphericalSurface::maxntheta = 39
+SphericalSurface::maxnphi   = 76
+
+SphericalSurface::ntheta      [0] = 39
+SphericalSurface::nphi        [0] = 76
+SphericalSurface::nghoststheta[0] = 2
+SphericalSurface::nghostsphi  [0] = 2
+
+
+
+ActiveThorns = "CarpetMask"
+
+CarpetMask::excluded_surface       [0] = 0
+CarpetMask::excluded_surface_factor[0] = 1.0
+
+
+
+ActiveThorns = "CarpetRegrid2 CarpetTracker"
+
+CarpetTracker::surface[0] = 0
+
+CarpetRegrid2::regrid_every            = 0
+CarpetRegrid2::freeze_unaligned_levels = yes
+CarpetRegrid2::symmetry_rotating180    = yes
+CarpetRegrid2::verbose                 = yes
+
+CarpetRegrid2::num_centres = 1
+
+CarpetRegrid2::num_levels_1 =  7
+CarpetRegrid2::radius_1[1]  =  64.0
+CarpetRegrid2::radius_1[2]  =  32.0
+CarpetRegrid2::radius_1[3]  =  16.0
+CarpetRegrid2::radius_1[4]  =   8.0
+CarpetRegrid2::radius_1[5]  =   4.0
+CarpetRegrid2::radius_1[6]  =   2.0
+CarpetRegrid2::radius_1[7]  =   1.0
+
+
+
+ActiveThorns = "MoL Time"
+
+MoL::ODE_Method             = "AB"
+MoL::AB_Type                = "4"
+MoL::MoL_Intermediate_Steps = 1
+
+Time::dtfac = 0.05   # 0.40
+
+
+
+ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase"
+
+ADMMacros::spatial_order = 4
+
+
+
+ActiveThorns = "Exact"
+
+ADMBase::metric_type = "physical"
+
+ADMBase::initial_data    = "exact"
+ADMBase::initial_lapse   = "exact"
+ADMBase::initial_shift   = "exact"
+ADMBase::initial_dtlapse = "zero"
+ADMBase::initial_dtshift = "zero"
+
+Exact::exact_model = "Schwarzschild/BL"
+Exact::Schwarzschild_BL__mass = 1.0
+
+
+
+ActiveThorns = "ML_BSSN ML_BSSN_Helper"
+
+ADMBase::evolution_method       = "ML_BSSN"
+ADMBase::lapse_evolution_method = "ML_BSSN"
+ADMBase::shift_evolution_method = "ML_BSSN"
+
+ML_BSSN::timelevels     = 3
+ML_BSSN::rhs_timelevels = 4
+
+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::my_boundary_condition = "Minkowski"
+
+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"
+
+Dissipation::order = 5
+Dissipation::vars  = "
+        ML_BSSN::ML_Gamma
+        ML_BSSN::ML_lapse
+        ML_BSSN::ML_shift
+        ML_BSSN::ML_log_confac
+        ML_BSSN::ML_metric
+        ML_BSSN::ML_trace_curv
+        ML_BSSN::ML_curv
+"
+
+
+
+ActiveThorns = "ADMConstraints"
+
+ADMConstraints::constraints_persist    = yes
+ADMConstraints::constraints_timelevels = 3
+
+
+
+ActiveThorns = "AHFinderDirect"
+
+AHFinderDirect::find_every = 1280
+
+AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no
+
+AHFinderDirect::move_origins            = yes
+AHFinderDirect::reshape_while_moving    = yes
+AHFinderDirect::predict_origin_movement = yes
+
+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::output_h_every = 0
+
+AHFinderDirect::N_horizons = 1
+
+AHFinderDirect::initial_guess__coord_sphere__radius[1] = 2.0
+AHFinderDirect::which_surface_to_store_info        [1] = 0
+AHFinderDirect::reset_horizon_after_not_finding    [1] = no
+
+
+
+ActiveThorns = "CarpetIOBasic"
+
+IOBasic::outInfo_every      = 1280
+IOBasic::outInfo_reductions = "norm2"
+IOBasic::outInfo_vars       = "
+        Carpet::physical_time_per_hour
+        ADMConstraints::ham
+"
+
+
+
+ActiveThorns = "CarpetIOScalar"
+
+IOScalar::one_file_per_group = yes
+
+IOScalar::outScalar_every = 1280
+IOScalar::outScalar_vars  = "
+        CarpetReduce::weight
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+        ADMConstraints::hamiltonian
+        ADMConstraints::momentum
+"
+
+
+
+ActiveThorns = "CarpetIOASCII"
+
+IOASCII::one_file_per_group = yes
+
+IOASCII::output_symmetry_points = no
+IOASCII::out3D_ghosts           = no
+
+IOASCII::out0D_every = 1280
+IOASCII::out0D_vars  = "
+        Carpet::timing
+        CarpetReduce::weight
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+        ADMConstraints::hamiltonian
+        ADMConstraints::momentum
+"
+
+IOASCII::out1D_every = 1280
+IOASCII::out1D_vars  = "
+        CarpetReduce::weight
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+        ADMConstraints::hamiltonian
+        ADMConstraints::momentum
+"
+
+IOASCII::out2D_every = 1280
+IOASCII::out2D_vars  = "
+        SphericalSurface::sf_radius
+"
+
+
+
+Activethorns = "CarpetIOHDF5"
+
+IOHDF5::out_every              = 5120
+IOHDF5::one_file_per_group     = yes
+IOHDF5::output_symmetry_points = no
+IOHDF5::out3D_ghosts           = no
+IOHDF5::compression_level      = 1
+IOHDF5::use_checksums          = yes
+IOHDF5::out_vars               = "
+        CarpetReduce::weight
+        ADMBase::lapse
+        ADMConstraints::hamiltonian
+"
+
+#IOHDF5::checkpoint                  = yes
+#IO::checkpoint_dir                  = $parfile
+#IO::checkpoint_ID                   = yes
+#IO::checkpoint_every_walltime_hours = 6.0
+#IO::checkpoint_on_terminate         = yes
+#
+#IO::recover     = "autoprobe"
+#IO::recover_dir = $parfile
+
+
+
+ActiveThorns = "Formaline"
+
+
+
+ActiveThorns = "TimerReport"
+
+TimerReport::out_every                  = 5120
+TimerReport::out_filename               = "TimerReport"
+TimerReport::output_all_timers_together = yes
+TimerReport::output_all_timers_readable = yes
+TimerReport::n_top_timers               = 20
diff --git a/par/bl_bssn_1.par b/par/old/bl_bssn_1.par
index 862e994..862e994 100644
--- a/par/bl_bssn_1.par
+++ b/par/old/bl_bssn_1.par
diff --git a/par/bl_bssn_2.par b/par/old/bl_bssn_2.par
index caefed4..caefed4 100644
--- a/par/bl_bssn_2.par
+++ b/par/old/bl_bssn_2.par
diff --git a/par/gw_bssn_1.par b/par/old/gw_bssn_1.par
index 10d9255..10d9255 100644
--- a/par/gw_bssn_1.par
+++ b/par/old/gw_bssn_1.par
diff --git a/par/gw_bssn_2.par b/par/old/gw_bssn_2.par
index ff09c74..ff09c74 100644
--- a/par/gw_bssn_2.par
+++ b/par/old/gw_bssn_2.par
diff --git a/par/head_on.par b/par/old/head_on.par
index e6c326d..e6c326d 100644
--- a/par/head_on.par
+++ b/par/old/head_on.par
diff --git a/par/ks_bssn_0_8th.par b/par/old/ks_bssn_0_8th.par
index fa8b6c6..fa8b6c6 100644
--- a/par/ks_bssn_0_8th.par
+++ b/par/old/ks_bssn_0_8th.par
diff --git a/par/ks_bssn_1.par b/par/old/ks_bssn_1.par
index 4f993f4..4f993f4 100644
--- a/par/ks_bssn_1.par
+++ b/par/old/ks_bssn_1.par
diff --git a/par/ks_bssn_1_8th.par b/par/old/ks_bssn_1_8th.par
index 64b0497..64b0497 100644
--- a/par/ks_bssn_1_8th.par
+++ b/par/old/ks_bssn_1_8th.par
diff --git a/par/ks_bssn_2.par b/par/old/ks_bssn_2.par
index 747dae3..747dae3 100644
--- a/par/ks_bssn_2.par
+++ b/par/old/ks_bssn_2.par
diff --git a/par/ks_bssn_2_8th.par b/par/old/ks_bssn_2_8th.par
index 4c757d8..4c757d8 100644
--- a/par/ks_bssn_2_8th.par
+++ b/par/old/ks_bssn_2_8th.par
diff --git a/par/ks_bssn_3.par b/par/old/ks_bssn_3.par
index 43a96bf..43a96bf 100644
--- a/par/ks_bssn_3.par
+++ b/par/old/ks_bssn_3.par
diff --git a/par/ks_bssn_evol.par b/par/old/ks_bssn_evol.par
index 0c8852c..0c8852c 100644
--- a/par/ks_bssn_evol.par
+++ b/par/old/ks_bssn_evol.par
diff --git a/par/old/ks_bssn_evol_mp.par b/par/old/ks_bssn_evol_mp.par
new file mode 100644
index 0000000..93cd92d
--- /dev/null
+++ b/par/old/ks_bssn_evol_mp.par
@@ -0,0 +1,309 @@
+ActiveThorns = "
+        Boundary
+        CartGrid3D
+        CoordBase
+        IOUtil
+        InitBase
+        MoL
+        SymBase
+	Time
+
+        Coordinates
+        Interpolate
+
+        ADMBase
+        ADMCoupling
+        ADMMacros
+        CoordGauge
+        SpaceMask
+        StaticConformal
+
+        NaNChecker
+
+        Carpet
+        CarpetIOASCII
+        CarpetIOBasic
+        CarpetIOScalar
+        CarpetInterp
+        CarpetLib
+        CarpetReduce
+	CarpetRegrid2
+        LoopControl
+
+        AEILocalInterp
+        AHFinderDirect
+	Dissipation
+        Exact
+        Formaline
+        GSL
+	SphericalSurface
+
+        Slab
+
+        GenericFD
+
+        ML_BSSN_MP_O2
+        ML_BSSN_MP_O2_Helper
+	TmunuBase
+"
+
+# output
+
+IO::out_dir = $parfile
+
+IOBasic::outInfo_every = 1
+IOBasic::outInfo_vars = "
+        ML_BSSN_MP_O2::H
+        ML_BSSN_MP_O2::trKrhs
+"
+
+IOScalar::one_file_per_group = yes
+IOScalar::outScalar_every    = 32
+IOScalar::outScalar_vars     = "
+        grid::coordinates
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ML_BSSN_MP_O2::ML_Ham
+        ML_BSSN_MP_O2::ML_mom
+        ML_BSSN_MP_O2::ML_cons_detg
+        ML_BSSN_MP_O2::ML_cons_Gamma
+        ML_BSSN_MP_O2::ML_cons_traceA
+        ML_BSSN_MP_O2::ML_log_confac
+        ML_BSSN_MP_O2::ML_metric
+        ML_BSSN_MP_O2::ML_trace_curv
+        ML_BSSN_MP_O2::ML_curv
+        ML_BSSN_MP_O2::ML_Gamma
+        ML_BSSN_MP_O2::ML_lapse
+        ML_BSSN_MP_O2::ML_dtlapse
+        ML_BSSN_MP_O2::ML_shift
+        ML_BSSN_MP_O2::ML_dtshift
+        ML_BSSN_MP_O2::ML_log_confacrhs
+        ML_BSSN_MP_O2::ML_metricrhs
+        ML_BSSN_MP_O2::ML_Gammarhs
+        ML_BSSN_MP_O2::ML_trace_curvrhs
+        ML_BSSN_MP_O2::ML_curvrhs
+        ML_BSSN_MP_O2::ML_curvrhs
+        ML_BSSN_MP_O2::ML_lapserhs
+        ML_BSSN_MP_O2::ML_dtlapserhs
+        ML_BSSN_MP_O2::ML_shiftrhs
+        ML_BSSN_MP_O2::ML_dtshiftrhs
+"
+
+IOASCII::one_file_per_group = yes
+IOASCII::out0D_every        = 32
+IOASCII::out0D_vars         = "
+"
+
+IOASCII::out1D_every        = 32
+IOASCII::out1D_vars         = "
+        grid::coordinates
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ML_BSSN_MP_O2::ML_Ham
+        ML_BSSN_MP_O2::ML_mom
+        ML_BSSN_MP_O2::ML_cons_detg
+        ML_BSSN_MP_O2::ML_cons_Gamma
+        ML_BSSN_MP_O2::ML_cons_traceA
+        ML_BSSN_MP_O2::ML_log_confac
+        ML_BSSN_MP_O2::ML_metric
+        ML_BSSN_MP_O2::ML_trace_curv
+        ML_BSSN_MP_O2::ML_curv
+        ML_BSSN_MP_O2::ML_Gamma
+        ML_BSSN_MP_O2::ML_lapse
+        ML_BSSN_MP_O2::ML_dtlapse
+        ML_BSSN_MP_O2::ML_shift
+        ML_BSSN_MP_O2::ML_dtshift
+        ML_BSSN_MP_O2::ML_log_confacrhs
+        ML_BSSN_MP_O2::ML_metricrhs
+        ML_BSSN_MP_O2::ML_Gammarhs
+        ML_BSSN_MP_O2::ML_trace_curvrhs
+        ML_BSSN_MP_O2::ML_curvrhs
+        ML_BSSN_MP_O2::ML_curvrhs
+        ML_BSSN_MP_O2::ML_lapserhs
+        ML_BSSN_MP_O2::ML_dtlapserhs
+        ML_BSSN_MP_O2::ML_shiftrhs
+        ML_BSSN_MP_O2::ML_dtshiftrhs
+"
+
+IOASCII::out2D_every        = 32
+IOASCII::out2D_vars         = "
+        grid::coordinates
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ML_BSSN_MP_O2::ML_Ham
+        ML_BSSN_MP_O2::ML_mom
+        ML_BSSN_MP_O2::ML_cons_detg
+        ML_BSSN_MP_O2::ML_cons_Gamma
+        ML_BSSN_MP_O2::ML_cons_traceA
+        ML_BSSN_MP_O2::ML_log_confac
+        ML_BSSN_MP_O2::ML_metric
+        ML_BSSN_MP_O2::ML_trace_curv
+        ML_BSSN_MP_O2::ML_curv
+        ML_BSSN_MP_O2::ML_Gamma
+        ML_BSSN_MP_O2::ML_lapse
+        ML_BSSN_MP_O2::ML_dtlapse
+        ML_BSSN_MP_O2::ML_shift
+        ML_BSSN_MP_O2::ML_dtshift
+        ML_BSSN_MP_O2::ML_log_confacrhs
+        ML_BSSN_MP_O2::ML_metricrhs
+        ML_BSSN_MP_O2::ML_Gammarhs
+        ML_BSSN_MP_O2::ML_trace_curvrhs
+        ML_BSSN_MP_O2::ML_curvrhs
+        ML_BSSN_MP_O2::ML_curvrhs
+        ML_BSSN_MP_O2::ML_lapserhs
+        ML_BSSN_MP_O2::ML_dtlapserhs
+        ML_BSSN_MP_O2::ML_shiftrhs
+        ML_BSSN_MP_O2::ML_dtshiftrhs
+"
+
+
+#--- driver
+
+Carpet::init_fill_timelevels = yes
+
+Cactus::cctk_itlast = 0 # 1024
+Carpet::use_buffer_zones      = yes
+Carpet::use_tapered_grids     = yes
+
+CarpetLib::print_memstats_every      = 1024
+CarpetLib::print_timestats_every     = 1024
+Carpet::output_timers_every          = 1024
+Carpet::print_timestats_every        = 1024
+
+carpet::verbose = yes
+carpet::veryverbose = no
+
+#--- Carpet prolongation order settings
+
+Carpet::prolongation_order_space                = 3 # 5
+Carpet::prolongation_order_time                 = 1
+
+
+
+# MoL time integration
+
+MoL::ODE_Method             = "RK4"
+MoL::MoL_Intermediate_Steps = 4
+MoL::MoL_Num_Scratch_Levels = 1
+
+Time::timestep_method = "given"
+Time::timestep        = 0.05
+
+#--- grid
+
+driver::ghost_size            = 2 # 3
+grid::type = "MultiPatch"
+
+Coordinates::coordinate_system          = "Thornburg04"
+Coordinates::h_cartesian                = 0.2
+Coordinates::h_radial                   = 0.2
+Coordinates::sphere_inner_radius        = 2.0
+Coordinates::sphere_outer_radius        = 5.0
+Coordinates::n_angular                  = 11
+
+Coordinates::patch_boundary_size        = 2
+Coordinates::additional_overlap_size    = 2
+Coordinates::outer_boundary_size        = 2
+
+Interpolate::interpolator_options       = "order=4"
+
+Carpet::domain_from_multipatch = yes
+
+Carpet::max_refinement_levels = 9
+CarpetRegrid2::num_centres    = 1
+CarpetRegrid2::num_levels_1   = 1 # 3
+CarpetRegrid2::radius_1[1]    = 0.1
+CarpetRegrid2::radius_1[2]    = 0.05
+
+#--- initial data
+
+ADMBase::initial_data      = "exact"
+ADMBase::initial_lapse     = "exact"
+ADMBase::initial_shift     = "exact"
+ADMBase::initial_dtlapse   = "exact"
+ADMBase::initial_dtshift   = "exact"
+
+Exact::exact_model = "Kerr/Kerr-Schild"
+Exact::Kerr_KerrSchild__mass = 1.0
+Exact::Kerr_KerrSchild__spin = 0.0
+#Exact::Kerr_KerrSchild__epsilon = 1.0e-4
+Exact::Kerr_KerrSchild__parabolic = yes
+Exact::Kerr_KerrSchild__power     = 8
+Exact::Kerr_KerrSchild__epsilon   = 1.0
+
+
+#--- ML_BSSN
+
+ADMBase::evolution_method       = "ML_BSSN_MP_O2"
+ADMBase::lapse_evolution_method = "ML_BSSN_MP_O2"
+ADMBase::shift_evolution_method = "ML_BSSN_MP_O2"
+
+ML_BSSN_MP_O2::timelevels = 2
+
+ML_BSSN_MP_O2::my_boundary_condition = "Minkowski"
+
+ML_BSSN_MP_O2::harmonicN       = 1      # 1+log
+ML_BSSN_MP_O2::harmonicF       = 2.0    # 1+log
+ML_BSSN_MP_O2::ShiftGammaCoeff = 0.75
+ML_BSSN_MP_O2::BetaDriver      = 0.5
+
+ML_BSSN_MP_O2::ML_log_confac_bound = "none"
+ML_BSSN_MP_O2::ML_metric_bound     = "none"
+ML_BSSN_MP_O2::ML_Gamma_bound      = "none"
+ML_BSSN_MP_O2::ML_trace_curv_bound = "none"
+ML_BSSN_MP_O2::ML_curv_bound       = "none"
+ML_BSSN_MP_O2::ML_lapse_bound      = "none"
+ML_BSSN_MP_O2::ML_dtlapse_bound    = "none"
+ML_BSSN_MP_O2::ML_shift_bound      = "none"
+ML_BSSN_MP_O2::ML_dtshift_bound    = "none"
+
+#--- Dissipation
+
+Dissipation::order  = 3 # 5
+Dissipation::epsdis = 0.3
+Dissipation::vars   = "
+	ML_BSSN_MP_O2::ML_Gamma
+  	ML_BSSN_MP_O2::ML_lapse
+  	ML_BSSN_MP_O2::ML_shift
+  	ML_BSSN_MP_O2::ML_log_confac
+  	ML_BSSN_MP_O2::ML_metric
+  	ML_BSSN_MP_O2::ML_trace_curv
+  	ML_BSSN_MP_O2::ML_curv
+"
+
+#--- Analysis
+
+nanchecker::check_every                         =       1
+nanchecker::check_vars                          =       "admbase::gxx"
+nanchecker::action_if_found                     =       "terminate"
+
+#--- Horizons
+
+SphericalSurface::nsurfaces = 1
+SphericalSurface::maxntheta = 39
+SphericalSurface::maxnphi   = 76
+
+SphericalSurface::ntheta      [0] = 39
+SphericalSurface::nphi        [0] = 76
+SphericalSurface::nghoststheta[0] = 2
+SphericalSurface::nghostsphi  [0] = 2
+
+AHFinderDirect::find_every = 32
+AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no
+
+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::N_horizons = 1
+
+AHFinderDirect::initial_guess__coord_sphere__radius  [1] = 2.0
+AHFinderDirect::which_surface_to_store_info          [1] = 0
+AHFinderDirect::reset_horizon_after_not_finding      [1] = no
diff --git a/par/old/ks_bssn_evol_mp_hi.par b/par/old/ks_bssn_evol_mp_hi.par
new file mode 100644
index 0000000..d4ed7cc
--- /dev/null
+++ b/par/old/ks_bssn_evol_mp_hi.par
@@ -0,0 +1,309 @@
+ActiveThorns = "
+        Boundary
+        CartGrid3D
+        CoordBase
+        IOUtil
+        InitBase
+        MoL
+        SymBase
+	Time
+
+        Coordinates
+        Interpolate
+
+        ADMBase
+        ADMCoupling
+        ADMMacros
+        CoordGauge
+        SpaceMask
+        StaticConformal
+
+        NaNChecker
+
+        Carpet
+        CarpetIOASCII
+        CarpetIOBasic
+        CarpetIOScalar
+        CarpetInterp
+        CarpetLib
+        CarpetReduce
+	CarpetRegrid2
+        LoopControl
+
+        AEILocalInterp
+        AHFinderDirect
+	Dissipation
+        Exact
+        Formaline
+        GSL
+	SphericalSurface
+
+        Slab
+
+        GenericFD
+
+        ML_BSSN_MP_O2
+        ML_BSSN_MP_O2_Helper
+	TmunuBase
+"
+
+# output
+
+IO::out_dir = $parfile
+
+IOBasic::outInfo_every = 1
+IOBasic::outInfo_vars = "
+        ML_BSSN_MP_O2::H
+        ML_BSSN_MP_O2::trKrhs
+"
+
+IOScalar::one_file_per_group = yes
+IOScalar::outScalar_every    = 64
+IOScalar::outScalar_vars     = "
+        grid::coordinates
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ML_BSSN_MP_O2::ML_Ham
+        ML_BSSN_MP_O2::ML_mom
+        ML_BSSN_MP_O2::ML_cons_detg
+        ML_BSSN_MP_O2::ML_cons_Gamma
+        ML_BSSN_MP_O2::ML_cons_traceA
+        ML_BSSN_MP_O2::ML_log_confac
+        ML_BSSN_MP_O2::ML_metric
+        ML_BSSN_MP_O2::ML_trace_curv
+        ML_BSSN_MP_O2::ML_curv
+        ML_BSSN_MP_O2::ML_Gamma
+        ML_BSSN_MP_O2::ML_lapse
+        ML_BSSN_MP_O2::ML_dtlapse
+        ML_BSSN_MP_O2::ML_shift
+        ML_BSSN_MP_O2::ML_dtshift
+        ML_BSSN_MP_O2::ML_log_confacrhs
+        ML_BSSN_MP_O2::ML_metricrhs
+        ML_BSSN_MP_O2::ML_Gammarhs
+        ML_BSSN_MP_O2::ML_trace_curvrhs
+        ML_BSSN_MP_O2::ML_curvrhs
+        ML_BSSN_MP_O2::ML_curvrhs
+        ML_BSSN_MP_O2::ML_lapserhs
+        ML_BSSN_MP_O2::ML_dtlapserhs
+        ML_BSSN_MP_O2::ML_shiftrhs
+        ML_BSSN_MP_O2::ML_dtshiftrhs
+"
+
+IOASCII::one_file_per_group = yes
+IOASCII::out0D_every        = 64
+IOASCII::out0D_vars         = "
+"
+
+IOASCII::out1D_every        = 64
+IOASCII::out1D_vars         = "
+        grid::coordinates
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ML_BSSN_MP_O2::ML_Ham
+        ML_BSSN_MP_O2::ML_mom
+        ML_BSSN_MP_O2::ML_cons_detg
+        ML_BSSN_MP_O2::ML_cons_Gamma
+        ML_BSSN_MP_O2::ML_cons_traceA
+        ML_BSSN_MP_O2::ML_log_confac
+        ML_BSSN_MP_O2::ML_metric
+        ML_BSSN_MP_O2::ML_trace_curv
+        ML_BSSN_MP_O2::ML_curv
+        ML_BSSN_MP_O2::ML_Gamma
+        ML_BSSN_MP_O2::ML_lapse
+        ML_BSSN_MP_O2::ML_dtlapse
+        ML_BSSN_MP_O2::ML_shift
+        ML_BSSN_MP_O2::ML_dtshift
+        ML_BSSN_MP_O2::ML_log_confacrhs
+        ML_BSSN_MP_O2::ML_metricrhs
+        ML_BSSN_MP_O2::ML_Gammarhs
+        ML_BSSN_MP_O2::ML_trace_curvrhs
+        ML_BSSN_MP_O2::ML_curvrhs
+        ML_BSSN_MP_O2::ML_curvrhs
+        ML_BSSN_MP_O2::ML_lapserhs
+        ML_BSSN_MP_O2::ML_dtlapserhs
+        ML_BSSN_MP_O2::ML_shiftrhs
+        ML_BSSN_MP_O2::ML_dtshiftrhs
+"
+
+IOASCII::out2D_every        = 64
+IOASCII::out2D_vars         = "
+        grid::coordinates
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ML_BSSN_MP_O2::ML_Ham
+        ML_BSSN_MP_O2::ML_mom
+        ML_BSSN_MP_O2::ML_cons_detg
+        ML_BSSN_MP_O2::ML_cons_Gamma
+        ML_BSSN_MP_O2::ML_cons_traceA
+        ML_BSSN_MP_O2::ML_log_confac
+        ML_BSSN_MP_O2::ML_metric
+        ML_BSSN_MP_O2::ML_trace_curv
+        ML_BSSN_MP_O2::ML_curv
+        ML_BSSN_MP_O2::ML_Gamma
+        ML_BSSN_MP_O2::ML_lapse
+        ML_BSSN_MP_O2::ML_dtlapse
+        ML_BSSN_MP_O2::ML_shift
+        ML_BSSN_MP_O2::ML_dtshift
+        ML_BSSN_MP_O2::ML_log_confacrhs
+        ML_BSSN_MP_O2::ML_metricrhs
+        ML_BSSN_MP_O2::ML_Gammarhs
+        ML_BSSN_MP_O2::ML_trace_curvrhs
+        ML_BSSN_MP_O2::ML_curvrhs
+        ML_BSSN_MP_O2::ML_curvrhs
+        ML_BSSN_MP_O2::ML_lapserhs
+        ML_BSSN_MP_O2::ML_dtlapserhs
+        ML_BSSN_MP_O2::ML_shiftrhs
+        ML_BSSN_MP_O2::ML_dtshiftrhs
+"
+
+
+#--- driver
+
+Carpet::init_fill_timelevels = yes
+
+Cactus::cctk_itlast = 0 # 1024
+Carpet::use_buffer_zones      = yes
+Carpet::use_tapered_grids     = yes
+
+CarpetLib::print_memstats_every      = 1024
+CarpetLib::print_timestats_every     = 1024
+Carpet::output_timers_every          = 1024
+Carpet::print_timestats_every        = 1024
+
+carpet::verbose = yes
+carpet::veryverbose = no
+
+#--- Carpet prolongation order settings
+
+Carpet::prolongation_order_space                = 3 # 5
+Carpet::prolongation_order_time                 = 1
+
+
+
+# MoL time integration
+
+MoL::ODE_Method             = "RK4"
+MoL::MoL_Intermediate_Steps = 4
+MoL::MoL_Num_Scratch_Levels = 1
+
+Time::timestep_method = "given"
+Time::timestep        = 0.05
+
+#--- grid
+
+driver::ghost_size            = 2 # 3
+grid::type = "MultiPatch"
+
+Coordinates::coordinate_system          = "Thornburg04"
+Coordinates::h_cartesian                = 0.1
+Coordinates::h_radial                   = 0.1
+Coordinates::sphere_inner_radius        = 2.0
+Coordinates::sphere_outer_radius        = 5.0
+Coordinates::n_angular                  = 21
+
+Coordinates::patch_boundary_size        = 2
+Coordinates::additional_overlap_size    = 2
+Coordinates::outer_boundary_size        = 2
+
+Interpolate::interpolator_options       = "order=4"
+
+Carpet::domain_from_multipatch = yes
+
+Carpet::max_refinement_levels = 9
+CarpetRegrid2::num_centres    = 1
+CarpetRegrid2::num_levels_1   = 1 # 3
+CarpetRegrid2::radius_1[1]    = 0.1
+CarpetRegrid2::radius_1[2]    = 0.05
+
+#--- initial data
+
+ADMBase::initial_data      = "exact"
+ADMBase::initial_lapse     = "exact"
+ADMBase::initial_shift     = "exact"
+ADMBase::initial_dtlapse   = "exact"
+ADMBase::initial_dtshift   = "exact"
+
+Exact::exact_model = "Kerr/Kerr-Schild"
+Exact::Kerr_KerrSchild__mass = 1.0
+Exact::Kerr_KerrSchild__spin = 0.0
+#Exact::Kerr_KerrSchild__epsilon = 1.0e-4
+Exact::Kerr_KerrSchild__parabolic = yes
+Exact::Kerr_KerrSchild__power     = 8
+Exact::Kerr_KerrSchild__epsilon   = 1.0
+
+
+#--- ML_BSSN
+
+ADMBase::evolution_method       = "ML_BSSN_MP_O2"
+ADMBase::lapse_evolution_method = "ML_BSSN_MP_O2"
+ADMBase::shift_evolution_method = "ML_BSSN_MP_O2"
+
+ML_BSSN_MP_O2::timelevels = 2
+
+ML_BSSN_MP_O2::my_boundary_condition = "Minkowski"
+
+ML_BSSN_MP_O2::harmonicN       = 1      # 1+log
+ML_BSSN_MP_O2::harmonicF       = 2.0    # 1+log
+ML_BSSN_MP_O2::ShiftGammaCoeff = 0.75
+ML_BSSN_MP_O2::BetaDriver      = 0.5
+
+ML_BSSN_MP_O2::ML_log_confac_bound = "none"
+ML_BSSN_MP_O2::ML_metric_bound     = "none"
+ML_BSSN_MP_O2::ML_Gamma_bound      = "none"
+ML_BSSN_MP_O2::ML_trace_curv_bound = "none"
+ML_BSSN_MP_O2::ML_curv_bound       = "none"
+ML_BSSN_MP_O2::ML_lapse_bound      = "none"
+ML_BSSN_MP_O2::ML_dtlapse_bound    = "none"
+ML_BSSN_MP_O2::ML_shift_bound      = "none"
+ML_BSSN_MP_O2::ML_dtshift_bound    = "none"
+
+#--- Dissipation
+
+Dissipation::order  = 3 # 5
+Dissipation::epsdis = 0.3
+Dissipation::vars   = "
+	ML_BSSN_MP_O2::ML_Gamma
+  	ML_BSSN_MP_O2::ML_lapse
+  	ML_BSSN_MP_O2::ML_shift
+  	ML_BSSN_MP_O2::ML_log_confac
+  	ML_BSSN_MP_O2::ML_metric
+  	ML_BSSN_MP_O2::ML_trace_curv
+  	ML_BSSN_MP_O2::ML_curv
+"
+
+#--- Analysis
+
+nanchecker::check_every                         =       1
+nanchecker::check_vars                          =       "admbase::gxx"
+nanchecker::action_if_found                     =       "terminate"
+
+#--- Horizons
+
+SphericalSurface::nsurfaces = 1
+SphericalSurface::maxntheta = 39
+SphericalSurface::maxnphi   = 76
+
+SphericalSurface::ntheta      [0] = 39
+SphericalSurface::nphi        [0] = 76
+SphericalSurface::nghoststheta[0] = 2
+SphericalSurface::nghostsphi  [0] = 2
+
+AHFinderDirect::find_every = 64
+AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no
+
+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::N_horizons = 1
+
+AHFinderDirect::initial_guess__coord_sphere__radius  [1] = 2.0
+AHFinderDirect::which_surface_to_store_info          [1] = 0
+AHFinderDirect::reset_horizon_after_not_finding      [1] = no
diff --git a/par/ks_ccatie_1.par b/par/old/ks_ccatie_1.par
index 739fabb..739fabb 100644
--- a/par/ks_ccatie_1.par
+++ b/par/old/ks_ccatie_1.par
diff --git a/par/m_bssn.par b/par/old/m_bssn.par
index f32af64..f32af64 100644
--- a/par/m_bssn.par
+++ b/par/old/m_bssn.par
diff --git a/par/qc0-mclachlan.par b/par/qc0-mclachlan.par
index 5ccdaed..ca2d880 100644
--- a/par/qc0-mclachlan.par
+++ b/par/qc0-mclachlan.par
@@ -4,8 +4,6 @@ Cactus::cctk_full_warnings         = yes
 Cactus::highlight_warning_messages = no
 
 Cactus::terminate       = "time"
-#Cactus::cctk_final_time = 1.0
-#Cactus::cctk_final_time = 10.0
 Cactus::cctk_final_time = 100.0
 
 
@@ -16,32 +14,19 @@ IO::out_dir = $parfile
 
 
 
-ActiveThorns = "ManualTermination"
-
-ManualTermination::max_walltime          = @WALLTIME_HOURS@   # hours
-ManualTermination::on_remaining_walltime = 60                 # minutes
-
-ManualTermination::termination_from_file   = yes
-ManualTermination::create_termination_file = yes
-ManualTermination::termination_file        = "../TERMINATE"
-
-
-
 ActiveThorns = "AEILocalInterp"
 
 ActiveThorns = "BLAS LAPACK"
 
 ActiveThorns = "Fortran"
 
-ActiveThorns = "LocalInterp"
-
 ActiveThorns = "GenericFD"
 
 ActiveThorns = "HDF5"
 
-ActiveThorns = "LoopControl"
+ActiveThorns = "LocalInterp"
 
-#ActiveThorns = "NaNCatcher"
+ActiveThorns = "LoopControl"
 
 ActiveThorns = "Slab"
 
@@ -64,8 +49,10 @@ ActiveThorns = "InitBase"
 ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab"
 
 Carpet::verbose           = no
-Carpet::schedule_barriers = no
 Carpet::veryverbose       = no
+Carpet::schedule_barriers = no
+Carpet::storage_verbose   = no
+Carpet::timers_verbose    = no
 CarpetLib::output_bboxes  = no
 
 Carpet::domain_from_coordbase = yes
@@ -79,29 +66,39 @@ Carpet::prolongation_order_time  = 2
 
 Carpet::convergence_level = 0
 
-Carpet::init_3_timelevels = yes
+Carpet::init_fill_timelevels = yes
 
 Carpet::poison_new_timelevels = yes
 CarpetLib::poison_new_memory  = yes
 
-Carpet::output_timers_every      = 512
-CarpetLib::print_timestats_every = 512
-CarpetLib::print_memstats_every  = 512
+Carpet::output_timers_every      = 5120
+CarpetLib::print_timestats_every = 5120
+CarpetLib::print_memstats_every  = 5120
 
 
 
 ActiveThorns = "NaNChecker"
 
-NaNChecker::check_every     = 512
-NaNChecker::action_if_found = "just warn"
+NaNChecker::check_every     = 1 # 512
+#NaNChecker::verbose         = "all"
+#NaNChecker::action_if_found = "just warn"
+NaNChecker::action_if_found = "terminate"
 NaNChecker::check_vars      = "
-	ML_BSSN::ML_Gamma
-  	ML_BSSN::ML_lapse
-  	ML_BSSN::ML_shift
-  	ML_BSSN::ML_log_confac
-  	ML_BSSN::ML_metric
-  	ML_BSSN::ML_trace_curv
-  	ML_BSSN::ML_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_shift
+        ML_BSSN::ML_dtlapse
+        ML_BSSN::ML_dtshift
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
 "
 
 
@@ -177,7 +174,7 @@ SphericalSurface::nghostsphi  [5] = 2
 
 ActiveThorns = "CarpetMask"
 
-CarpetMask::verbose = yes
+CarpetMask::verbose = no
 
 CarpetMask::excluded_surface       [0] = 0
 CarpetMask::excluded_surface_factor[0] = 1.0
@@ -196,9 +193,10 @@ CarpetTracker::surface[0] = 0
 CarpetTracker::surface[1] = 1
 CarpetTracker::surface[2] = 2
 
-CarpetRegrid2::regrid_every = 128
-
-CarpetRegrid2::symmetry_rotating180 = yes
+CarpetRegrid2::regrid_every            = 128
+CarpetRegrid2::freeze_unaligned_levels = yes
+CarpetRegrid2::symmetry_rotating180    = yes
+CarpetRegrid2::verbose                 = yes
 
 CarpetRegrid2::num_centres = 3
 
@@ -240,7 +238,9 @@ MoL::ODE_Method             = "RK4"
 MoL::MoL_Intermediate_Steps = 4
 MoL::MoL_Num_Scratch_Levels = 1
 
-Time::dtfac = 0.40
+Carpet::time_refinement_factors = "[1, 1, 2, 4, 8, 16, 32, 64, 128, 256]"
+
+Time::dtfac = 0.25
 
 
 
@@ -268,7 +268,6 @@ TwoPunctures::par_P_minus[1] = -0.3331917498
 
 #TODO# TwoPunctures::grid_setup_method = "evaluation"
 
-#TwoPunctures::TP_epsilon = 1.0e-4
 TwoPunctures::TP_epsilon = 1.0e-2
 TwoPunctures::TP_Tiny    = 1.0e-2
 
@@ -276,11 +275,13 @@ TwoPunctures::verbose = yes
 
 
 
-ActiveThorns = "ML_BSSN ML_BSSN_Helper"
+ActiveThorns = "ML_BSSN ML_BSSN_Helper NewRad"
 
-ADMBase::evolution_method       = "ML_BSSN"
-ADMBase::lapse_evolution_method = "ML_BSSN"
-ADMBase::shift_evolution_method = "ML_BSSN"
+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::harmonicN           = 1      # 1+log
 ML_BSSN::harmonicF           = 2.0    # 1+log
@@ -289,17 +290,11 @@ ML_BSSN::BetaDriver          = 1.0
 ML_BSSN::LapseAdvectionCoeff = 1.0
 ML_BSSN::ShiftAdvectionCoeff = 1.0
 
-#ML_BSSN::ML_log_confac_bound = "radiative"
-#ML_BSSN::ML_metric_bound     = "radiative"
-#ML_BSSN::ML_Gamma_bound      = "radiative"
-#ML_BSSN::ML_trace_curv_bound = "radiative"
-#ML_BSSN::ML_curv_bound       = "radiative"
-#ML_BSSN::ML_lapse_bound      = "radiative"
-#ML_BSSN::ML_dtlapse_bound    = "radiative"
-#ML_BSSN::ML_shift_bound      = "radiative"
-#ML_BSSN::ML_dtshift_bound    = "radiative"
+ML_BSSN::MinimumLapse        = 1.0e-8
 
-ML_BSSN::my_boundary_condition = "Minkowski"
+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"
@@ -317,13 +312,15 @@ ActiveThorns = "Dissipation"
 
 Dissipation::order = 5
 Dissipation::vars  = "
-	ML_BSSN::ML_Gamma
-  	ML_BSSN::ML_lapse
-  	ML_BSSN::ML_shift
-  	ML_BSSN::ML_log_confac
-  	ML_BSSN::ML_metric
-  	ML_BSSN::ML_trace_curv
-  	ML_BSSN::ML_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_shift
+        ML_BSSN::ML_dtlapse
+        ML_BSSN::ML_dtshift
 "
 
 
@@ -394,116 +391,13 @@ AHFinderDirect::find_after_individual_time           [6] = 15.0
 
 
 
-ActiveThorns = "IsolatedHorizon"
-
-IsolatedHorizon::verbose              = yes
-IsolatedHorizon::interpolator         = "Lagrange polynomial interpolation"
-IsolatedHorizon::interpolator_options = "order=4"
-IsolatedHorizon::spatial_order        = 6
-
-IsolatedHorizon::num_horizons = 6
-
-IsolatedHorizon::surface_index  [0] = 0
-IsolatedHorizon::companion_index[0] = 4
-
-IsolatedHorizon::surface_index  [1] = 1
-
-IsolatedHorizon::surface_index  [2] = 2
-
-IsolatedHorizon::surface_index  [3] = 3
-IsolatedHorizon::companion_index[3] = 5
-
-IsolatedHorizon::surface_index  [4] = 4
-
-IsolatedHorizon::surface_index  [5] = 5
-
-
-
-ActiveThorns = "ProperDistance ProperTime"
-
-ProperDistance::number_geodesics      = 512
-ProperDistance::direction             = "x"
-ProperDistance::plane                 = "xy"
-ProperDistance::opening_angle         = 180.0
-ProperDistance::step_size             =   0.04
-ProperDistance::horizon_number        = 1
-ProperDistance::interpolation_order   = 3
-ProperDistance::integration_method    = "rk4"
-ProperDistance::eps                   = 0.00001
-ProperDistance::max_proper_distance   = 15
-ProperDistance::calc_every            = 128
-ProperDistance::use_second_horizon    = yes
-ProperDistance::second_horizon_number = 2
-
-
-
-ActiveThorns = "PsiKadelia"
-
-PsiKadelia::psikadelia_persists     = yes
-PsiKadelia::calc_every              = 256   # only on levels 0 and 1
-PsiKadelia::ricci_prolongation_type = "none"
-PsiKadelia::weyl_timelevels         = 3
-PsiKadelia::psif_vec                = "standard-radial"
-
-
-
-ActiveThorns = "ModeDecomposition"
-
-ModeDecomposition::lmax   =  8
-ModeDecomposition::ntheta = 35
-ModeDecomposition::nphi   = 72
-
-ModeDecomposition::num_vars       = 5
-ModeDecomposition::variable_re[0] = "PsiKadelia::psi0re"
-ModeDecomposition::variable_im[0] = "PsiKadelia::psi0im"
-ModeDecomposition::spin_weight[0] = +2
-ModeDecomposition::variable_re[1] = "PsiKadelia::psi1re"
-ModeDecomposition::variable_im[1] = "PsiKadelia::psi1im"
-ModeDecomposition::spin_weight[1] = +1
-ModeDecomposition::variable_re[2] = "PsiKadelia::psi2re"
-ModeDecomposition::variable_im[2] = "PsiKadelia::psi2im"
-ModeDecomposition::spin_weight[2] =  0
-ModeDecomposition::variable_re[3] = "PsiKadelia::psi3re"
-ModeDecomposition::variable_im[3] = "PsiKadelia::psi3im"
-ModeDecomposition::spin_weight[3] = -1
-ModeDecomposition::variable_re[4] = "PsiKadelia::psi4re"
-ModeDecomposition::variable_im[4] = "PsiKadelia::psi4im"
-ModeDecomposition::spin_weight[4] = -2
-
-ModeDecomposition::num_radii = 5
-ModeDecomposition::radius[0] = 30
-ModeDecomposition::radius[1] = 40
-ModeDecomposition::radius[2] = 50
-ModeDecomposition::radius[3] = 60
-ModeDecomposition::radius[4] = 70
-
-
-
-ActiveThorns = "WaveExtract"
-
-WaveExtract::out_every               = 128
-WaveExtract::maximum_detector_number = 5
-WaveExtract::switch_output_format    = 100
-WaveExtract::rsch2_computation       = "average Schwarzschild metric"
-WaveExtract::l_mode                  = 8
-WaveExtract::m_mode                  = 8
-WaveExtract::detector_radius[0]      = 30
-WaveExtract::detector_radius[1]      = 40
-WaveExtract::detector_radius[2]      = 50
-WaveExtract::detector_radius[3]      = 60
-WaveExtract::detector_radius[4]      = 70
-
-
-
 ActiveThorns = "CarpetIOBasic"
 
 IOBasic::outInfo_every      = 128
 IOBasic::outInfo_reductions = "norm2"
 IOBasic::outInfo_vars       = "
-        Carpet::timing
+        Carpet::physical_time_per_hour
         ADMConstraints::ham
-        IsolatedHorizon::ih_spin[0]
-        IsolatedHorizon::ih_radius[0]
 "
 
 
@@ -524,21 +418,6 @@ IOScalar::outScalar_vars  = "
         ADMConstraints::hamiltonian
         ADMConstraints::momentum
         SphericalSurface::sf_radius
-        IsolatedHorizon::ih_shapes
-        IsolatedHorizon::ih_coordinates
-        IsolatedHorizon::ih_tetrad_l
-        IsolatedHorizon::ih_tetrad_n
-        IsolatedHorizon::ih_tetrad_m
-        IsolatedHorizon::ih_newman_penrose
-        IsolatedHorizon::ih_weyl_scalars
-        IsolatedHorizon::ih_twometric
-        IsolatedHorizon::ih_killing_vector
-        IsolatedHorizon::ih_killed_twometric
-        IsolatedHorizon::ih_invariant_coordinates
-        IsolatedHorizon::ih_multipole_moments
-        IsolatedHorizon::ih_3determinant
-        PsiKadelia::WeylComponents
-        ModeDecomposition::modes
 "
 
 
@@ -568,28 +447,6 @@ IOASCII::out0D_vars  = "
         SphericalSurface::sf_radius
         SphericalSurface::sf_origin
         SphericalSurface::sf_coordinate_descriptors
-        IsolatedHorizon::ih_shapes
-        IsolatedHorizon::ih_state
-        IsolatedHorizon::ih_grid_int
-        IsolatedHorizon::ih_grid_real
-        IsolatedHorizon::ih_shapes
-        IsolatedHorizon::ih_coordinates
-        IsolatedHorizon::ih_tetrad_l
-        IsolatedHorizon::ih_tetrad_n
-        IsolatedHorizon::ih_tetrad_m
-        IsolatedHorizon::ih_newman_penrose
-        IsolatedHorizon::ih_weyl_scalars
-        IsolatedHorizon::ih_twometric
-        IsolatedHorizon::ih_killing_vector
-        IsolatedHorizon::ih_killed_twometric
-        IsolatedHorizon::ih_scalars
-        IsolatedHorizon::ih_invariant_coordinates
-        IsolatedHorizon::ih_multipole_moments
-        IsolatedHorizon::ih_3determinant
-        PsiKadelia::WeylComponents
-        ModeDecomposition::modes
-        ProperDistance::pdistance
-        ProperTime::ptime
 "
 
 IOASCII::out1D_every = 128
@@ -604,29 +461,24 @@ IOASCII::out1D_vars  = "
         ADMConstraints::hamiltonian
         ADMConstraints::momentum
         SphericalSurface::sf_radius
-        IsolatedHorizon::ih_shapes
-        IsolatedHorizon::ih_weyl_scalars
-        IsolatedHorizon::ih_killing_vector
-        IsolatedHorizon::ih_killed_twometric
-        IsolatedHorizon::ih_3determinant
-        PsiKadelia::WeylComponents
-        ModeDecomposition::modes
 "
 
 IOASCII::out2D_every = 128
 IOASCII::out2D_vars  = "
         SphericalSurface::sf_radius
-        PsiKadelia::WeylComponents
 "
 
 
 
 Activethorns = "CarpetIOHDF5"
 
-IOHDF5::out_every          = 512
-IOHDF5::one_file_per_group = yes
-IOHDF5::compression_level  = 1
-IOHDF5::out_vars  = "
+IOHDF5::out_every              = 512
+IOHDF5::one_file_per_group     = yes
+IOHDF5::output_symmetry_points = no
+IOHDF5::out3D_ghosts           = no
+IOHDF5::compression_level      = 1
+IOHDF5::use_checksums          = yes
+IOHDF5::out_vars               = "
         CarpetReduce::weight
         ADMBase::metric
         ADMBase::curv
@@ -636,14 +488,13 @@ IOHDF5::out_vars  = "
         ADMBase::dtshift
         ADMConstraints::hamiltonian
         ADMConstraints::momentum
-        PsiKadelia::WeylComponents
 "
 
-IOHDF5::checkpoint          = yes
-IO::checkpoint_dir          = $parfile
-IO::checkpoint_ID           = yes
-IO::checkpoint_every        = 512
-IO::checkpoint_on_terminate = yes
+IOHDF5::checkpoint                  = yes
+IO::checkpoint_dir                  = $parfile
+IO::checkpoint_ID                   = yes
+IO::checkpoint_every_walltime_hours = 6.0
+IO::checkpoint_on_terminate         = yes
 
 IO::recover     = "autoprobe"
 IO::recover_dir = $parfile
@@ -652,14 +503,12 @@ IO::recover_dir = $parfile
 
 ActiveThorns = "Formaline"
 
-#Formaline::publish_level   = 2
-#Formaline::send_as_rdf     = yes
-#Formaline::rdf_hostname[0] = "buran.aei.mpg.de"
-#Formaline::rdf_port    [0] = 24997
-
 
 
 ActiveThorns = "TimerReport"
 
-TimerReport::out_every    = 512
-TimerReport::out_filename = "TimerReport"
+TimerReport::out_every                  = 5120
+TimerReport::out_filename               = "TimerReport"
+TimerReport::output_all_timers_together = yes
+TimerReport::output_all_timers_readable = yes
+TimerReport::n_top_timers               = 20
diff --git a/par/tov-mclachlan.par b/par/tov-mclachlan.par
new file mode 100644
index 0000000..afb169b
--- /dev/null
+++ b/par/tov-mclachlan.par
@@ -0,0 +1,404 @@
+Cactus::cctk_run_title = "Tolman-Oppenheimer-Volkoff"
+
+Cactus::cctk_full_warnings         = yes
+Cactus::highlight_warning_messages = no
+
+Cactus::terminate       = "time"
+Cactus::cctk_final_time = 1500.0
+
+
+
+ActiveThorns = "IOUtil"
+
+IO::out_dir = $parfile
+
+
+
+ActiveThorns = "AEILocalInterp"
+
+ActiveThorns = "Constants"
+
+ActiveThorns = "Fortran"
+
+ActiveThorns = "GSL"
+
+ActiveThorns = "GenericFD"
+
+ActiveThorns = "HDF5"
+
+ActiveThorns = "LocalInterp"
+
+ActiveThorns = "LoopControl"
+
+ActiveThorns = "Slab"
+
+
+
+ActiveThorns = "InitBase"
+
+
+
+ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab"
+
+Carpet::verbose           = no
+Carpet::veryverbose       = no
+Carpet::schedule_barriers = no
+Carpet::storage_verbose   = no
+Carpet::timers_verbose    = no
+CarpetLib::output_bboxes  = no
+
+Carpet::domain_from_coordbase = yes
+Carpet::max_refinement_levels = 10
+
+driver::ghost_size       = 3
+Carpet::use_buffer_zones = yes
+
+Carpet::prolongation_order_space = 5
+Carpet::prolongation_order_time  = 2
+
+Carpet::convergence_level = 0
+
+Carpet::init_fill_timelevels = yes
+
+Carpet::poison_new_timelevels = yes
+CarpetLib::poison_new_memory  = yes
+
+Carpet::output_timers_every      = 5120
+CarpetLib::print_timestats_every = 5120
+CarpetLib::print_memstats_every  = 5120
+
+
+
+ActiveThorns = "NaNChecker"
+
+NaNChecker::check_every     = 1 # 512
+NaNChecker::verbose         = "all"
+NaNChecker::action_if_found = "just warn"
+NaNChecker::check_vars      = "
+        ML_BSSN_M::ML_log_confac
+        ML_BSSN_M::ML_metric
+        ML_BSSN_M::ML_trace_curv
+        ML_BSSN_M::ML_curv
+        ML_BSSN_M::ML_Gamma
+        ML_BSSN_M::ML_lapse
+        ML_BSSN_M::ML_shift
+        ML_BSSN_M::ML_dtlapse
+        ML_BSSN_M::ML_dtshift
+        Whisky::dens
+        Whisky::tau
+        Whisky::scon
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+        HydroBase::rho
+        HydroBase::press
+        HydroBase::eps
+        HydroBase::vel
+"
+
+
+
+ActiveThorns = "Boundary CartGrid3D CoordBase ReflectionSymmetry RotatingSymmetry90 SymBase"
+
+CoordBase::domainsize = "minmax"
+
+CoordBase::xmin =     0.00
+CoordBase::ymin =     0.00
+CoordBase::zmin =     0.00
+CoordBase::xmax = +1800.00
+CoordBase::ymax = +1800.00
+CoordBase::zmax = +1800.00
+CoordBase::dx   =    30.00
+CoordBase::dy   =    30.00
+CoordBase::dz   =    30.00
+
+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
+
+CartGrid3D::type = "coordbase"
+
+ReflectionSymmetry::reflection_z   = yes
+ReflectionSymmetry::avoid_origin_z = no
+
+
+
+ActiveThorns = "SphericalSurface"
+
+
+
+ActiveThorns = "CarpetRegrid2"
+
+CarpetRegrid2::regrid_every            = 0
+CarpetRegrid2::freeze_unaligned_levels = yes
+CarpetRegrid2::symmetry_rotating90     = yes
+CarpetRegrid2::verbose                 = yes
+
+CarpetRegrid2::num_centres = 1
+
+CarpetRegrid2::num_levels_1 = 4
+CarpetRegrid2::radius_1[1]  = 900.0
+CarpetRegrid2::radius_1[2]  = 450.0
+CarpetRegrid2::radius_1[3]  = 225.0
+
+
+
+ActiveThorns = "MoL Time"
+
+MoL::ODE_Method             = "RK3"
+MoL::MoL_Intermediate_Steps = 3
+MoL::MoL_Num_Scratch_Levels = 2
+
+Time::dtfac = 0.25
+
+
+
+ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge HydroBase SpaceMask StaticConformal TmunuBase"
+
+ADMMacros::spatial_order = 4
+
+HydroBase::timelevels        = 3
+HydroBase::prolongation_type = "ENO"
+
+TmunuBase::stress_energy_storage = yes
+TmunuBase::stress_energy_at_RHS  = yes
+TmunuBase::timelevels            = 1
+TmunuBase::prolongation_type     = "none"
+
+SpaceMask::use_mask = yes
+
+
+
+ActiveThorns = "EOS_Base EOS_Polytrope"
+
+EOS_Polytrope::eos_gamma = 1.3333333333
+EOS_Polytrope::eos_k     = 0.4640517
+
+
+
+ActiveThorns = "ML_BSSN_M ML_BSSN_M_Helper NewRad"
+
+ADMBase::evolution_method         = "ML_BSSN_M"
+ADMBase::lapse_evolution_method   = "ML_BSSN_M"
+ADMBase::shift_evolution_method   = "ML_BSSN_M"
+ADMBase::dtlapse_evolution_method = "ML_BSSN_M"
+ADMBase::dtshift_evolution_method = "ML_BSSN_M"
+
+ML_BSSN_M::timelevels = 3
+
+ML_BSSN_M::harmonicN           = 1      # 1+log
+ML_BSSN_M::harmonicF           = 2.0    # 1+log
+ML_BSSN_M::ShiftGammaCoeff     = 0.75
+ML_BSSN_M::AlphaDriver         = 0.0
+ML_BSSN_M::BetaDriver          = 1.0
+ML_BSSN_M::LapseAdvectionCoeff = 1.0
+ML_BSSN_M::ShiftAdvectionCoeff = 1.0
+
+ML_BSSN_M::MinimumLapse        = 1.0e-8
+
+ML_BSSN_M::my_initial_boundary_condition = "extrapolate-gammas"
+ML_BSSN_M::my_rhs_boundary_condition     = "NewRad"
+Boundary::radpower                     = 2
+
+ML_BSSN_M::ML_log_confac_bound = "none"
+ML_BSSN_M::ML_metric_bound     = "none"
+ML_BSSN_M::ML_Gamma_bound      = "none"
+ML_BSSN_M::ML_trace_curv_bound = "none"
+ML_BSSN_M::ML_curv_bound       = "none"
+ML_BSSN_M::ML_lapse_bound      = "none"
+ML_BSSN_M::ML_dtlapse_bound    = "none"
+ML_BSSN_M::ML_shift_bound      = "none"
+ML_BSSN_M::ML_dtshift_bound    = "none"
+
+
+
+ActiveThorns = "Whisky"
+
+HydroBase::evolution_method = "whisky"
+
+Whisky::riemann_solver            = "Marquina"
+Whisky::whisky_eos_type           = "Polytype"
+Whisky::whisky_eos_table          = "2D_Polytrope"
+Whisky::recon_method              = "ppm"
+Whisky::whisky_stencil            = 3
+Whisky::bound                     = "none"
+Whisky::rho_abs_min               = 1.e-12
+Whisky::initial_atmosphere_factor = 0.999
+
+Whisky::ppm_detect        = no
+Whisky::ppm_epsilon       = 0.0
+Whisky::ppm_omega1        = 0.52
+Whisky::ppm_omega2        = 0.0
+Whisky::ppm_epsilon_shock = 0.1
+Whisky::ppm_eta1          = 5.0
+Whisky::ppm_k0            = 0.2
+
+
+
+ActiveThorns = "Whisky_Init_Data Whisky_TOVSolverC"
+
+ADMBase::initial_data    = "tov"
+ADMBase::initial_lapse   = "tov"
+ADMBase::initial_shift   = "tov"
+ADMBase::initial_dtlapse = "zero"
+ADMBase::initial_dtshift = "zero"
+
+WhiskyTOVSolver::TOV_Rho_Central[0] = 1.61930347e-08
+WhiskyTOVSolver::TOV_Gamma[0]       = 1.3333333333
+WhiskyTOVSolver::TOV_K[0]           = 0.4640517
+WhiskyTOVSolver::TOV_Combine_Method = "maximum"
+WhiskyTOVSolver::TOV_Num_Radial     = 200000
+WhiskyTOVSolver::TOV_dr[0]          = 0.05
+
+
+
+ActiveThorns = "Dissipation"
+
+Dissipation::order = 5
+Dissipation::vars  = "
+        ML_BSSN_M::ML_log_confac
+        ML_BSSN_M::ML_metric
+        ML_BSSN_M::ML_trace_curv
+        ML_BSSN_M::ML_curv
+        ML_BSSN_M::ML_Gamma
+        ML_BSSN_M::ML_lapse
+        ML_BSSN_M::ML_shift
+        ML_BSSN_M::ML_dtlapse
+        ML_BSSN_M::ML_dtshift
+"
+
+
+
+ActiveThorns = "ADMConstraints"
+
+ADMConstraints::constraints_persist    = yes
+ADMConstraints::constraints_timelevels = 3
+
+
+
+ActiveThorns = "CarpetIOBasic"
+
+IOBasic::outInfo_every      = 1280
+IOBasic::outInfo_reductions = "norm2"
+IOBasic::outInfo_vars       = "
+        Carpet::physical_time_per_hour
+        ADMConstraints::ham
+        HydroBase::rho
+"
+
+
+
+ActiveThorns = "CarpetIOScalar"
+
+IOScalar::one_file_per_group = yes
+
+IOScalar::outScalar_every = 1280
+IOScalar::outScalar_vars  = "
+        CarpetReduce::weight
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+        HydroBase::rho
+        HydroBase::press
+        HydroBase::eps
+        HydroBase::vel
+        ADMConstraints::hamiltonian
+        ADMConstraints::momentum
+"
+
+
+
+ActiveThorns = "CarpetIOASCII"
+
+IOASCII::one_file_per_group = yes
+
+IOASCII::output_symmetry_points = no
+IOASCII::out3D_ghosts           = no
+
+IOASCII::out0D_every = 1280
+IOASCII::out0D_vars  = "
+        Carpet::timing
+        CarpetReduce::weight
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+        HydroBase::rho
+        HydroBase::press
+        HydroBase::eps
+        HydroBase::vel
+        ADMConstraints::hamiltonian
+        ADMConstraints::momentum
+"
+
+IOASCII::out1D_every = 1280
+IOASCII::out1D_vars  = "
+        CarpetReduce::weight
+        ADMBase::metric
+        ADMBase::curv
+        ADMBase::lapse
+        ADMBase::shift
+        ADMBase::dtlapse
+        ADMBase::dtshift
+        HydroBase::rho
+        HydroBase::press
+        HydroBase::eps
+        HydroBase::vel
+        ADMConstraints::hamiltonian
+        ADMConstraints::momentum
+"
+
+
+
+Activethorns = "CarpetIOHDF5"
+
+IOHDF5::out_every              = 5120
+IOHDF5::one_file_per_group     = yes
+IOHDF5::output_symmetry_points = no
+IOHDF5::out3D_ghosts           = no
+IOHDF5::compression_level      = 1
+IOHDF5::use_checksums          = yes
+IOHDF5::out_vars               = "
+        CarpetReduce::weight
+        ADMBase::lapse
+        HydroBase::rho
+        ADMConstraints::hamiltonian
+"
+
+#IOHDF5::checkpoint                  = yes
+#IO::checkpoint_dir                  = $parfile
+#IO::checkpoint_ID                   = yes
+#IO::checkpoint_every_walltime_hours = 6.0
+#IO::checkpoint_on_terminate         = yes
+#
+#IO::recover     = "autoprobe"
+#IO::recover_dir = $parfile
+
+
+
+ActiveThorns = "Formaline"
+
+
+
+ActiveThorns = "TimerReport"
+
+TimerReport::out_every                  = 5120
+TimerReport::out_filename               = "TimerReport"
+TimerReport::output_all_timers_together = yes
+TimerReport::output_all_timers_readable = yes
+TimerReport::n_top_timers               = 20