Update example parameter files:

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

1 files changed, 73 insertions, 224 deletions

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