[Users] McLachlan and Cartoon2D conflicting?
Federico Maria Guercilena
fguercilena at theorie.ikp.physik.tu-darmstadt.de
Mon Sep 14 13:10:11 CDT 2020
Hello everyone,
I have found a very weird issue in trying to run 2D axysimmetric
simulations with Cartoon2D and McLachlan, I'm really puzzled on what
could be the cause and how to address it.
In brief: in running a test TOV simulation with the above setup, the
constraints have NaNs around the z-axis although all hydro and spacetime
variables are fine, and the simulation crashes after a few iterations
due to the NaNs propagating. However, this behavior is
resolution-dependent (nothing strange at low resolution). Also Lean
doesn't show any problems at any resolution, and the simulations run
just fine in this case.
More in detail, the setup that show this problem (the parfile is
attached) is just a test TOV. It is a 2D axisymmetric run, using
Cartoon2D. The spacetime is handled by ML_CCZ4, and I'm using David
Radice's THC code for the hydro evolution (which strictly speaking is
not part of the ET, but as you'll see the hydro is not at fault here).
The initial data is handled by TOVSolver. When running at the resolution
labelled as "mid" in the parfile, everything goes as planned, the ID is
set up properly (in particular constraints violations are negligible),
and the evolution proceeds nicely (the amplitude of the star
oscillations seems a bit higher than it would be in a 3D simulation, but
that might just be noise from the Cartoon2D boundary conditions).
However when using the "high" or "fine" resolutions in the parfile
(respectively 1.25 and 1.25^2=1.5625 times higher than mid), there are
NaNs in the constraints, especially the Hamiltonian constraint, around
the z-axis, in a sort of band (see the attached plot). The NaNs are on
all reflevels, and the width of the band on the xz-plane varies,
seemingly not just because the resolution varies between reflevels, but
the NaNs band is sometimes a 1 or 2 points wider or narrower. Trying to
evolve this data the NaNs propagate, the con2prim kills the hydro and
progressively I end up with a "vacuum" simulation. ML_BSSN behaves
similarly.
All variables (primitives, conserved, spacetime, stress-energy tensor)
are properly set as far as I can tell (no NaNs or weird values, and
correct symmetry), even ADMMass::ADMMass_VolumeMass_GF. The fact that
this last one is ok gave me the hunch that the problem is in the
spacetime solver, and in fact switching to Lean did the trick: proper
evolution at all resolutions. This also tells me the hydro code is not
to blame, nor the initial data code (tried with a different, custom
written one, no change).
For various reasons (Lean is slower and doesn't implement CCZ4) I'd like
to stick to ML_CCZ4. At first I thought that the way ML_CCZ4 handles
symmetries and boundaries doesn't play well with Cartoon2D (some
ordering problem), so I tried to modify its scheduler to make it as
close to Lean as possible, but no change. It might be a problem with my
grid, but then I don't see why Lean should work. I'm starting to think
there's an issue with McLachlan internals (... the derivative
stencils?), but than it should show up in 3D runs as well.
In brief, I'm baffled.
I'd be grateful for any insights, and of course I'll open a ticket if
need be. But I'm still harboring a glimpse of hope that this is just a
trivial mistake on my part...
Thank you very much,
Federico Guercilena
PS Unrelated: in a couple of previous mails I wrote about the GCC10
issues, and how on my laptop I couldn't compile despite setting the
necessary flags. That is still the case and I still have no idea why,
but thanks to the latest efforts (mainly by Roland, as I understand?)
the ET has been made "GCC10-proof" and that is not an issue anymore. I
just wanted to say thanks.
--
Dr. Federico Maria Guercilena
Technische Universität Darmstadt
Institut für Kernphysik
Theoriezentrum
S2|11
Schlossgartenstraße 12
64289 Darmstadt
Room 302
+49 6151 16 21562
fguercilena at theorie.ikp.physik.tu-darmstadt.de
-------------- next part --------------
#!/usr/bin/env python2
from math import *
import sys
import re
from string import Template
from subprocess import call
import os
###########################
# Outer boundary
outer_boundary = 5120.
###########################
# Resolution
# The resolutions increase in steps of 1.25, a constant factor,
# which should make it easy to perform a convergence study.
# The resolution is set by setting directly dx on the coarsest level.
# The values [12.8, 10.24, 8.192] result in resolutions of
# [0.2, 0.16, 0.128] on level 6
try:
resolution = str("@RESOLUTION@")
except:
resolution = "mid"
if resolution == "@"+"RESOLUTION"+"@":
resolution = "mid"
if resolution=="mid":
dx = 12.8
elif resolution=="high":
dx = 10.24
elif resolution=="fine":
dx = 8.192
else:
print("Unknown RESOLUTION keyword")
sys.exit()
###########################
# CIDRFILE
try :
cidr_file = str("@CIDRF@")
except :
cidr_file = "nosuchthing"
if cidr_file == "@"+"CIDRF"+"@":
cidr_file = "nosuchthing"
###########################
# EVOLVE
try :
evolve = str("@EVOLVE@")
except :
evolve = "no"
if evolve == "@"+"EVOLVE"+"@":
evolve = "no"
if evolve == "no":
end = "iteration"
elif evolve == "yes":
end = "never"
###########################
# betaDriver
try:
betaDriver = float(str("@BDRIVER@"))
except:
betaDriver = 1.0
if betaDriver == "@"+"BDRIVER"+"@":
betaDriver = 1.0
lines = """
ActiveThorns = "
ADMBase
ADMCoupling
ADMMass
ADMMacros
AEILocalInterp
Boundary
Carpet
CarpetInterp
CarpetIOASCII
CarpetIOBasic
CarpetIOHDF5
CarpetIOScalar
CarpetLib
CarpetReduce
CarpetRegrid2
CarpetSlab
CartGrid3D
Cartoon2D
CoordBase
CoordGauge
Constants
hwloc
EOS_Barotropic
EOS_Thermal
EOS_Thermal_IdealGas
Fortran
GenericFD
HydroBase
HRSCCore
InitBase
IOUtil
LoopControl
ML_CCZ4
ML_CCZ4_Helper
MoL
NaNChecker
NewRad
Nice
PizzaBase
PizzaNumUtils
ReflectionSymmetry
Slab
SpaceMask
SphericalSurface
StaticConformal
SymBase
SystemStatistics
SystemTopology
TerminationTrigger
THC_Core
THC_Tracer
Time
TimerReport
TmunuBase
TOVSolver
WatchDog
"
# =============================================================================
# General
# =============================================================================
Cactus::cctk_run_title = "CIDR"
Cactus::cctk_show_banners = "no"
Cactus::recovery_mode = "strict"
Cactus::schedule_sort_mode = "ascending"
Cactus::terminate = "$end"
Cactus::cctk_itlast = 0
TerminationTrigger::max_walltime = @WALLTIME_HOURS@
TerminationTrigger::on_remaining_walltime = 10
TerminationTrigger::create_termination_file = "yes"
TerminationTrigger::termination_from_file = "yes"
TerminationTrigger::termination_file = "../TERMINATE"
TerminationTrigger::output_remtime_every_minutes = 120
TimerReport::out_every = 0
TimerReport::output_all_timers_readable = "yes"
Nice::Nice_nice = 0
# =============================================================================
# Grid
# =============================================================================
Grid::avoid_origin = "no"
Grid::domain = "full"
Grid::type = "coordbase"
ReflectionSymmetry::reflection_x = "no"
ReflectionSymmetry::reflection_y = "no"
ReflectionSymmetry::reflection_z = "yes"
ReflectionSymmetry::avoid_origin_x = "no"
ReflectionSymmetry::avoid_origin_y = "no"
ReflectionSymmetry::avoid_origin_z = "no"
CoordBase::xmin = 0.
CoordBase::xmax = +$outer_boundary
CoordBase::ymin = 0.
CoordBase::ymax = 0.
CoordBase::zmin = 0.
CoordBase::zmax = +$outer_boundary
CoordBase::spacing = "gridspacing"
CoordBase::dx = $dx
CoordBase::dy = $dx
CoordBase::dz = $dx
CoordBase::boundary_size_x_lower = 3
CoordBase::boundary_size_x_upper = 3
CoordBase::boundary_shiftout_x_lower = 1
CoordBase::boundary_shiftout_x_upper = 1
CoordBase::boundary_staggered_x_lower = "no"
CoordBase::boundary_staggered_x_upper = "no"
CoordBase::boundary_size_y_lower = 3
CoordBase::boundary_size_y_upper = 3
CoordBase::boundary_shiftout_y_lower = 1
CoordBase::boundary_shiftout_y_upper = 1
CoordBase::boundary_staggered_y_lower = "no"
CoordBase::boundary_staggered_y_upper = "no"
CoordBase::boundary_size_z_lower = 3
CoordBase::boundary_size_z_upper = 3
CoordBase::boundary_shiftout_z_lower = 1
CoordBase::boundary_shiftout_z_upper = 1
CoordBase::boundary_staggered_z_lower = "no"
CoordBase::boundary_staggered_z_upper = "no"
Driver::ghost_size = 3
Driver::ghost_size_x = 3
Driver::ghost_size_y = 3
Driver::ghost_size_z = 3
Cartoon2D::cartoon_active = "yes"
Cartoon2D::order = 3
Cartoon2D::eno_order = 3
Carpet::processor_topology = "automatic"
Carpet::no_split_direction = 1
Carpet::domain_from_coordbase = "yes"
InitBase::initial_data_setup_method = "init_some_levels"
LoopControl::settle_after_iteration = 0
Carpet::max_refinement_levels = 7
Carpet::prolongation_order_space = 3
Carpet::prolongation_order_time = 2
CarpetLib::use_higher_order_restriction = "no"
Carpet::use_buffer_zones = "no"
Carpet::enable_all_storage = "no"
Carpet::init_fill_timelevels = "yes"
Carpet::time_refinement_factors = "[1,1,2,4,8,16,32,64,128,256,512]"
Carpet::grid_coordinates_filename = "grid.carpet"
CarpetLib::poison_new_memory = "yes"
CarpetRegrid2::num_centres = 1
CarpetRegrid2::regrid_every = -1
CarpetRegrid2::ensure_proper_nesting = "yes"
CarpetRegrid2::snap_to_coarse = "yes"
CarpetRegrid2::freeze_unaligned_levels = "yes"
CarpetRegrid2::freeze_unaligned_parent_levels = "yes"
CarpetRegrid2::active_1 = "yes"
CarpetRegrid2::num_levels_1 = 7
CarpetRegrid2::position_x_1 = 0.0
CarpetRegrid2::radius_1[1] = 512.0
CarpetRegrid2::radius_1[2] = 256.0
CarpetRegrid2::radius_1[3] = 128.0
CarpetRegrid2::radius_1[4] = 64.0
CarpetRegrid2::radius_1[5] = 32.0
CarpetRegrid2::radius_1[6] = 16.0
CarpetRegrid2::symmetry_rotating180 = "no"
CarpetRegrid2::verbose = "no"
Carpet::print_grid_info = "no"
NaNChecker::check_every = 4096
NaNChecker::check_vars = "
HydroBase::rho
HydroBase::w_lorentz
"
NaNChecker::action_if_found = "abort"
# =============================================================================
# Time integration
# =============================================================================
Carpet::num_integrator_substeps = 3
MoL::ode_method = "RK3"
MoL::MoL_Intermediate_Steps = 3
MoL::MoL_Num_Scratch_Levels = 0
MoL::verbose = "register"
MoL::disable_prolongation = "no"
HydroBase::timelevels = 3
Time::timestep_method = "courant_static"
Time::dtfac = 0.15
# =============================================================================
# Initial data
# =============================================================================
ADMBase::initial_data = "tov"
ADMBase::initial_lapse = "tov"
ADMBase::initial_shift = "tov"
ADMBase::initial_dtlapse = "zero"
ADMBase::initial_dtshift = "zero"
HydroBase::initial_hydro = "tov"
HydroBase::initial_entropy = "zero"
HydroBase::initial_temperature = "zero"
HydroBase::initial_Y_e = "one"
TOVSolver::TOV_Num_TOVs = 1
TOVSolver::TOV_Enforce_Interpolation = "yes"
TOVSolver::TOV_Rho_Central[0] = 0.00128
TOVSolver::TOV_Fast_Interpolation = "no"
TOVSolver::TOV_Combine_Method = "maximum"
TOVSolver::TOV_save_to_datafile = ""
# =============================================================================
# Templated hydrodynamics code
# =============================================================================
HydroBase::evolution_method = "THCode"
HydroBase::temperature_evolution_method = "none"
HydroBase::entropy_evolution_method = "none"
HydroBase::Y_e_evolution_method = "none"
HydroBase::Abar_evolution_method = "none"
THC_Core::verbose = "no"
THC_Core::eos_type = "ideal"
THC_Core::physics = "GRHD"
TmunuBase::prolongation_type = "none"
TmunuBase::stress_energy_storage = "yes"
TmunuBase::stress_energy_at_RHS = "yes"
TmunuBase::support_old_CalcTmunu_mechanism = "no"
THC_Core::bc_type = "none"
HRSCCore::scheme = "FV"
HRSCCore::pplim = "yes"
HRSCCore::reconstruction = "MP5"
HRSCCore::flux_split = "LLF"
HRSCCore::riemann_solver = "HLLE"
HRSCCore::system_split = "components"
HRSCCore::speed_eps = 0.05
HRSCCore::refluxing = "no"
THC_Core::atmo_rho = 1e-16
THC_Core::atmo_eps = 1e-8
THC_Core::atmo_tolerance = 1e-3
THC_Core::atmo_I_know_what_I_am_doing = "yes"
THC_Core::c2a_BH_alp = 0.15
THC_Core::c2a_rho_strict = 2e-5
THC_Core::c2a_set_to_nan_on_failure = "no"
THC_Core::c2a_fix_conservatives = "yes"
THC_Core::c2a_kill_on_failure = "no"
THC_Core::c2a_max_num_msg = 0
EOS_Thermal::evol_eos_name = "IdealGas"
EOS_Thermal_IdealGas::index_n = 1.0
EOS_Thermal_IdealGas::rho_max = 10.0
EOS_Thermal_IdealGas::rho_min = 0.0
EOS_Thermal_IdealGas::eps_max = 100.0
EOS_Thermal_IdealGas::eps_min = 1e-8
# =============================================================================
# Spacetime evolution
# =============================================================================
ADMBase::evolution_method = "ML_CCZ4"
ADMBase::lapse_evolution_method = "ML_CCZ4"
ADMBase::shift_evolution_method = "ML_CCZ4"
ADMBase::dtlapse_evolution_method = "ML_CCZ4"
ADMBase::dtshift_evolution_method = "ML_CCZ4"
ML_CCZ4::timelevels = 3
ML_CCZ4::harmonicN = 1.0 # 1+log
ML_CCZ4::harmonicF = 2.0 # 1+log
ML_CCZ4::alphaDriver = 0.0 # not more than 0.5
ML_CCZ4::betaDriver = $betaDriver # ~1/M (eta)
ML_CCZ4::advectLapse = 1
ML_CCZ4::advectShift = 1
ML_CCZ4::evolveA = 0
ML_CCZ4::evolveB = 1
ML_CCZ4::shiftGammaCoeff = 0.75
ML_CCZ4::shiftFormulation = 0 # Gamma driver
ML_CCZ4::fixAdvectionTerms = 0
ML_CCZ4::dampk1 = 0.02
ML_CCZ4::dampk2 = 0.0
ML_CCZ4::GammaShift = 0.5
ML_CCZ4::MinimumLapse = 1.0e-8
ML_CCZ4::conformalMethod = 1 # 1 for W
ML_CCZ4::dt_lapse_shift_method = "noLapseShiftAdvection"
ML_CCZ4::initial_boundary_condition = "extrapolate-gammas"
ML_CCZ4::rhs_boundary_condition = "scalar"
Boundary::radpower = 2
ML_CCZ4::fdOrder = 4
THC_Core::fd_order = 4
ML_CCZ4::epsDiss = 0.2
ML_CCZ4::ML_CCZ4_EvolutionAnalysisInit_calc_every = 1
ML_CCZ4::ML_CCZ4_EvolutionAnalysisInit_calc_offset = 2
ML_CCZ4::ML_CCZ4_EvolutionAnalysisInterior_calc_every = 1
ML_CCZ4::ML_CCZ4_EvolutionAnalysisInterior_calc_offset = 2
ML_CCZ4::ML_CCZ4_ConstraintsEverywhere_calc_every = 1
ML_CCZ4::ML_CCZ4_ConstraintsEverywhere_calc_offset = 2
ML_CCZ4::ML_CCZ4_ConstraintsInterior_calc_every = 1
ML_CCZ4::ML_CCZ4_ConstraintsInterior_calc_offset = 0
ML_CCZ4::ML_log_confac_bound = "none"
ML_CCZ4::ML_log_confac_bound_speed = 1.0
ML_CCZ4::ML_log_confac_bound_limit = 0.0
ML_CCZ4::ML_log_confac_bound_scalar = 0.0
ML_CCZ4::ML_metric_bound = "none"
ML_CCZ4::ML_metric_bound_speed = 1.0
ML_CCZ4::ML_metric_bound_limit = 0.0
ML_CCZ4::ML_metric_bound_scalar = 0.0
ML_CCZ4::ML_Gamma_bound = "none"
ML_CCZ4::ML_Gamma_bound_speed = 1.0
ML_CCZ4::ML_Gamma_bound_limit = 0.0
ML_CCZ4::ML_Gamma_bound_scalar = 0.0
ML_CCZ4::ML_trace_curv_bound = "none"
ML_CCZ4::ML_trace_curv_bound_speed = 1.0
ML_CCZ4::ML_trace_curv_bound_limit = 0.0
ML_CCZ4::ML_trace_curv_bound_scalar = 0.0
ML_CCZ4::ML_curv_bound = "none"
ML_CCZ4::ML_curv_bound_speed = 1.0
ML_CCZ4::ML_curv_bound_limit = 0.0
ML_CCZ4::ML_curv_bound_scalar = 0.0
ML_CCZ4::ML_Theta_bound = "none"
ML_CCZ4::ML_Theta_bound_speed = 1.0
ML_CCZ4::ML_Theta_bound_limit = 0.0
ML_CCZ4::ML_Theta_bound_scalar = 0.0
ML_CCZ4::ML_lapse_bound = "none"
ML_CCZ4::ML_lapse_bound_speed = 1.0
ML_CCZ4::ML_lapse_bound_limit = 0.0
ML_CCZ4::ML_lapse_bound_scalar = 0.0
ML_CCZ4::ML_dtlapse_bound = "none"
ML_CCZ4::ML_dtlapse_bound_speed = 1.0
ML_CCZ4::ML_dtlapse_bound_limit = 0.0
ML_CCZ4::ML_dtlapse_bound_scalar = 0.0
ML_CCZ4::ML_shift_bound = "none"
ML_CCZ4::ML_shift_bound_speed = 1.0
ML_CCZ4::ML_shift_bound_limit = 0.0
ML_CCZ4::ML_shift_bound_scalar = 0.0
ML_CCZ4::ML_dtshift_bound = "none"
ML_CCZ4::ML_dtshift_bound_speed = 1.0
ML_CCZ4::ML_dtshift_bound_limit = 0.0
ML_CCZ4::ML_dtshift_bound_scalar = 0.0
ML_CCZ4::phi_bound = "skip"
ML_CCZ4::phi_bound_speed = 1.0
ML_CCZ4::phi_bound_limit = 0.0
ML_CCZ4::phi_bound_scalar = 0.0
ML_CCZ4::gt11_bound = "skip"
ML_CCZ4::gt11_bound_speed = 1.0
ML_CCZ4::gt11_bound_limit = 0.0
ML_CCZ4::gt11_bound_scalar = 0.0
ML_CCZ4::gt12_bound = "skip"
ML_CCZ4::gt12_bound_speed = 1.0
ML_CCZ4::gt12_bound_limit = 0.0
ML_CCZ4::gt12_bound_scalar = 0.0
ML_CCZ4::gt13_bound = "skip"
ML_CCZ4::gt13_bound_speed = 1.0
ML_CCZ4::gt13_bound_limit = 0.0
ML_CCZ4::gt13_bound_scalar = 0.0
ML_CCZ4::gt22_bound = "skip"
ML_CCZ4::gt22_bound_speed = 1.0
ML_CCZ4::gt22_bound_limit = 0.0
ML_CCZ4::gt22_bound_scalar = 0.0
ML_CCZ4::gt23_bound = "skip"
ML_CCZ4::gt23_bound_speed = 1.0
ML_CCZ4::gt23_bound_limit = 0.0
ML_CCZ4::gt23_bound_scalar = 0.0
ML_CCZ4::gt33_bound = "skip"
ML_CCZ4::gt33_bound_speed = 1.0
ML_CCZ4::gt33_bound_limit = 0.0
ML_CCZ4::gt33_bound_scalar = 0.0
ML_CCZ4::Xt1_bound = "skip"
ML_CCZ4::Xt1_bound_speed = 1.0
ML_CCZ4::Xt1_bound_limit = 0.0
ML_CCZ4::Xt1_bound_scalar = 0.0
ML_CCZ4::Xt2_bound = "skip"
ML_CCZ4::Xt2_bound_speed = 1.0
ML_CCZ4::Xt2_bound_limit = 0.0
ML_CCZ4::Xt2_bound_scalar = 0.0
ML_CCZ4::Xt3_bound = "skip"
ML_CCZ4::Xt3_bound_speed = 1.0
ML_CCZ4::Xt3_bound_limit = 0.0
ML_CCZ4::Xt3_bound_scalar = 0.0
ML_CCZ4::trK_bound = "skip"
ML_CCZ4::trK_bound_speed = 1.0
ML_CCZ4::trK_bound_limit = 0.0
ML_CCZ4::trK_bound_scalar = 0.0
ML_CCZ4::At11_bound = "skip"
ML_CCZ4::At11_bound_speed = 1.0
ML_CCZ4::At11_bound_limit = 0.0
ML_CCZ4::At11_bound_scalar = 0.0
ML_CCZ4::At12_bound = "skip"
ML_CCZ4::At12_bound_speed = 1.0
ML_CCZ4::At12_bound_limit = 0.0
ML_CCZ4::At12_bound_scalar = 0.0
ML_CCZ4::At13_bound = "skip"
ML_CCZ4::At13_bound_speed = 1.0
ML_CCZ4::At13_bound_limit = 0.0
ML_CCZ4::At13_bound_scalar = 0.0
ML_CCZ4::At22_bound = "skip"
ML_CCZ4::At22_bound_speed = 1.0
ML_CCZ4::At22_bound_limit = 0.0
ML_CCZ4::At22_bound_scalar = 0.0
ML_CCZ4::At23_bound = "skip"
ML_CCZ4::At23_bound_speed = 1.0
ML_CCZ4::At23_bound_limit = 0.0
ML_CCZ4::At23_bound_scalar = 0.0
ML_CCZ4::At33_bound = "skip"
ML_CCZ4::At33_bound_speed = 1.0
ML_CCZ4::At33_bound_limit = 0.0
ML_CCZ4::At33_bound_scalar = 0.0
ML_CCZ4::Theta_bound = "skip"
ML_CCZ4::Theta_bound_speed = 1.0
ML_CCZ4::Theta_bound_limit = 0.0
ML_CCZ4::Theta_bound_scalar = 0.0
ML_CCZ4::alpha_bound = "skip"
ML_CCZ4::alpha_bound_speed = 1.0
ML_CCZ4::alpha_bound_limit = 0.0
ML_CCZ4::alpha_bound_scalar = 0.0
ML_CCZ4::A_bound = "skip"
ML_CCZ4::A_bound_speed = 1.0
ML_CCZ4::A_bound_limit = 0.0
ML_CCZ4::A_bound_scalar = 0.0
ML_CCZ4::beta1_bound = "skip"
ML_CCZ4::beta1_bound_speed = 1.0
ML_CCZ4::beta1_bound_limit = 0.0
ML_CCZ4::beta1_bound_scalar = 0.0
ML_CCZ4::beta2_bound = "skip"
ML_CCZ4::beta2_bound_speed = 1.0
ML_CCZ4::beta2_bound_limit = 0.0
ML_CCZ4::beta2_bound_scalar = 0.0
ML_CCZ4::beta3_bound = "skip"
ML_CCZ4::beta3_bound_speed = 1.0
ML_CCZ4::beta3_bound_limit = 0.0
ML_CCZ4::beta3_bound_scalar = 0.0
ML_CCZ4::B1_bound = "skip"
ML_CCZ4::B1_bound_speed = 1.0
ML_CCZ4::B1_bound_limit = 0.0
ML_CCZ4::B1_bound_scalar = 0.0
ML_CCZ4::B2_bound = "skip"
ML_CCZ4::B2_bound_speed = 1.0
ML_CCZ4::B2_bound_limit = 0.0
ML_CCZ4::B2_bound_scalar = 0.0
ML_CCZ4::B3_bound = "skip"
ML_CCZ4::B3_bound_speed = 1.0
ML_CCZ4::B3_bound_limit = 0.0
ML_CCZ4::B3_bound_scalar = 0.0
# =============================================================================
# Analysis
# =============================================================================
ADMMass::ADMMass_number = 1
ADMMAss::ADMMass_use_surface_distance_as_volume_radius = "no"
ADMMass::ADMMass_use_all_volume_as_volume_radius = "yes"
ADMMass::ADMMass_Excise_Horizons = "no"
ADMMass::ADMMass_debug = "no"
# =============================================================================
# Checkpoint
# =============================================================================
CarpetIOHDF5::checkpoint = "yes"
CarpetIOHDF5::use_reflevels_from_checkpoint = "yes"
IOUtil::abort_on_io_errors = "yes"
IOUtil::checkpoint_on_terminate = "yes"
IOUtil::checkpoint_every_walltime_hours = 6.
IOUtil::checkpoint_keep = 1
IOUtil::recover = "autoprobe"
IOUtil::checkpoint_dir = "../checkpoints"
IOUtil::recover_dir = "../checkpoints"
# =============================================================================
# Output
# =============================================================================
IOUtil::out_dir = "data"
CarpetIOHDF5::out1D_dir = "data/out_1D"
CarpetIOHDF5::out2D_dir = "data/out_2D"
CarpetIOHDF5::out3D_dir = "data/out_3D"
IOUtil::strict_io_parameter_check = "yes"
IOUtil::parfile_write = "copy"
CarpetIOScalar::one_file_per_group = "yes"
CarpetIOScalar::all_reductions_in_one_file = "yes"
CarpetIOASCII::one_file_per_group = "yes"
CarpetIOHDF5::one_file_per_group = "yes"
IO::out_single_precision = "yes"
CarpetIOHDF5::output_symmetry_points = "yes"
CarpetIOHDF5::output_ghost_points = "yes"
CarpetIOHDF5::output_boundary_points = "yes"
CarpetIOHDF5::output_buffer_points = "no"
CarpetIOHDF5::compression_level = 0
CarpetIOHDF5::minimum_size_for_compression = 1000000
CarpetIOBasic::outinfo_vars = "
Carpet::physical_time_per_hour
ADMBase::lapse
HydroBase::rho
"
CarpetIOScalar::outscalar_reductions = "
minimum
maximum
norm_inf
norm1
norm2
"
CarpetIOScalar::outscalar_vars = "
ADMBase::lapse{ reductions = 'minimum' }
ADMBase::shift{ reductions = 'norm_inf' }
HydroBase::rho{ reductions = 'minimum maximum' }
HydroBase::w_lorentz{ reductions = 'minimum maximum' }
ML_CCZ4::ML_Ham{ reductions = 'norm_inf' }
ML_CCZ4::ML_Mom{ reductions = 'norm_inf' }
ML_CCZ4::ML_log_confac{ reductions = 'maximum minimum' }
"
CarpetIOASCII::out0D_vars = "
Carpet::timing
ADMMass::ADMMass_Masses
"
CarpetIOHDF5::out1d_d = "no"
CarpetIOHDF5::out1d_x = "yes"
CarpetIOHDF5::out1d_y = "yes"
CarpetIOHDF5::out1d_z = "yes"
CarpetIOHDF5::out1d_vars = "
ADMBase::lapse
ADMBase::shift
ADMBase::metric
ADMBase::curv
HydroBase::rho
HydroBase::vel
HydroBase::w_lorentz
HydroBase::press
HydroBase::eps
THC_Core::dens
THC_Core::scon
THC_Core::tau
TmunuBase::stress_energy_scalar
TmunuBase::stress_energy_vector
TmunuBase::stress_energy_tensor
ML_CCZ4::ML_Ham
ML_CCZ4::ML_Mom
ML_CCZ4::ML_log_confac
ADMMass::ADMMass_VolumeMass_GF
CarpetReduce::weight
"
CarpetIOHDF5::out2d_xy = "no"
CarpetIOHDF5::out2d_xz = "yes"
CarpetIOHDF5::out2d_yz = "no"
CarpetIOHDF5::out2d_vars = "
ADMBase::lapse
ADMBase::shift
ADMBase::metric
ADMBase::curv
HydroBase::rho
HydroBase::vel
HydroBase::w_lorentz
HydroBase::press
HydroBase::eps
THC_Core::dens
THC_Core::scon
THC_Core::tau
TmunuBase::stress_energy_scalar
TmunuBase::stress_energy_vector
TmunuBase::stress_energy_tensor
ML_CCZ4::ML_Ham
ML_CCZ4::ML_Mom
ML_CCZ4::ML_log_confac
ADMMass::ADMMass_VolumeMass_GF
CarpetReduce::weight
"
CarpetIOBasic::outinfo_every = 1024
CarpetIOScalar::outscalar_every = 256
CarpetIOASCII::out0d_every = 256
CarpetIOHDF5::out1d_every = 2048
CarpetIOHDF5::out2d_every = 1024
CarpetIOHDF5::out_every = -1
"""
open(re.sub(r'(.*)\.rpar$', r'\1.par', sys.argv[0]), 'w').write(re.sub(r'\n *',r'\n',Template(Template(lines).substitute(locals())).substitute(locals())))
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Ham_nans.pdf
Type: application/pdf
Size: 193804 bytes
Desc: not available
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20200914/85c70a36/attachment-0001.pdf
More information about the Users
mailing list