Hi all
I'm currently trying to run a BHNS simulation with initial data generated
by Kadath elliptic solver (Fuka branch).
After initial data import, I'm getting an error from THC:
*WARNING level 0 from host nid00462 process 0 while executing schedule bin
HydroBase_Select_Boundaries, routine THC_Core::THC_SelectBC in thorn
THC_Core, file
/lustre/tetyda/home/topolski/Cactus/arrangements/THCCore/THC_Core/src/thc_select_bc.cc:76:
-> Failed to select the BCs*
The simulation uses the following thorns NOT included with the official
release of ET:
HRSCCore
Boost
KadathImporter
THCBase/CPPUtils, THCBase/FDCore ,THCBase/HRSCCore
THCExtra/EOS_Barotropic, THCExtra/EOS_Thermal, THCExtra/EOS_Thermal_Idealgas
Here's the parameter file. If I can provide more info, please let me know
(also, if to upload it to a somewhat more convenient place).
Best regards
Konrad
# The following is a modification of the ETK Gallery par file to utilize
# Kadath initial data for a binary system of BH-NS
################################################################################
# Active thorns
################################################################################
ActiveThorns = "
ADMBase ML_ADMConstraints ADMCoupling ADMMacros
AEILocalInterp AHFinderDirect Boundary Carpet CarpetIOASCII
CarpetIOBasic CarpetIOHDF5 CarpetIOScalar CarpetInterp
CarpetInterp2 CarpetLib CarpetReduce CarpetRegrid2
CarpetTracker CartGrid3D CoordBase Coordinates
CoordinatesSymmetry GlobalDerivative hwloc
IOUtil InitBase Interpolate2 QuasiLocalMeasures
KadathImporter LocalInterp LoopControl MoL
NaNChecker PunctureTracker Slab SpaceMask
SphericalSurface StaticConformal SummationByParts
SymBase SystemStatistics SystemTopology TerminationTrigger
TensorTypes Time TmunuBase Vectors ML_BSSN
ML_BSSN_Helper NewRad GenericFD WeylScal4
Multipole ADMDerivatives CarpetSlab CoordGauge
Dissipation EOS_Barotropic EOS_Thermal EOS_Thermal_Idealgas
HydroBase HRSCCore ML_ADMConstraints THC_Core
TimerReport BLAS LAPACK GSL Boost CPPUtils
FDCore HDF5 MPI Timers CycleClock TGRtensor
zlib Fortran VolumeIntegrals_GRMHD
CarpetMask ReflectionSymmetry
# currently turned off
#GRHydro
#PizzaBase
#PizzaIDBase
#PizzaNumUtils
#NSTracker
#Hydro_Analysis
"
################################################################################
# Grid structure
################################################################################
#------------------------------------------------------------------------------
# Diagnostic parameters:
#------------------------------------------------------------------------------
Carpet::output_timers_every = 0
Carpet::storage_verbose = "no"
Carpet::verbose = "no"
Carpet::veryverbose = "no"
Carpet::grid_structure_filename = "carpet-grid-structure"
Carpet::grid_coordinates_filename = "carpet-grid-coordinates"
CarpetLib::output_bboxes = "no"
CarpetMask::verbose = "no"
CarpetReduce::verbose = "no"
CarpetRegrid2::verbose = "no"
CarpetRegrid2::veryverbose = "no"
CarpetTracker::verbose = "no"
TimerReport::out_every = 4096
TimerReport::out_filename = "TimerReport"
TimerReport::output_all_timers = "yes"
TimerReport::output_all_timers_together = "yes"
TimerReport::output_all_timers_readable = "yes"
TimerReport::n_top_timers = 40
QuasiLocalMeasures::verbose = "no"
SphericalSurface::verbose = "no"
Carpet::domain_from_multipatch = yes
CartGrid3D::type = "multipatch"
CartGrid3D::set_coordinate_ranges_on = "all maps"
Coordinates::coordinate_system = "Thornburg04"
Coordinates::h_cartesian = 1.5574825174818182
Coordinates::h_radial = 1.5574825174818182
Coordinates::sphere_inner_radius = 51.3969230769
Coordinates::sphere_outer_radius = 2194.15912088
Coordinates::n_angular = 28
Driver::ghost_size = 5
Coordinates::patch_boundary_size = 5
Coordinates::additional_overlap_size = 3
Coordinates::outer_boundary_size = 5
CoordinatesSymmetry::reflection_z = yes
CoordinatesSymmetry::stagger = no
Coordinates::symmetry = "+z bitant"
Coordinates::additional_symmetry_size = 1
Coordinates::verbose = no
Time::timestep_method = "given"
Time::timestep = 0.3504335664334091
Carpet::time_refinement_factors = "[1,1,2,4,8,16,32,64,128,256]"
# General Carpet parameters:
Carpet::enable_all_storage = "no"
Carpet::use_buffer_zones = "yes"
Carpet::schedule_barriers = "no"
Carpet::poison_new_timelevels = "yes"
Carpet::check_for_poison = "no"
Carpet::init_3_timelevels = "no"
Carpet::init_fill_timelevels = "yes"
CarpetLib::poison_new_memory = "yes"
CarpetLib::poison_value = 114
CarpetLib::check_bboxes = "no"
CarpetLib::interleave_communications = "yes"
CarpetLib::combine_sends = "yes"
CarpetInterp::tree_search = "yes"
CarpetInterp::check_tree_search = "no"
CarpetRegrid2::freeze_unaligned_levels = "yes"
CarpetRegrid2::freeze_unaligned_parent_levels = "yes"
CarpetRegrid2::ensure_proper_nesting = "yes"
CarpetRegrid2::snap_to_coarse = "yes"
CarpetRegrid2::symmetry_rotating180 = "no"
###############################################################################
# SpaceMask - only needed for GRHydro
###############################################################################
#SpaceMask::use_mask = "yes"
################################################################################
# Mesh refinement
################################################################################
Carpet::max_refinement_levels = 9
CarpetRegrid2::num_centres = 2
CarpetRegrid2::num_levels_1 = 7
CarpetRegrid2::position_x_1 = +17.37737
CarpetRegrid2::radius_1 =
[0,21.267692,10.633846,5.316923,2.658462,1.329231,0.664615]
CarpetRegrid2::num_levels_2 = 7
CarpetRegrid2::position_x_2 = -13.76685
CarpetRegrid2::radius_2 =
[0,17.132308,8.566154,4.283077,2.141538,1.070769,0.535385]
Carpet::prolongation_order_space = 5
Carpet::prolongation_order_time = 2
CarpetRegrid2::regrid_every = 128
################################################################################
# Initial data
################################################################################
HydroBase::timelevels = 3
HydroBase::prolongation_type ="None"
ADMBase::metric_type = "physical"
ADMBase::initial_data = "Kadath"
ADMBase::initial_lapse = "Kadath"
ADMBase::initial_shift = "Kadath"
ADMBase::initial_dtlapse = "Kadath"
ADMBase::initial_dtshift = "Kadath"
HydroBase::initial_hydro = "Kadath"
# Kadath BNS file
KadathImporter::type = "BHNS"
#PizzaIDBase::eos_file = "gam2.polytrope"
# Geometric unit system for initial data, specified by length unit.
# use CACTUS units
#PizzaBase::length_unit = 1476.7161818921163
KadathImporter::filename = "/home/topolski/Cactus/BHNS_par/BHNS.info"
KadathImporter::puncture_lapse = yes
InitBase::initial_data_setup_method = "init_single_level"
#
============================================================================
# Templated hydrodynamics code
#
=============================================================================
HydroBase::evolution_method = "THCode"
HydroBase::initial_entropy = "THCode"
THC_Core::physics = "GRHD"
THC_Core::eos_type = "ideal"
THC_Core::bc_type = "none"
THC_Core::atmo_rho = 1e-13
HRSCCore::scheme = "FV"
HRSCCore::reconstruction = "MP5"
HRSCCore::riemann_solver = "HLLE"
THC_Core::c2a_BH_alp = 0.15
THC_Core::c2a_rho_strict = 2.0e-5
THC_Core::c2a_set_to_nan_on_failure = "no"
THC_Core::c2a_fix_conservatives = "yes"
THC_Core::c2a_kill_on_failure = "no"
EOS_Thermal::evol_eos_name = "IdealGas"
EOS_Thermal_IdealGas::index_n = 1
EOS_Thermal_IdealGas::eps_min = 1e-8
################################################################################
# Evolution and boundary
################################################################################
GenericFD::jacobian_group = "Coordinates::jacobian"
GenericFD::jacobian_derivative_group = "Coordinates::jacobian2"
GenericFD::jacobian_identity_map = 0
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"
ADMBase::lapse_timelevels = 3
ADMBase::shift_timelevels = 3
ADMBase::metric_timelevels = 3
TmunuBase::prolongation_type = "none"
TmunuBase::stress_energy_storage = "yes"
TmunuBase::stress_energy_at_RHS = "yes"
TmunuBase::support_old_CalcTmunu_mechanism = "no"
ML_BSSN::timelevels = 3
ML_BSSN::fdOrder = 8
THC_Core::fd_order = 4
# 1+log slicing requires harmonicN = 1 and harmonicF = 2.0
ML_BSSN::harmonicN = 1
ML_BSSN::harmonicF = 2.0
ML_BSSN::ShiftGammaCoeff = 0.75
ML_BSSN::evolveA = 1
ML_BSSN::evolveB = 1
ML_BSSN::AlphaDriver = 0.0
ML_BSSN::BetaDriver = 0.7 # ~ 1/M
ML_BSSN::advectLapse = 1
ML_BSSN::advectShift = 1
ML_BSSN::MinimumLapse = 1.0e-8
# conformalaMethod = 1 for W, 0 for phi
ML_BSSN::conformalMethod = 1
# We apply dissipation using GlobalDerivaitive so disable it here
ML_BSSN::epsDiss = 0.0
ML_BSSN::dt_lapse_shift_method = "noLapseShiftAdvection"
ML_BSSN::initial_boundary_condition = "extrapolate-gammas"
ML_BSSN::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"
Dissipation::order = 5
Dissipation::epsdis = 0.1
Dissipation::vars = "
ML_BSSN::ML_log_confac
ML_BSSN::ML_metric
ML_BSSN::ML_curv
ML_BSSN::ML_trace_curv
ML_BSSN::ML_Gamma
ML_BSSN::ML_lapse
ML_BSSN::ML_shift
"
################################################################################
# BH and NS tracking (currently using VolumeIntegrals_GRMHD for tracking
the
# neutron star, will switch to NSTracker some time in the future)
# BH is at x>0, NS at x<0
################################################################################
CarpetTracker::surface [0] = 0
CarpetTracker::surface [1] = 1
#PunctureTracker::track [0] = yes
#PunctureTracker::initial_x [0] = -13.76685
#PunctureTracker::which_surface_to_store_info[0] = 0
PunctureTracker::track [0] = yes
PunctureTracker::initial_x [0] = +17.37737
PunctureTracker::which_surface_to_store_info[0] = 1
VolumeIntegrals_GRMHD::verbose = 1
VolumeIntegrals_GRMHD::VolIntegral_out_every = 32
VolumeIntegrals_GRMHD::NumIntegrals = 4
# Integrate the entire volume with an integrand of 1. (used for
testing/validation purposes only).
VolumeIntegrals_GRMHD::Integration_quantity_keyword[1] = "one"
# To compute the center of mass in an integration volume originally
centered at (x,y,z) = (-13.76685,0,0)
#with a coordinate radius of 13.5. Also use the center of mass integral
result to set the ZEROTH AMR center.
VolumeIntegrals_GRMHD::Integration_quantity_keyword[2] = "centerofmass"
VolumeIntegrals_GRMHD::volintegral_sphere__center_x_initial [2] =
-13.76685
VolumeIntegrals_GRMHD::volintegral_inside_sphere__radius [2] =
9.10722
VolumeIntegrals_GRMHD::amr_centre__tracks__volintegral_inside_sphere[2] = 0
# Same as above, except use the integrand=1 (for validation purposes, to
ensure the integration volume is approximately 4/3*pi*13.5^3).
VolumeIntegrals_GRMHD::Integration_quantity_keyword[3] = "one"
VolumeIntegrals_GRMHD::volintegral_sphere__center_x_initial [3] =
-13.76685
VolumeIntegrals_GRMHD::volintegral_inside_sphere__radius [3] =
9.10722
# Perform rest-mass integrals over entire volume.
VolumeIntegrals_GRMHD::Integration_quantity_keyword[4] = "restmass"
VolumeIntegrals_GRMHD::enable_file_output = 1
VolumeIntegrals_GRMHD::outVolIntegral_dir = "volume_integration"
VolumeIntegrals_GRMHD::CoM_integrand_GAMMA_SPEED_LIMIT = 1e4
########################
# Setting up NS Tracker and Hydro_Analysis
#########################
#Hydro_Analysis::Hydro_Analysis_comp_rho_max = "true"
#Hydro_Analysis::Hydro_Analysis_rho_max_loc_only_positive_x = "true"
#Hydro_Analysis::Hydro_Analysis_comp_rho_max_origin_distance = "yes"
#Hydro_Analysis::Hydro_Analysis_average_multiple_maxima_locations = "yes"
#Hydro_Analysis::Hydro_Analysis_interpolator_name = "Lagrange polynomial
interpolation (tensor product)"
#NSTracker::NSTracker_SF_Name = "Lefthand NS"
#NSTracker::NSTracker_SF_Index = 0
#NSTracker::NSTracker_max_distance = 3
#NSTracker::NSTracker_tracked_location =
"Hydro_Analysis::Hydro_Analysis_rho_max_loc" # location of maximum
#position_x[NSTracker_SF_Index] = Hydro_Analysis_rho_max_loc
#position_x[NSTracker_SF_Index_Opposite] = - Hydro_Analysis_rho_max_loc
#NSTracker::NSTracker_tracked_location =
"Hydro_Analysis::Hydro_Analysis_rho_core_center_volume_weighted" # center
of mass
#NSTracker::NSTracker_stop_time = -1
################################################################################
# Spatial finite differencing
################################################################################
SummationByParts::order = 8
# Drop order instead of using upwinded stencils, only for advection
derivatives
SummationByParts::sbp_upwind_deriv = no
SummationByParts::sbp_1st_deriv = yes
SummationByParts::sbp_2nd_deriv = no
SummationByParts::onesided_interpatch_boundaries = no
SummationByParts::onesided_outer_boundaries = yes
SummationByParts::use_dissipation = no
GlobalDerivative::use_dissipation = yes
SummationByParts::scale_with_h = yes
SummationByParts::dissipation_type = "Kreiss-Oliger"
SummationByParts::epsdis = 0.15
# Because dt/dx is half that of the other levels we also need half the
dissipation
GlobalDerivative::epsdis_for_level [0] = 0.075
# Variables for dissipation
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
"
################################################################################
# Time integration
################################################################################
MoL::ODE_Method = "rk4"
MoL::MoL_Intermediate_Steps = 4
MoL::MoL_Num_Scratch_Levels = 1
################################################################################
# Interpolation
################################################################################
# Use 5-th order interpatch interpolation on the Llama grid
Interpolate::interpolator_order = 5
################################################################################
# Psi4 computation
################################################################################
WeylScal4::fdOrder = 8
WeylScal4::calc_scalars = "psis"
WeylScal4::calc_invariants = "always"
################################################################################
# Psi4 mode decomposition
################################################################################
# Radii are chosen to be evenly spaced in 1/r as that is the variable
# extrapolation is performed in
Multipole::nradii = 7
Multipole::radius[0] = 100
Multipole::radius[1] = 115
Multipole::radius[2] = 136
Multipole::radius[3] = 167
Multipole::radius[4] = 214
Multipole::radius[5] = 300
Multipole::radius[6] = 500
Multipole::ntheta = 120
Multipole::nphi = 240
Multipole::variables = "WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i'
name='psi4'}"
Multipole::out_every = 256
Multipole::l_max = 8
Multipole::output_hdf5 = yes
# Disable ASCII output to avoid creating a large number of files
Multipole::output_ascii = no
################################################################################
# Apparent Horizons
################################################################################
AHFinderDirect::N_horizons = 2
AHFinderDirect::find_every = 256
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] = +17.37737
AHFinderDirect::initial_guess__coord_sphere__x_center[1] = +17.37737
AHFinderDirect::initial_guess__coord_sphere__radius [1] = 0.73
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] = -5.53846153846
#AHFinderDirect::initial_guess__coord_sphere__x_center[2] = -5.53846153846
#AHFinderDirect::initial_guess__coord_sphere__radius [2] = 0.535384615385
#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 [2] = 0
AHFinderDirect::find_after_individual [2] = 371902
AHFinderDirect::initial_guess__coord_sphere__x_center[2] = 0
AHFinderDirect::initial_guess__coord_sphere__radius [2] = 1.67
AHFinderDirect::which_surface_to_store_info [2] = 4
AHFinderDirect::set_mask_for_individual_horizon [2] = no
AHFinderDirect::max_allowable_horizon_radius [2] = 6
################################################################################
# Spherical surfaces
################################################################################
SphericalSurface::nsurfaces = 4
SphericalSurface::maxntheta = 66
SphericalSurface::maxnphi = 124
# 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 [3] = 41
SphericalSurface::nphi [3] = 80
SphericalSurface::nghoststheta [3] = 2
SphericalSurface::nghostsphi [3] = 2
################################################################################
# Isolated Horizons
################################################################################
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
QuasiLocalMeasures::output_vtk_every = 0
################################################################################
# Correctness checking
################################################################################
NaNChecker::check_every = 256
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
HydroBase::w_lorentz
ADMBase::curv
ADMBase::metric
ADMBase::lapse
ADMBase::shift
HydroBase::rho
HydroBase::eps
HydroBase::press
HydroBase::vel
"
################################################################################
# Timers
################################################################################
Carpet::output_timer_tree_every = 1024
Carpet::output_initialise_timer_tree = yes
################################################################################
# Output
################################################################################
IO::out_dir = "data"
CarpetIOScalar::all_reductions_in_one_file = "no"
CarpetIOScalar::one_file_per_group = "no"
CarpetIOASCII::one_file_per_group = "no"
CarpetIOHDF5::one_file_per_group = "no"
IOBasic::outInfo_every = 1
IOBasic::outInfo_reductions = "minimum maximum"
IOBasic::outInfo_vars = "
ML_BSSN::ML_Ham
ML_BSSN::ML_log_confac
Carpet::physical_time_per_hour
SystemStatistics::maxrss_mb
SystemStatistics::swap_used_mb
Carpet::physical_time_per_hour
HydroBase::rho
#GRHydro::dens{reductions = 'sum maximum'}
HydroBase::w_lorentz
"
IOScalar::outScalar_every = 256
IOScalar::outScalar_reductions = "minimum maximum average"
IOScalar::outScalar_vars = "
SystemStatistics::process_memory_mb
ML_BSSN::ML_Ham{ reductions = 'maximum' }
"
IOASCII::out0D_every = 256
IOASCII::out0D_vars = "
Carpet::physical_time_per_hour
Carpet::timing
PunctureTracker::pt_loc
QuasiLocalMeasures::qlm_scalars{out_every = 256}
"
IOASCII::out1D_every = 0
IOASCII::out1D_d = no
IOASCII::out1D_x = yes
IOASCII::out1D_y = no
IOASCII::out1D_z = yes
IOASCII::out1D_vars = "
ML_BSSN::ML_log_confac
ML_BSSN::ML_trace_curv
WeylScal4::Psi4r
"
IOASCII::out2D_every = 0
IOASCII::out2D_vars = "
"
IOHDF5::out_every = 0
IOHDF5::out_vars = "
Grid::Coordinates{out_every=1000000000 refinement_levels={0}}
ML_BSSN::ML_log_confac
WeylScal4::Psi4r
WeylScal4::Psi4i
WeylScal4::curvIr{refinement_levels={3 5}}
WeylScal4::curvIi{refinement_levels={3 5}}
WeylScal4::curvJr{refinement_levels={3 5}}
WeylScal4::curvJi{refinement_levels={3 5}}
"
IOHDF5::out2D_every = 64
IOHDF5::out2D_vars = "
Grid::Coordinates{out_every=1000000000 refinement_levels={0}}
ML_BSSN::ML_Ham
ML_BSSN::alpha
ML_BSSN::phi
ML_BSSN::trK
WeylScal4::Psi4r
WeylScal4::Psi4i
HydroBase::rho
HydroBase::press
HydroBase::w_lorentz
"
################################################################################
# Checkpointing and recovery
################################################################################
#CarpetIOHDF5::checkpoint = yes
IOHDF5::checkpoint = "yes"
IO::checkpoint_ID = "yes"
IO::recover = "autoprobe"
IO::checkpoint_every_walltime_hours = 1.5
IO::checkpoint_keep=2
IO::checkpoint_on_terminate = "yes"
IO::out_mode = "onefile"
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 = 0
################################################################################
# Run termination
################################################################################
TerminationTrigger::max_walltime = 48
# Trigger termination 30 minutes before the walltime is reached
TerminationTrigger::on_remaining_walltime = 30
TerminationTrigger::output_remtime_every_minutes = 30
TerminationTrigger::termination_from_file = yes
TerminationTrigger::termination_file = "terminate.txt"
TerminationTrigger::create_termination_file = yes
Cactus::terminate = time
Cactus::cctk_final_time = 1900.0
