Add some sample CCZ4 parameter files.

3 files changed, 1006 insertions, 0 deletions

diff --git a/par/CCZ4/bbh_ccz4.par b/par/CCZ4/bbh_ccz4.par
new file mode 100644
index 0000000..e009fd5
--- /dev/null
+++ b/par/CCZ4/bbh_ccz4.par
@@ -0,0 +1,493 @@
+################################################################################
+  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
diff --git a/par/CCZ4/gw1d_ccz4.par b/par/CCZ4/gw1d_ccz4.par
new file mode 100644
index 0000000..e74f951
--- /dev/null
+++ b/par/CCZ4/gw1d_ccz4.par
@@ -0,0 +1,283 @@
+Cactus::cctk_run_title = "CCZ4 1D Gauge Wave"
+
+ActiveThorns = "
+  ADMBase
+  ADMCoupling
+  ADMMacros
+  Boundary
+  Carpet
+  CarpetInterp
+  CarpetInterp2
+  CarpetIOASCII
+  CarpetIOBasic
+  CarpetIOHDF5
+  CarpetIOScalar
+  CarpetLib
+  CarpetReduce
+  CarpetRegrid
+  CarpetSlab
+  CartGrid3d
+  CoordBase
+  CoordGauge
+  Dissipation
+  Exact
+  GenericFD
+  IOUtil
+  LoopControl
+  ML_BSSN
+  ML_BSSN_Helper
+  MoL
+  NanChecker
+  Periodic
+  Slab
+  SpaceMask
+  SphericalSurface
+  StaticConformal
+  SummationByParts
+  SymBase
+  Time
+  TimerReport
+  TMuNuBase
+"
+
+
+#############################################################
+# 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          = 0.0
+ML_BSSN::dampk1              = 1.0
+ML_BSSN::harmonicN           = 2     
+ML_BSSN::harmonicF           = 1.0   
+ML_BSSN::ShiftGammaCoeff     = 0.0 # beta to zero
+ML_BSSN::ShiftBCoeff         = 0.0 # B to zero
+ML_BSSN::BetaDriver          = 0.0 # damp in beta
+ML_BSSN::LapseACoeff         = 0.0 # A to zero
+ML_BSSN::LapseAdvectionCoeff = 0.0 # no advection terms
+ML_BSSN::ShiftAdvectionCoeff = 0.0 # no advection terms
+
+ML_BSSN::MinimumLapse        = 1.0e-8
+ML_BSSN::conformalMethod     = 1 # 1 for W, 0 for phi
+ML_BSSN::formulation         = 1 # 1 for CCZ4
+
+#############################################################
+# Boundary conditions
+#############################################################
+
+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"
+
+#############################################################
+# Spatial derivatives
+#############################################################
+
+SummationByParts::order                              = 4
+SummationByParts::sbp_upwind_deriv = no # Drop order instead, only for
+                                        # advection derivatives
+SummationByParts::sbp_1st_deriv                      = no
+SummationByParts::sbp_2nd_deriv                      = no
+SummationByParts::onesided_interpatch_boundaries     = no
+SummationByParts::onesided_outer_boundaries          = no
+SummationByParts::use_dissipation                    = yes
+SummationByParts::scale_with_h                       = no
+SummationByParts::dissipation_type                   = "Kreiss-Oliger"
+SummationByParts::epsdis                             = 0.1
+SummationByParts::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
+"
+
+#############################################################
+# Grid
+#############################################################
+
+CoordBase::domainsize                   = minmax
+
+CoordBase::boundary_size_x_lower        = 3
+CoordBase::boundary_size_y_lower        = 3
+CoordBase::boundary_size_z_lower        = 3
+CoordBase::boundary_shiftout_x_lower    = 1
+CoordBase::boundary_shiftout_y_lower    = 1
+CoordBase::boundary_shiftout_z_lower    = 1
+
+CoordBase::boundary_size_x_upper        = 3
+CoordBase::boundary_size_y_upper        = 3
+CoordBase::boundary_size_z_upper        = 3
+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.005
+CoordBase::dy                           = 0.1
+CoordBase::dz                           = 0.1
+
+#############################################################
+# 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                      = 3
+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/gauge wave"
+Exact::Minkowski_gauge_wave__what_fn   = "sin"
+Exact::Minkowski_gauge_wave__amplitude = 0.1
+Exact::Minkowski_gauge_wave__omega     = 1
+Exact::Minkowski_gauge_wave__lambda    = 1
+Exact::Minkowski_gauge_wave__diagonal  = "no"
+
+################################################################################
+# Constraints
+################################################################################
+
+################################################################################
+# Timers
+################################################################################
+
+Cactus::cctk_timer_output               = full
+TimerReport::out_every                  = 4000
+TimerReport::n_top_timers               = 40
+TimerReport::out_filename               = "TimerReport"
+TimerReport::output_all_timers          = yes
+
+#############################################################
+# 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        = 1
+CarpetIOBasic::outInfo_vars         = "
+  ADMBase::alp
+"
+
+IOASCII::out1D_every          = 4000
+IOASCII::out1D_x              = "yes"
+IOASCII::out1D_y              = "yes"
+IOASCII::out1D_z              = "yes"
+IOASCII::out1D_d              = "no"
+IOASCII::out1D_vars           = "
+  ADMBase::alp
+  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_Theta
+  ML_BSSN::ML_Ham
+  ML_BSSN::ML_mom
+  ML_BSSN::ML_cons_detg
+  ML_BSSN::ML_cons_traceA
+  ML_BSSN::ML_cons_Gamma
+"
+
+CarpetIOASCII::out_precision = 19
+
diff --git a/par/CCZ4/rbst_ccz4.par b/par/CCZ4/rbst_ccz4.par
new file mode 100644
index 0000000..3801ba4
--- /dev/null
+++ b/par/CCZ4/rbst_ccz4.par
@@ -0,0 +1,230 @@
+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
+