[Users] `tl >= 0 and tl < timelevels' failure in recovering with ML_BSSN

Jens Mahlmann jensmahl at alumni.uv.es
Wed Oct 3 15:51:28 CDT 2018


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

-------------- next part --------------
#########################################################################
# 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


More information about the Users mailing list