Skip to content

Instantly share code, notes, and snippets.

@mzilhao
Last active June 26, 2019 09:11
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save mzilhao/6eadbee1abd6ffb04faa1857e84d72a8 to your computer and use it in GitHub Desktop.
Save mzilhao/6eadbee1abd6ffb04faa1857e84d72a8 to your computer and use it in GitHub Desktop.
Llama multipatch boundary
reflevel = GetRefinementLevel(cctkGH)
map = MultiPatch_GetMap(cctkGH)
write(*,*) 'reflevel = ', reflevel
write(*,*) 'map = ', map
if (reflevel /= 0 .or. map == 0) return
do j = 1, cctk_lsh(2)
do i = 1, cctk_lsh(1)
do k = cctk_lsh(3)-cctk_nghostzones(3)+2, cctk_lsh(3)
if (Sn(i,j,k) == -2) then
write(*,*) 'ijk = ', i,j,k
write(*,*) 'rr = ', rr
write(*,*) 'map = ', map
end if
end do
end do
end do
INFO (Coordinates): Setting up volume form for thornburg04.
INFO (Interpolate2): Setting up interpolation for level 0
INFO (Time): Timestep set to 0.4 (courant_static)
INFO (Time): Timestep set to 0.0224399 (courant_static)
INFO (Time): Timestep set to 0.0224399 (courant_static)
INFO (Time): Timestep set to 0.0224399 (courant_static)
INFO (Time): Timestep set to 0.0224399 (courant_static)
INFO (Time): Timestep set to 0.0224399 (courant_static)
reflevel = 0
map = 1
ijk = 1 1 34
rr = 37.000000000000000
map = 1
ijk = 1 1 35
rr = 38.000000000000000
map = 1
ijk = 2 1 34
rr = 37.000000000000000
map = 1
ijk = 2 1 35
rr = 38.000000000000000
map = 1
ijk = 3 1 34
rr = 37.000000000000000
map = 1
ijk = 3 1 35
rr = 38.000000000000000
map = 1
ijk = 4 1 34
rr = 37.000000000000000
map = 1
ijk = 4 1 35
rr = 38.000000000000000
map = 1
ijk = 5 1 34
rr = 37.000000000000000
map = 1
ijk = 5 1 35
rr = 38.000000000000000
map = 1
ijk = 6 1 34
rr = 37.000000000000000
map = 1
ijk = 6 1 35
rr = 38.000000000000000
map = 1
ijk = 7 1 34
rr = 37.000000000000000
map = 1
ijk = 7 1 35
rr = 38.000000000000000
map = 1
ijk = 8 1 34
rr = 37.000000000000000
map = 1
ijk = 8 1 35
rr = 38.000000000000000
map = 1
ijk = 9 1 34
rr = 37.000000000000000
map = 1
ijk = 9 1 35
rr = 38.000000000000000
map = 1
ijk = 10 1 34
rr = 37.000000000000000
map = 1
ijk = 10 1 35
rr = 38.000000000000000
map = 1
ijk = 11 1 34
rr = 36.999999999999993
map = 1
ijk = 11 1 35
rr = 38.000000000000000
map = 1
ijk = 12 1 34
rr = 37.000000000000007
map = 1
ijk = 12 1 35
rr = 38.000000000000007
map = 1
ijk = 13 1 34
rr = 37.000000000000000
map = 1
ijk = 13 1 35
rr = 37.999999999999993
map = 1
ijk = 14 1 34
rr = 37.000000000000000
map = 1
ijk = 14 1 35
rr = 38.000000000000000
map = 1
ijk = 15 1 34
rr = 37.000000000000007
map = 1
ijk = 15 1 35
rr = 38.000000000000000
map = 1
ijk = 16 1 34
rr = 37.000000000000000
map = 1
ijk = 16 1 35
rr = 38.000000000000000
map = 1
ijk = 17 1 34
rr = 37.000000000000000
map = 1
ijk = 17 1 35
rr = 38.000000000000000
map = 1
ijk = 18 1 34
rr = 37.000000000000000
map = 1
ijk = 18 1 35
rr = 38.000000000000000
map = 1
ijk = 19 1 34
rr = 37.000000000000000
map = 1
ijk = 19 1 35
rr = 38.000000000000000
map = 1
ijk = 20 1 34
rr = 37.000000000000000
map = 1
ijk = 20 1 35
rr = 38.000000000000000
map = 1
ijk = 21 1 34
rr = 37.000000000000000
map = 1
ijk = 21 1 35
rr = 38.000000000000000
map = 1
ijk = 22 1 34
rr = 37.000000000000000
map = 1
ijk = 22 1 35
rr = 38.000000000000000
map = 1
ijk = 23 1 34
rr = 37.000000000000000
map = 1
ijk = 23 1 35
rr = 38.000000000000000
map = 1
ijk = 24 1 34
rr = 37.000000000000000
map = 1
ijk = 24 1 35
rr = 37.999999999999993
map = 1
ijk = 1 2 34
rr = 36.999999999999993
map = 1
ijk = 1 2 35
rr = 37.999999999999993
map = 1
ijk = 2 2 34
rr = 37.000000000000000
map = 1
ijk = 2 2 35
rr = 38.000000000000000
map = 1
ijk = 3 2 34
rr = 36.999999999999993
map = 1
ijk = 3 2 35
rr = 37.999999999999993
map = 1
ijk = 4 2 34
rr = 37.000000000000000
map = 1
ijk = 4 2 35
rr = 38.000000000000000
map = 1
ijk = 5 2 34
rr = 36.999999999999993
map = 1
ijk = 5 2 35
rr = 37.999999999999993
map = 1
ijk = 6 2 34
rr = 37.000000000000000
map = 1
ijk = 6 2 35
rr = 38.000000000000000
map = 1
ijk = 7 2 34
rr = 36.999999999999993
map = 1
ijk = 7 2 35
rr = 37.999999999999993
map = 1
ijk = 8 2 34
rr = 37.000000000000007
map = 1
ijk = 8 2 35
rr = 38.000000000000000
map = 1
ijk = 9 2 34
rr = 37.000000000000000
map = 1
ijk = 9 2 35
rr = 38.000000000000000
map = 1
ijk = 10 2 34
rr = 37.000000000000000
map = 1
ijk = 10 2 35
rr = 38.000000000000000
map = 1
ijk = 11 2 34
rr = 37.000000000000000
map = 1
ijk = 11 2 35
rr = 38.000000000000000
map = 1
ijk = 12 2 34
rr = 37.000000000000000
map = 1
ijk = 12 2 35
rr = 38.000000000000000
map = 1
ijk = 13 2 34
rr = 37.000000000000000
map = 1
ijk = 13 2 35
rr = 38.000000000000000
map = 1
ijk = 14 2 34
rr = 37.000000000000000
map = 1
ijk = 14 2 35
rr = 38.000000000000000
map = 1
ijk = 15 2 34
rr = 37.000000000000007
map = 1
ijk = 15 2 35
rr = 38.000000000000000
map = 1
ijk = 16 2 34
rr = 37.000000000000007
map = 1
ijk = 16 2 35
rr = 38.000000000000000
map = 1
ijk = 17 2 34
rr = 37.000000000000000
map = 1
ijk = 17 2 35
rr = 38.000000000000000
map = 1
ijk = 18 2 34
rr = 37.000000000000000
map = 1
ijk = 18 2 35
rr = 38.000000000000000
map = 1
ijk = 19 2 34
rr = 37.000000000000000
map = 1
ijk = 19 2 35
rr = 38.000000000000000
map = 1
ijk = 20 2 34
rr = 37.000000000000000
map = 1
ijk = 20 2 35
rr = 38.000000000000000
map = 1
ijk = 21 2 34
rr = 37.000000000000000
map = 1
ijk = 21 2 35
rr = 38.000000000000000
map = 1
ijk = 22 2 34
rr = 37.000000000000000
map = 1
ijk = 22 2 35
rr = 38.000000000000000
map = 1
ijk = 23 2 34
rr = 37.000000000000000
map = 1
ijk = 23 2 35
rr = 38.000000000000000
map = 1
ijk = 24 2 34
rr = 37.000000000000000
map = 1
ijk = 24 2 35
rr = 38.000000000000000
map = 1
ijk = 1 3 34
rr = 37.000000000000000
map = 1
ijk = 1 3 35
rr = 38.000000000000000
map = 1
ijk = 2 3 34
rr = 37.000000000000000
map = 1
ijk = 2 3 35
rr = 38.000000000000000
map = 1
ijk = 3 3 34
rr = 37.000000000000000
map = 1
ijk = 3 3 35
rr = 38.000000000000000
map = 1
ijk = 4 3 34
rr = 37.000000000000000
map = 1
ijk = 4 3 35
rr = 38.000000000000000
map = 1
ijk = 5 3 34
rr = 37.000000000000000
map = 1
ijk = 5 3 35
rr = 38.000000000000000
map = 1
ijk = 6 3 34
rr = 37.000000000000000
map = 1
ijk = 6 3 35
rr = 38.000000000000000
map = 1
ijk = 7 3 34
rr = 37.000000000000000
map = 1
ijk = 7 3 35
rr = 38.000000000000000
map = 1
ijk = 8 3 34
rr = 37.000000000000000
map = 1
ijk = 8 3 35
rr = 38.000000000000000
map = 1
ijk = 9 3 34
rr = 36.999999999999993
map = 1
ijk = 9 3 35
rr = 38.000000000000000
map = 1
ijk = 10 3 34
rr = 37.000000000000000
map = 1
ijk = 10 3 35
rr = 38.000000000000000
map = 1
ijk = 11 3 34
rr = 37.000000000000000
map = 1
ijk = 11 3 35
rr = 38.000000000000000
map = 1
ijk = 12 3 34
rr = 37.000000000000000
map = 1
ijk = 12 3 35
rr = 38.000000000000000
map = 1
ijk = 13 3 34
rr = 37.000000000000000
map = 1
ijk = 13 3 35
rr = 38.000000000000000
map = 1
ijk = 14 3 34
rr = 37.000000000000000
map = 1
ijk = 14 3 35
rr = 38.000000000000007
map = 1
ijk = 15 3 34
rr = 37.000000000000000
map = 1
ijk = 15 3 35
rr = 38.000000000000000
map = 1
ijk = 16 3 34
rr = 37.000000000000000
map = 1
ijk = 16 3 35
rr = 38.000000000000000
map = 1
ijk = 17 3 34
rr = 37.000000000000000
map = 1
ijk = 17 3 35
rr = 37.999999999999993
map = 1
ijk = 18 3 34
rr = 37.000000000000007
map = 1
ijk = 18 3 35
rr = 38.000000000000000
map = 1
ijk = 19 3 34
rr = 37.000000000000000
map = 1
ijk = 19 3 35
rr = 38.000000000000000
map = 1
ijk = 20 3 34
rr = 37.000000000000000
map = 1
ijk = 20 3 35
rr = 38.000000000000000
map = 1
ijk = 21 3 34
rr = 37.000000000000000
map = 1
ijk = 21 3 35
rr = 38.000000000000000
map = 1
ijk = 22 3 34
rr = 36.999999999999993
map = 1
ijk = 22 3 35
rr = 37.999999999999993
map = 1
ijk = 23 3 34
rr = 37.000000000000000
map = 1
ijk = 23 3 35
rr = 38.000000000000000
map = 1
ijk = 24 3 34
rr = 37.000000000000000
map = 1
ijk = 24 3 35
rr = 38.000000000000000
map = 1
ijk = 1 4 34
rr = 37.000000000000000
map = 1
ijk = 1 4 35
rr = 38.000000000000000
map = 1
ijk = 2 4 34
rr = 37.000000000000000
map = 1
ijk = 2 4 35
rr = 38.000000000000000
map = 1
ijk = 3 4 34
rr = 37.000000000000000
map = 1
ijk = 3 4 35
rr = 37.999999999999993
map = 1
ijk = 4 4 34
rr = 37.000000000000000
map = 1
ijk = 4 4 35
rr = 38.000000000000000
map = 1
ijk = 5 4 34
rr = 37.000000000000000
map = 1
ijk = 5 4 35
rr = 37.999999999999993
map = 1
ijk = 6 4 34
rr = 37.000000000000000
map = 1
ijk = 6 4 35
rr = 38.000000000000000
map = 1
ijk = 7 4 34
#------------------------------------------------------------------------------
ActiveThorns = "
ADMBase
ADMCoupling
ADMMacros
AEILocalInterp
# AHFinderDirect
BBH_metric
Boundary
Carpet
CarpetInterp
CarpetInterp2
CarpetIOASCII
CarpetIOBasic
CarpetIOHDF5
CarpetIOScalar
CarpetLib
CarpetMask
CarpetReduce
CarpetRegrid2
CarpetSlab
CarpetTracker
CartGrid3D
Coordinates
CoordinatesSymmetry
CoordBase
CoordGauge
Dissipation
# Formaline
Fortran
GenericFD
# GlobalDerivative
GSL
HDF5
InitBase
Interpolate2
IOUtil
# KillingQuantities
# LeanBSSNMoL
IST_IDScalar
IST_PQ
IST_ScalarBase
IST_ScalarEvolve
LocalInterp
LoopControl
# ML_ADMConstraints
MoL
Multipole
NaNChecker
NewRad
# NPScalars
# PunctureTracker
# QuasiLocalMeasures
# ReflectionSymmetry
# RotatingSymmetry90
# RotatingSymmetry180
Slab
SpaceMask
# Sphere
SphericalSurface
StaticConformal
# SummationByParts
SymBase
SystemStatistics
TerminationTrigger
Time
TimerReport
TmunuBase
# Trigger
# TwoPunctures
# Vectors
# WeylScal4
"
#------------------------------------------------------------------------------
# Grid setup
#------------------------------------------------------------------------------
CartGrid3D::type = "multipatch"
CartGrid3D::set_coordinate_ranges_on = "all maps"
Carpet::domain_from_multipatch = yes
Coordinates::coordinate_system = "Thornburg04"
Coordinates::h_radial = 1.0
Coordinates::h_cartesian = 1.0
Coordinates::sphere_inner_radius = 10
Coordinates::sphere_outer_radius = 128.0
Coordinates::n_angular = 28
driver::ghost_size = 3
Coordinates::patch_boundary_size = 3
Coordinates::additional_overlap_size = 3
Coordinates::outer_boundary_size = 3
CoordinatesSymmetry::reflection_z = yes
CoordinatesSymmetry::stagger = no
Coordinates::symmetry = "+z bitant"
Coordinates::additional_symmetry_size = 1
Coordinates::verbose = no
Coordinates::store_volume_form = yes
# Mesh refinement
#------------------------------------------------------------------------------
Carpet::max_refinement_levels = 4
CarpetRegrid2::num_centres = 1
# CarpetRegrid2::num_levels_1 = 4
CarpetRegrid2::num_levels_1 = 1
# CarpetRegrid2::radius_1[ 1] = 4.0
# CarpetRegrid2::radius_1[ 2] = 2.0
# CarpetRegrid2::radius_1[ 3] = 1.0
# CarpetRegrid2::movement_threshold_1 = 0.16
Carpet::use_buffer_zones = yes
Carpet::prolongation_order_space = 5
Carpet::prolongation_order_time = 2
#CarpetRegrid2::freeze_unaligned_levels = yes
#CarpetRegrid2::regrid_every = 32
#CarpetRegrid2::verbose = no
Carpet::grid_structure_filename = "carpet-grid-structure"
Carpet::grid_coordinates_filename = "carpet-grid-coordinates"
# Carpet::time_refinement_factors = "[1]"
Time::dtfac = 0.4
# Initial Data
#------------------------------------------------------------------------------
ADMBase::initial_data = "Cartesian Minkowski"
#ADMBase::initial_lapse = "BBH_metric"
#ADMBase::initial_shift = "BBH_metric"
#ADMBase::initial_dtlapse = "BBH_metric"
#ADMBase::initial_dtshift = "BBH_metric"
ADMBase::lapse_timelevels = 3
ADMBase::shift_timelevels = 3
ADMBase::metric_timelevels = 3
#BBH_metric::m_bh1 = 0.5
#BBH_metric::m_bh2 = 0.5
#BBH_metric::initial_bbh_separation = 20.0002
#BBH_metric::t_shrink_bbh = 1.e6
IST_ScalarBase::mu = 0.0
IST_ScalarBase::initial_scalar = "gaussian"
IST_ScalarEvolve::z_is_radial = yes
IST_IDScalar::amp_g = 1.
IST_IDScalar::amp_c22_g = 0.0
IST_IDScalar::sigma_g = 10.0
IST_IDScalar::R0pert_g = 0.0
IST_IDScalar::x0_g = 0.0001
InitBase::initial_data_setup_method = "init_some_levels"
Carpet::init_fill_timelevels = yes
Carpet::init_3_timelevels = no
# Evolution
#------------------------------------------------------------------------------
TmunuBase::stress_energy_storage = yes
TmunuBase::stress_energy_at_RHS = yes
TmunuBase::timelevels = 1
TmunuBase::prolongation_type = "none"
IST_ScalarBase::evolution_method = "IST_ScalarEvolve"
ADMBase::evolution_method = "static"
IST_ScalarEvolve::use_excision = no
#IST_ScalarEvolve::outer_bound = "radiative_sph"
IST_ScalarEvolve::outer_bound = "reflecting"
# Spatial finite differencing
#------------------------------------------------------------------------------
Dissipation::epsdis = 0.2
Dissipation::order = 5
Dissipation::vars = "
IST_ScalarBase::phi
IST_ScalarBase::Kphi
"
# Integration method
#------------------------------------------------------------------------------
MoL::ODE_Method = "RK4"
MoL::MoL_Intermediate_Steps = 4
MoL::MoL_Num_Scratch_Levels = 1
Carpet::num_integrator_substeps = 4
# Interpolation
#------------------------------------------------------------------------------
CarpetInterp::check_tree_search = no
CarpetInterp::tree_search = yes
# Use 5-th order interpatch interpolation on the Llama grid
Interpolate::interpolator_order = 5
# Spherical surfaces
#------------------------------------------------------------------------------
#CarpetTracker::surface[0] = 0
#CarpetTracker::surface[1] = 1
#SphericalSurface::nsurfaces = 2
#SphericalSurface::maxntheta = 66
#SphericalSurface::maxnphi = 124
# SphericalSurface::verbose = no
#SphericalSurface::verbose = yes
# Horizon 1
#BBH_metric::which_surface_to_store_info[0] = 0
#SphericalSurface::ntheta [0] = 41
#SphericalSurface::nphi [0] = 80
#SphericalSurface::nghoststheta [0] = 2
#SphericalSurface::nghostsphi [0] = 2
#CarpetMask::excluded_surface [0] = 0
#CarpetMask::excluded_surface_factor [0] = 0.8
# Horizon 2
#BBH_metric::which_surface_to_store_info[1] = 1
#SphericalSurface::ntheta [1] = 41
#SphericalSurface::nphi [1] = 80
#SphericalSurface::nghoststheta [1] = 2
#SphericalSurface::nghostsphi [1] = 2
#CarpetMask::excluded_surface [1] = 1
#CarpetMask::excluded_surface_factor [1] = 0.8
# # for AHFinderDirect
# SphericalSurface::ntheta [2] = 41
# SphericalSurface::nphi [2] = 80
# SphericalSurface::nghoststheta [2] = 2
# SphericalSurface::nghostsphi [2] = 2
# SphericalSurface::ntheta [3] = 41
# SphericalSurface::nphi [3] = 80
# SphericalSurface::nghoststheta [3] = 2
# SphericalSurface::nghostsphi [3] = 2
# # Horizons
# #------------------------------------------------------------------------------
# AHFinderDirect::verbose_level = "algorithm highlights"
# # AHFinderDirect::verbose_level = "physics details"
# AHFinderDirect::output_BH_diagnostics = "true"
# AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no
# AHFinderDirect::N_horizons = 2
# AHFinderDirect::find_every = 4
# 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::move_origins = yes
# AHFinderDirect::reshape_while_moving = yes
# AHFinderDirect::predict_origin_movement = yes
# AHFinderDirect::origin_x [1] = -10
# AHFinderDirect::initial_guess__coord_sphere__x_center[1] = -10
# AHFinderDirect::initial_guess__coord_sphere__radius [1] = 0.6
# 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] = "BBH_metric::bh1_loc_x"
# AHFinderDirect::track_origin_source_y [1] = "BBH_metric::bh1_loc_y"
# AHFinderDirect::track_origin_source_z [1] = "BBH_metric::bh1_loc_z"
# AHFinderDirect::max_allowable_horizon_radius [1] = 3
# #AHFinderDirect::dont_find_after_individual_time [1] = 30.0
# AHFinderDirect::origin_x [2] = 10
# AHFinderDirect::initial_guess__coord_sphere__x_center[2] = 10
# AHFinderDirect::initial_guess__coord_sphere__radius [2] = 0.6
# 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] = "BBH_metric::bh2_loc_x"
# AHFinderDirect::track_origin_source_y [2] = "BBH_metric::bh2_loc_y"
# AHFinderDirect::track_origin_source_z [2] = "BBH_metric::bh2_loc_z"
# AHFinderDirect::max_allowable_horizon_radius [2] = 3
# #AHFinderDirect::dont_find_after_individual_time [2] = 30.0
# Wave extraction
#------------------------------------------------------------------------------
Multipole::nradii = 3
Multipole::radius[0] = 5
Multipole::radius[1] = 10
Multipole::radius[2] = 15
#Multipole::radius[3] = 20
Multipole::ntheta = 120
Multipole::nphi = 240
Multipole::variables = "
IST_ScalarBase::phi{sw=0 name='Phi'}
"
Multipole::integration_method = "Simpson"
Multipole::out_every = 64
Multipole::l_max = 4
Multipole::output_hdf5 = no
Multipole::output_ascii = yes
# Check for NaNs
#-------------------------------------------------------------------------------
Carpet::poison_new_timelevels = yes
CarpetLib::poison_new_memory = yes
Carpet::check_for_poison = no
NaNChecker::check_every = 512
NanChecker::check_after = 0
NaNChecker::report_max = 10
# NaNChecker::verbose = "all"
NaNChecker::action_if_found = "terminate"
NaNChecker::out_NaNmask = yes
NaNChecker::check_vars = "
IST_ScalarBase::phi
IST_ScalarBase::Kphi
"
# Timers
#-------------------------------------------------------------------------------
Cactus::cctk_timer_output = "full"
TimerReport::out_every = 5120
TimerReport::n_top_timers = 40
TimerReport::output_all_timers_together = yes
TimerReport::output_all_timers_readable = yes
TimerReport::output_schedule_timers = no
# I/O thorns
#-------------------------------------------------------------------------------
Cactus::cctk_run_title = $parfile
IO::out_dir = $parfile
IOScalar::one_file_per_group = yes
IOASCII::one_file_per_group = yes
IOHDF5::use_checksums = no
IOHDF5::one_file_per_group = no
IOBasic::outInfo_every = 1
IOBasic::outInfo_reductions = "minimum maximum"
IOBasic::outInfo_vars = "
Carpet::physical_time_per_hour
# ADMBase::lapse
IST_ScalarBase::phi
SystemStatistics::maxrss_mb
"
# for scalar reductions of 3D grid functions
IOScalar::outScalar_every = 1
IOScalar::outScalar_reductions = "minimum maximum average norm2"
IOScalar::outScalar_vars = "
# ADMBase::lapse
# ML_ADMConstraints::ML_mom
# ML_ADMConstraints::ML_Ham
IST_ScalarBase::phi
"
# output just at one point (0D)
IOASCII::out0D_point_x = 0.0
IOASCII::out0D_point_y = 0.0
IOASCII::out0D_every = 1
IOASCII::out0D_vars = "
IST_ScalarBase::phi
IST_ScalarBase::Kphi
# SphericalSurface::sf_info
# SphericalSurface::sf_origin
#BBH_metric::bh1_loc
#BBH_metric::bh2_loc
# ML_ADMConstraints::ML_mom
# ML_ADMConstraints::ML_Ham
IST_PQ::komarenergy
IST_PQ::energy
"
IOASCII::output_symmetry_points = no
IOASCII::out3D_ghosts = no
# 1D text output
IOASCII::out1D_every = 1
IOASCII::out1D_d = no
IOASCII::out1D_x = yes
IOASCII::out1D_y = yes
IOASCII::out1D_z = yes
IOASCII::out1D_vars = "
#ADMBase::lapse
#ADMBase::shift
#ADMBase::metric
#ADMBase::curv
IST_ScalarBase::phi
IST_ScalarBase::Kphi
# SpaceMask::space_mask
"
# 2D HDF5 output
IOHDF5::out2D_every = 256
IOHDF5::out2D_xy = yes
IOHDF5::out2D_xz = no
IOHDF5::out2D_yz = no
IOHDF5::out2D_vars = "
Grid::Coordinates{out_every=1000000000 refinement_levels={0}}
IST_ScalarBase::phi
# IST_ScalarBase::Kphi
# ADMBase::lapse
# ADMBase::shift
# ADMBase::metric
# ADMBase::curv
# ML_ADMConstraints::ML_mom
# ML_ADMConstraints::ML_Ham
# SpaceMask::space_mask
"
# 3D HDF5 output
IOHDF5::out_every = 256
IOHDF5::out_vars = "
#ADMBase::lapse
#ScalarBase::phi1
#ScalarBase::Kphi1
Grid::Coordinates{out_every=1000000000 refinement_levels={0}}
IST_ScalarBase::phi
#Interpolate::Sn
"
Carpet::verbose = no
Carpet::veryverbose = no
Carpet::schedule_barriers = no
Carpet::storage_verbose = no
CarpetLib::output_bboxes = no
Cactus::cctk_full_warnings = yes
Cactus::highlight_warning_messages = no
# Checkpointing and recovery
#-------------------------------------------------------------------------------
CarpetIOHDF5::checkpoint = yes
IO::checkpoint_dir = "../checkpoints"
IO::checkpoint_ID = no
IO::checkpoint_every_walltime_hours = 23
IO::checkpoint_on_terminate = yes
IO::out_proc_every = 2
IO::checkpoint_keep = 1
IO::recover = "autoprobe"
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 = 72 # hours
TerminationTrigger::on_remaining_walltime = 60 # minutes
TerminationTrigger::output_remtime_every_minutes = 60
TerminationTrigger::termination_from_file = yes
TerminationTrigger::termination_file = "terminate.txt"
TerminationTrigger::create_termination_file = yes
Cactus::terminate = "time"
Cactus::cctk_final_time = 1
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment