Last active
June 26, 2019 09:11
-
-
Save mzilhao/6eadbee1abd6ffb04faa1857e84d72a8 to your computer and use it in GitHub Desktop.
Llama multipatch boundary
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
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 | |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
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 |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#------------------------------------------------------------------------------ | |
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