Dear all,

I have encountered an `tl >= 0 and tl < timelevels' error when recovering checkpoints while using ML_BSSN and helpers.

I am currently running a minimal example in which I merely activate

ML_BSSN ML_BSSN_Helper ML_ADMConstraints

without giving any further options (nor using them for spacetime evolution). The errors occur after I want to to the recovering.

I am sure there is something very basic I am missing in my parameter file. I have included the example in the following.

Thank you for your help and input!

My best from Tel Aviv

Jens

#########################################################################
# Parameter file for the testing of a force-free evolver
# Author: Jens Mahlmann
#########################################################################

#####################
# Thorn specification
#####################
# Kerr setup thorns
#####################
ActiveThorns                            = " Carpet CarpetLib  CarpetInterp 
CarpetReduce "
ActiveThorns                            = " CartGrid3D CoordBase CarpetRegrid2 
SymBase Boundary ReflectionSymmetry  "
ActiveThorns                            = " ADMBase ADMCoupling ADMMacros 
InitBase MoL Time "
ActiveThorns                            = " ML_BSSN ML_BSSN_Helper 
ML_ADMConstraints NewRad"
ActiveThorns                            = " CoordGauge StaticConformal 
LoopControl AEILocalInterp Formaline  "
ActiveThorns                            = " TmunuBase "
ActiveThorns                            = "  SphericalSurface NaNChecker 
SpaceMask "
ActiveThorns                            = " SummationByParts "
ActiveThorns                            = " IOUtil CarpetIOBasic CarpetIOASCII 
CarpetIOScalar CarpetIOHDF5 "
ActiveThorns                            = " CarpetTracker CarpetMask "
ActiveThorns                            = " LocalInterp Slab CarpetSlab Exact  
Fortran "
ActiveThorns                            = "  Dissipation "

#######################
# Simulation management
#######################

Cactus::terminate                               = time
Cactus::cctk_final_time                         = 1000.0

Carpet::output_timers_every                     = 512
Carpet::print_timestats_every                   = 512

CarpetLib::interleave_communications            = yes
CarpetLib::combine_sends                        = yes
CarpetLib::print_memstats_every                 = 1024
CarpetLib::print_timestats_every                = 1024

#####################
# Grid specifications
#####################


driver::ghost_size = 4

CartGrid3D::type                    = "CoordBase"
CartGrid3D::domain                  = "full"


CoordBase::domainsize                       = "minmax"
CoordBase::xmax                             = 128.0 #16.0 #32.0 #40.0 #332.8 
#40.0 #96.0
CoordBase::ymax                             = 128.0 #16.0 #32.0 #40.0 #332.8 
#40.0 #96.0
CoordBase::zmax                             = 128.0 #16.0 #32.0 #40.0 #332.8 
#40.0 #96.0
CoordBase::xmin                             = -128.0 # 0.0 #-32.0 #-40.0 
#-332.8 #-40.0 #-24.0 #-24.0
CoordBase::ymin                             = -128.0 #0.0 #-32.0 #-40.0 #-332.8 
#-40.0 #-24.0 #-24.0
CoordBase::zmin                             = -128.0 # -16.0 #-32.0 #-40.0 
#-332.8 #-40.0 #-24.0 #-24.0
CoordBase::dx                               = 4.0 #1.6 #5.12 #1.6 #2.0 #62.5 
#3.2 #0.8
CoordBase::dy                               = 4.0 #1.6 #5.12 #1.6 #2.0 #62.5 
#3.2 #0.8
CoordBase::dz                               = 4.0 #1.6 #5.12 #1.6 #2.0 #62.5 
#3.2 #0.8
CoordBase::boundary_size_x_lower            = 4
CoordBase::boundary_size_y_lower            = 4
CoordBase::boundary_size_z_lower            = 4
CoordBase::boundary_size_x_upper            = 4
CoordBase::boundary_size_y_upper            = 4
CoordBase::boundary_size_z_upper            = 4
CoordBase::boundary_shiftout_x_lower        = 0
CoordBase::boundary_shiftout_y_lower        = 0
CoordBase::boundary_shiftout_z_lower        = 0
CoordBase::boundary_shiftout_x_upper        = 0
CoordBase::boundary_shiftout_y_upper        = 0
CoordBase::boundary_shiftout_z_upper        = 0
ReflectionSymmetry::reflection_x            = "no"
ReflectionSymmetry::reflection_y            = "no"
ReflectionSymmetry::reflection_z            = "no"
ReflectionSymmetry::avoid_origin_x          = "no"
ReflectionSymmetry::avoid_origin_y          = "no"
ReflectionSymmetry::avoid_origin_z          = "no"

#################
# Grid management
#################

Carpet::domain_from_coordbase                           = "yes"
Carpet::enable_all_storage                              = "yes"
Carpet::use_buffer_zones                                = "yes"
Carpet::schedule_barriers                               = "no"
Carpet::max_refinement_levels                           = 8
Carpet::poison_new_timelevels                           = "yes"
Carpet::check_for_poison                                = "no"
Carpet::poison_value                                    = 113
Carpet::init_3_timelevels                               = "no"
Carpet::init_fill_timelevels                            = "yes"
Carpet::refinement_centering                            = vertex
Carpet::use_tapered_grids                               = no
Carpet::verbose                                         = no
Carpet::veryverbose                                     = no
Carpet::prolongation_order_space                        = 5
Carpet::prolongation_order_time                         = 2

CarpetRegrid2::freeze_unaligned_levels     = "yes"
CarpetRegrid2::regrid_every                = 0
CarpetRegrid2::ensure_proper_nesting       = "yes"
CarpetRegrid2::snap_to_coarse              = "yes"
CarpetRegrid2::num_centres                 = 1
CarpetRegrid2::num_levels_1                = 8
CarpetRegrid2::position_x_1                = 0.0
CarpetRegrid2::position_y_1                = 0.0
CarpetRegrid2::position_z_1                = 0.0
CarpetRegrid2::radius_1[1]                 = 64.0 #4.0 #5.0 #41.6 #0.64 #0.8
CarpetRegrid2::radius_1[2]                 = 32.0 #4.0 #5.0 #41.6 #0.64 #0.8
CarpetRegrid2::radius_x_1[3]               = 24.0 #4.0 #5.0 #41.6 #0.64 #0.8
CarpetRegrid2::radius_y_1[3]               = 24.0 #4.0 #5.0 #41.6 #0.64 #0.8
CarpetRegrid2::radius_z_1[3]               = 16.0 #4.0 #5.0 #41.6 #0.64 #0.8
CarpetRegrid2::radius_x_1[4]               = 16.0 #4.0 #5.0 #41.6 #0.64 #0.8
CarpetRegrid2::radius_y_1[4]               = 16.0 #4.0 #5.0 #41.6 #0.64 #0.8
CarpetRegrid2::radius_z_1[4]               = 8.0 #4.0 #5.0 #41.6 #0.64 #0.8
CarpetRegrid2::radius_1[5]                 = 4.0 #4.0 #5.0 #41.6 #0.64 #0.8
CarpetRegrid2::radius_1[6]                 = 2.0 #4.0 #5.0 #41.6 #0.64 #0.8
CarpetRegrid2::radius_1[7]                 = 1.0 #4.0 #5.0 #41.6 #0.64 #0.8

SpaceMask::use_mask                        = "yes"

#########################################
# Physical quantities evolution procedure
#########################################

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

InitBase::initial_data_setup_method                     = "init_all_levels"

ADMBase::initial_dtlapse                                = "zero"
ADMBase::initial_dtshift                                = "zero"

###############################
# Spacetime evolution procedure
###############################


Carpet::num_integrator_substeps             = 4

ADMMacros::spatial_order                = 4

SummationByParts::order                = 4

Boundary::radpower                              = 2

MoL::ODE_Method                     = "rk4"
MoL::MoL_Intermediate_Steps             = 4
MoL::MoL_Num_Scratch_Levels             = 3
MOL::initial_data_is_crap          = "yes"
Time::dtfac                                     = 0.1


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

IO::out_dir = $parfile

CarpetIOHDF5::one_file_per_group            = yes
CarpetIOHDF5::out3D_every                   = 512
CarpetIOHDF5::out3D_vars                        = "
"

IOScalar::outScalar_every                               = 128
IOScalar::one_file_per_group                    = no
IOScalar::outscalar_reductions                  = "average"
IOScalar::outScalar_vars                                = "
"

IOBasic::outInfo_every              = 1
IOBasic::outInfo_vars               = "Carpet::time_total ADMBase::lapse"


##################
# Runtime analysis
##################
nanchecker::check_every                 = 1
nanchecker::action_if_found             = "terminate"
nanchecker::check_vars                              = "
"

#############
# Checkpoints
#############

Carpet::regrid_during_recovery                  = no

CarpetIOHDF5::checkpoint                    = yes
IO::checkpoint_ID                           = yes
IO::recover                                 = "autoprobe"
IO::checkpoint_every_walltime_hours         = 1
IO::out_proc_every                          = 2
IO::checkpoint_keep                         = 3
IO::checkpoint_on_terminate                 = yes
IO::abort_on_io_errors                      = yes
CarpetIOHDF5::open_one_input_file_at_a_time = yes
_______________________________________________
Users mailing list
[email protected]
http://lists.einsteintoolkit.org/mailman/listinfo/users

Reply via email to