From 18fe14998dac0391740f14a2bc5a8a2f014eb36b Mon Sep 17 00:00:00 2001 From: Frank Löffler Date: Tue, 29 Oct 2013 03:52:57 -0500 Subject: par/CCZ4: remove parfiles because they don't work and never worked. See trac #1387. --- par/CCZ4/bbh_ccz4.par | 493 ------------------------------------------------- par/CCZ4/gw1d_ccz4.par | 283 ---------------------------- par/CCZ4/rbst_ccz4.par | 230 ----------------------- 3 files changed, 1006 deletions(-) delete mode 100644 par/CCZ4/bbh_ccz4.par delete mode 100644 par/CCZ4/gw1d_ccz4.par delete mode 100644 par/CCZ4/rbst_ccz4.par diff --git a/par/CCZ4/bbh_ccz4.par b/par/CCZ4/bbh_ccz4.par deleted file mode 100644 index e009fd5..0000000 --- a/par/CCZ4/bbh_ccz4.par +++ /dev/null @@ -1,493 +0,0 @@ -################################################################################ - Cactus::cctk_run_title = "CCZ4 Binary Black Hole" -################################################################################ - -ActiveThorns = " - ADMBase - ADMCoupling - ADMMacros - AEILocalInterp - AHFinderDirect - BLAS - Boundary - Carpet - CarpetInterp - CarpetIOASCII - CarpetIOBasic - CarpetIOHDF5 - CarpetIOScalar - CarpetLib - CarpetMask - CarpetReduce - CarpetRegrid2 - CarpetSlab - CarpetTracker - CartGrid3D - CoordBase - CoordGauge - Dissipation - Formaline - Fortran - GenericFD - GSL - HDF5 - InitBase - IOUtil - LAPACK - LocalInterp - LoopControl - ML_ADMConstraints - ML_BSSN - ML_BSSN_Helper - MoL - Multipole - NaNChecker - NewRad - PunctureTracker - QuasiLocalMeasures - ReflectionSymmetry - RotatingSymmetry180 - Slab - SpaceMask - SphericalSurface - StaticConformal - SummationByParts - SymBase - SystemStatistics - TerminationTrigger - Time - TimerReport - TmunuBase - TwoPunctures - Vectors - WeylScal4 -" - -################################################################################ -# Grid structure -################################################################################ - -CartGrid3D::type = "coordbase" -Carpet::domain_from_coordbase = "yes" -CoordBase::domainsize = "minmax" - -CoordBase::xmin = 0.00 -CoordBase::ymin = -120.00 -CoordBase::zmin = 0.00 -CoordBase::xmax = +120.00 -CoordBase::ymax = +120.00 -CoordBase::zmax = +120.00 -CoordBase::dx = 2.00 -CoordBase::dy = 2.00 -CoordBase::dz = 2.00 - -Driver::ghost_size = 3 - -CoordBase::boundary_size_x_lower = 3 -CoordBase::boundary_size_y_lower = 3 -CoordBase::boundary_size_z_lower = 3 -CoordBase::boundary_size_x_upper = 3 -CoordBase::boundary_size_y_upper = 3 -CoordBase::boundary_size_z_upper = 3 - -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -CarpetRegrid2::symmetry_rotating180 = yes - -ReflectionSymmetry::reflection_z = yes -ReflectionSymmetry::avoid_origin_z = no - -Time::dtfac = 0.25 - -################################################################################ -# Mesh refinement -################################################################################ - -Carpet::max_refinement_levels = 7 -CarpetRegrid2::num_centres = 2 -CarpetRegrid2::num_levels_1 = 7 -CarpetRegrid2::position_x_1 = 3 -CarpetRegrid2::radius_1 [1] = 24.0 -CarpetRegrid2::radius_1 [2] = 12.0 -CarpetRegrid2::radius_1 [3] = 6.0 -CarpetRegrid2::radius_1 [4] = 3.0 -CarpetRegrid2::radius_1 [5] = 1.5 -CarpetRegrid2::radius_1 [6] = 0.6 -CarpetRegrid2::num_levels_2 = 7 -CarpetRegrid2::position_x_2 = -3 -CarpetRegrid2::radius_2 [1] = 24.0 -CarpetRegrid2::radius_2 [2] = 12.0 -CarpetRegrid2::radius_2 [3] = 6.0 -CarpetRegrid2::radius_2 [4] = 3.0 -CarpetRegrid2::radius_2 [5] = 1.5 -CarpetRegrid2::radius_2 [6] = 0.6 - -Carpet::use_buffer_zones = yes -Carpet::prolongation_order_space = 5 -Carpet::prolongation_order_time = 2 - -CarpetRegrid2::regrid_every = 64 -CarpetRegrid2::verbose = yes -Carpet::grid_coordinates_filename = "carpet-grid.asc" - -Carpet::time_refinement_factors = "[1,1,2,4,8,16,32,64,128,256]" - -############################################################################### -# Initial Data -############################################################################### - -ADMBase::initial_data = "twopunctures" -ADMBase::initial_lapse = "twopunctures-averaged" -ADMBase::initial_shift = "zero" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" - -TwoPunctures::target_M_plus = 0.5 -TwoPunctures::target_M_minus = 0.5 - -TwoPunctures::par_m_plus = 0.476534633024028 -TwoPunctures::par_m_minus = 0.476534633024028 - -TwoPunctures::par_b = 3 -TwoPunctures::center_offset[0] = 0 - -TwoPunctures::par_P_plus[0] = -0.0058677669328272 -TwoPunctures::par_P_plus[1] = 0.138357448824906 -TwoPunctures::par_P_plus[2] = 0. - -TwoPunctures::par_P_minus[0] = 0.0058677669328272 -TwoPunctures::par_P_minus[1] = -0.138357448824906 -TwoPunctures::par_P_minus[2] = 0. - -TwoPunctures::par_S_plus[0] = 0. -TwoPunctures::par_S_plus[1] = 0. -TwoPunctures::par_S_plus[2] = 0. - -TwoPunctures::par_S_minus[0] = 0. -TwoPunctures::par_S_minus[1] = 0. -TwoPunctures::par_S_minus[2] = 0. - -TwoPunctures::give_bare_mass = yes -TwoPunctures::TP_epsilon = 1.0e-6 - -Carpet::init_fill_timelevels = yes -InitBase::initial_data_setup_method = "init_single_level" - -############################################################################### -# Evolution -############################################################################### - -ADMBase::evolution_method = "ML_BSSN" -ADMBase::lapse_evolution_method = "ML_BSSN" -ADMBase::shift_evolution_method = "ML_BSSN" -ADMBase::dtlapse_evolution_method = "ML_BSSN" -ADMBase::dtshift_evolution_method = "ML_BSSN" - -ML_BSSN::GammaShift = 1.0 -ML_BSSN::dampk1 = 1.0 -ML_BSSN::harmonicN = 1 # 1+log -ML_BSSN::harmonicF = 2.0 # 1+log -ML_BSSN::ShiftGammaCoeff = 0.75 -ML_BSSN::BetaDriver = 1.0 -ML_BSSN::LapseAdvectionCoeff = 1.0 -ML_BSSN::ShiftAdvectionCoeff = 1.0 - -ML_BSSN::MinimumLapse = 1.0e-8 -ML_BSSN::conformalMethod = 1 # 1 for W -ML_BSSN::formulation = 1 # 1 for CCZ4 -ML_BSSN::apply_dissipation = never -ML_BSSN::dt_lapse_shift_method = "noLapseShiftAdvection" - -################################################################################ -# Boundary conditions -################################################################################ - -ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" -ML_BSSN::my_rhs_boundary_condition = "NewRad" -Boundary::radpower = 2 - -ML_BSSN::ML_log_confac_bound = "none" -ML_BSSN::ML_metric_bound = "none" -ML_BSSN::ML_Gamma_bound = "none" -ML_BSSN::ML_trace_curv_bound = "none" -ML_BSSN::ML_curv_bound = "none" -ML_BSSN::ML_lapse_bound = "none" -ML_BSSN::ML_dtlapse_bound = "none" -ML_BSSN::ML_shift_bound = "none" -ML_BSSN::ML_dtshift_bound = "none" -ML_BSSN::ML_Theta_bound = "none" - -################################################################################ -# BH tracking -################################################################################ - -CarpetTracker::surface [0] = 0 -CarpetTracker::surface [1] = 1 -PunctureTracker::track [0] = yes -PunctureTracker::initial_x [0] = 3 -PunctureTracker::which_surface_to_store_info[0] = 0 -PunctureTracker::track [1] = yes -PunctureTracker::initial_x [1] = -3 -PunctureTracker::which_surface_to_store_info[1] = 1 - -############################################################################### -# Spatial finite differencing -############################################################################### - -SummationByParts::order = 4 -ML_BSSN::fdOrder = 4 - -Dissipation::order = 5 -Dissipation::vars = " - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ML_BSSN::ML_lapse - ML_BSSN::ML_shift - ML_BSSN::ML_dtlapse - ML_BSSN::ML_dtshift - ML_BSSN::ML_Theta -" - -############################################################################### -# Time Integration -############################################################################### - -MoL::ODE_Method = "RK4" -MoL::MoL_Intermediate_Steps = 4 -MoL::MoL_Num_Scratch_Levels = 1 - -################################################################################ -# Psi4 mode decomposition by Multipole -################################################################################ - -WeylScal4::fd_order = "4th" - -Multipole::nradii = 4 -Multipole::radius[0] = 30 -Multipole::radius[1] = 50 -Multipole::radius[2] = 80 -Multipole::radius[3] = 100 -Multipole::ntheta = 120 -Multipole::nphi = 240 -Multipole::variables = " - WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i' name='psi4'} -" - -Multipole::out_every = 32 -Multipole::l_max = 8 -Multipole::output_hdf5 = yes -Multipole::output_ascii = no - -################################################################################ -# Apparent Horizons -################################################################################ - - -AHFinderDirect::N_horizons = 3 -AHFinderDirect::find_every = 128 -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::verbose_level = "physics details" -AHFinderDirect::move_origins = yes - -AHFinderDirect::origin_x [1] = 3 -AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 3 -AHFinderDirect::initial_guess__coord_sphere__radius [1] = 0.25 -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] = "PunctureTracker::pt_loc_x[0]" -AHFinderDirect::track_origin_source_y [1] = "PunctureTracker::pt_loc_y[0]" -AHFinderDirect::track_origin_source_z [1] = "PunctureTracker::pt_loc_z[0]" -AHFinderDirect::max_allowable_horizon_radius [1] = 3 - -AHFinderDirect::origin_x [2] = -3 -AHFinderDirect::initial_guess__coord_sphere__x_center[2] = -3 -AHFinderDirect::initial_guess__coord_sphere__radius [2] = 0.25 -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] = "PunctureTracker::pt_loc_x[1]" -AHFinderDirect::track_origin_source_y [2] = "PunctureTracker::pt_loc_y[1]" -AHFinderDirect::track_origin_source_z [2] = "PunctureTracker::pt_loc_z[1]" -AHFinderDirect::max_allowable_horizon_radius [2] = 3 - -AHFinderDirect::origin_x [3] = 0 -AHFinderDirect::find_after_individual [3] = 0 -AHFinderDirect::initial_guess__coord_sphere__x_center[3] = 0 -AHFinderDirect::initial_guess__coord_sphere__radius [3] = 1.0 -AHFinderDirect::which_surface_to_store_info [3] = 4 -AHFinderDirect::set_mask_for_individual_horizon [3] = no -AHFinderDirect::max_allowable_horizon_radius [3] = 6 - -################################################################################ -# Spherical surfaces -################################################################################ - -SphericalSurface::nsurfaces = 5 -SphericalSurface::maxntheta = 66 -SphericalSurface::maxnphi = 124 -SphericalSurface::verbose = no - -# Surfaces 0 and 1 are used by PunctureTracker - -# Horizon 1 -SphericalSurface::ntheta [2] = 41 -SphericalSurface::nphi [2] = 80 -SphericalSurface::nghoststheta [2] = 2 -SphericalSurface::nghostsphi [2] = 2 - -# Horizon 2 -SphericalSurface::ntheta [3] = 41 -SphericalSurface::nphi [3] = 80 -SphericalSurface::nghoststheta [3] = 2 -SphericalSurface::nghostsphi [3] = 2 - -# Common horizon -SphericalSurface::ntheta [4] = 41 -SphericalSurface::nphi [4] = 80 -SphericalSurface::nghoststheta [4] = 2 -SphericalSurface::nghostsphi [4] = 2 - -################################################################################ -# Isolated Horizons -################################################################################ - -QuasiLocalMeasures::verbose = yes -QuasiLocalMeasures::veryverbose = no -QuasiLocalMeasures::interpolator = "Lagrange polynomial interpolation" -QuasiLocalMeasures::interpolator_options = "order=4" -QuasiLocalMeasures::spatial_order = 4 -QuasiLocalMeasures::num_surfaces = 3 -QuasiLocalMeasures::surface_index [0] = 2 -QuasiLocalMeasures::surface_index [1] = 3 -QuasiLocalMeasures::surface_index [2] = 4 - -################################################################################ -# Correctness Checking -################################################################################ - -Carpet::poison_new_timelevels = yes -Carpet::poison_value = 100 -CarpetLib::poison_new_memory = yes -CarpetLib::poison_value = 200 -Carpet::check_for_poison = no - -NaNChecker::check_every = 128 -NanChecker::check_after = 0 -NaNChecker::report_max = 10 -NaNChecker::verbose = "all" -NaNChecker::action_if_found = "terminate" -NaNChecker::out_NaNmask = yes -NaNChecker::check_vars = " - ML_BSSN::ML_log_confac -" - -################################################################################ -# Timers -################################################################################ - -Cactus::cctk_timer_output = full -TimerReport::out_every = 256 -TimerReport::n_top_timers = 40 -TimerReport::output_all_timers_together = yes -TimerReport::output_schedule_timers = no - -################################################################################ -# Output -################################################################################ - -IO::out_dir = $parfile -IOScalar::one_file_per_group = yes -IOASCII::one_file_per_group = yes - -IOBasic::outInfo_every = 1 -IOBasic::outInfo_reductions = "minimum maximum" -IOBasic::outInfo_vars = " - Carpet::physical_time_per_hour - ML_BSSN::ML_log_confac - SystemStatistics::maxrss_mb -" - -IOScalar::outScalar_every = 128 -IOScalar::outScalar_reductions = "minimum maximum average" -IOScalar::outScalar_vars = "SystemStatistics::process_memory_mb" - -IOASCII::out0D_every = 32 -IOASCII::out0D_vars = " - Carpet::timing - PunctureTracker::pt_loc - QuasiLocalMeasures::qlm_scalars{out_every = 128} -" - -IOHDF5::out1D_every = 256 -IOHDF5::out1D_d = no -IOHDF5::out1D_x = yes -IOHDF5::out1D_y = no -IOHDF5::out1D_z = no -IOHDF5::out1D_vars = " - ML_BSSN::ML_log_confac -" - -IOHDF5::out2D_every = 1024 -IOHDF5::out2D_xy = yes -IOHDF5::out2D_xz = no -IOHDF5::out2D_yz = no -IOHDF5::out2D_vars = " - ML_BSSN::ML_log_confac -" - -IOHDF5::out_every = 8192 -IOHDF5::out_vars = " - ML_BSSN::ML_log_confac -" - -Carpet::verbose = no -Carpet::veryverbose = no -Carpet::schedule_barriers = no -Carpet::storage_verbose = no -Carpet::timers_verbose = no -CarpetLib::output_bboxes = no - -Cactus::cctk_full_warnings = yes -Cactus::highlight_warning_messages = no - -############################################################################### -# Checkpointing and recovery -############################################################################### - -CarpetIOHDF5::checkpoint = yes -IO::checkpoint_ID = no -IO::recover = "autoprobe" -IO::checkpoint_every_walltime_hours = 3 -IO::out_proc_every = 2 -IO::checkpoint_keep = 1 -IO::checkpoint_on_terminate = yes -IO::checkpoint_dir = "../checkpoints" -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 = @WALLTIME_HOURS@ # hours -TerminationTrigger::on_remaining_walltime = 30 # minutes -TerminationTrigger::output_remtime_every_minutes = 30 -Cactus::terminate = "time" -Cactus::cctk_final_time = 300.0 diff --git a/par/CCZ4/gw1d_ccz4.par b/par/CCZ4/gw1d_ccz4.par deleted file mode 100644 index e74f951..0000000 --- a/par/CCZ4/gw1d_ccz4.par +++ /dev/null @@ -1,283 +0,0 @@ -Cactus::cctk_run_title = "CCZ4 1D Gauge Wave" - -ActiveThorns = " - ADMBase - ADMCoupling - ADMMacros - Boundary - Carpet - CarpetInterp - CarpetInterp2 - CarpetIOASCII - CarpetIOBasic - CarpetIOHDF5 - CarpetIOScalar - CarpetLib - CarpetReduce - CarpetRegrid - CarpetSlab - CartGrid3d - CoordBase - CoordGauge - Dissipation - Exact - GenericFD - IOUtil - LoopControl - ML_BSSN - ML_BSSN_Helper - MoL - NanChecker - Periodic - Slab - SpaceMask - SphericalSurface - StaticConformal - SummationByParts - SymBase - Time - TimerReport - TMuNuBase -" - - -############################################################# -# Evolution -############################################################# - -ADMBase::evolution_method = "ML_BSSN" -ADMBase::lapse_evolution_method = "ML_BSSN" -ADMBase::shift_evolution_method = "ML_BSSN" -ADMBase::dtlapse_evolution_method = "ML_BSSN" -ADMBase::dtshift_evolution_method = "ML_BSSN" - -ML_BSSN::GammaShift = 0.0 -ML_BSSN::dampk1 = 1.0 -ML_BSSN::harmonicN = 2 -ML_BSSN::harmonicF = 1.0 -ML_BSSN::ShiftGammaCoeff = 0.0 # beta to zero -ML_BSSN::ShiftBCoeff = 0.0 # B to zero -ML_BSSN::BetaDriver = 0.0 # damp in beta -ML_BSSN::LapseACoeff = 0.0 # A to zero -ML_BSSN::LapseAdvectionCoeff = 0.0 # no advection terms -ML_BSSN::ShiftAdvectionCoeff = 0.0 # no advection terms - -ML_BSSN::MinimumLapse = 1.0e-8 -ML_BSSN::conformalMethod = 1 # 1 for W, 0 for phi -ML_BSSN::formulation = 1 # 1 for CCZ4 - -############################################################# -# Boundary conditions -############################################################# - -Boundary::radpower = 2 - -ML_BSSN::ML_log_confac_bound = "none" -ML_BSSN::ML_metric_bound = "none" -ML_BSSN::ML_Gamma_bound = "none" -ML_BSSN::ML_trace_curv_bound = "none" -ML_BSSN::ML_curv_bound = "none" -ML_BSSN::ML_lapse_bound = "none" -ML_BSSN::ML_dtlapse_bound = "none" -ML_BSSN::ML_shift_bound = "none" -ML_BSSN::ML_dtshift_bound = "none" -ML_BSSN::ML_Theta_bound = "none" - -############################################################# -# Spatial derivatives -############################################################# - -SummationByParts::order = 4 -SummationByParts::sbp_upwind_deriv = no # Drop order instead, only for - # advection derivatives -SummationByParts::sbp_1st_deriv = no -SummationByParts::sbp_2nd_deriv = no -SummationByParts::onesided_interpatch_boundaries = no -SummationByParts::onesided_outer_boundaries = no -SummationByParts::use_dissipation = yes -SummationByParts::scale_with_h = no -SummationByParts::dissipation_type = "Kreiss-Oliger" -SummationByParts::epsdis = 0.1 -SummationByParts::vars = " - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ML_BSSN::ML_lapse - ML_BSSN::ML_shift - ML_BSSN::ML_dtlapse - ML_BSSN::ML_dtshift - ML_BSSN::ML_Theta -" - -############################################################# -# Grid -############################################################# - -CoordBase::domainsize = minmax - -CoordBase::boundary_size_x_lower = 3 -CoordBase::boundary_size_y_lower = 3 -CoordBase::boundary_size_z_lower = 3 -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -CoordBase::boundary_size_x_upper = 3 -CoordBase::boundary_size_y_upper = 3 -CoordBase::boundary_size_z_upper = 3 -CoordBase::boundary_shiftout_x_upper = 0 -CoordBase::boundary_shiftout_y_upper = 0 -CoordBase::boundary_shiftout_z_upper = 0 - -CartGrid3D::type = "coordbase" -CartGrid3D::domain = "full" -CartGrid3D::avoid_origin = "no" - -Periodic::periodic = "yes" - -CoordBase::xmin = -0.5 -CoordBase::ymin = -0.5 -CoordBase::zmin = -0.5 - -CoordBase::xmax = 0.5 -CoordBase::ymax = 0.5 -CoordBase::zmax = 0.5 - -CoordBase::dx = 0.005 -CoordBase::dy = 0.1 -CoordBase::dz = 0.1 - -############################################################# -# CarpetRegrid -############################################################# - -CarpetRegrid::refinement_levels = 1 -CarpetRegrid::refined_regions = "manual-coordinate-list" -CarpetRegrid::smart_outer_boundaries = yes -CarpetRegrid::keep_same_grid_structure = "yes" - -############################################################# -# Carpet -############################################################# - -Carpet::ghost_size = 3 -Carpet::domain_from_coordbase = "yes" -Carpet::max_refinement_levels = 1 -Carpet::init_3_timelevels = "no" -Carpet::prolongation_order_space = 3 -Carpet::prolongation_order_time = 2 -Carpet::use_buffer_zones = yes -Carpet::num_integrator_substeps = 4 - -############################################################# -# CarpetLib -############################################################# - -############################################################# -# Time integration -############################################################# - -Cactus::terminate = "time" -Cactus::cctk_final_time = 1000 -# Use 0.5 for RK4 and 0.25 for ICN -Time::dtfac = 0.5 -MethodOfLines::ode_method = "RK4" -MethodOfLines::MoL_Intermediate_Steps = 4 -MethodOfLines::MoL_Num_Scratch_Levels = 1 - -MethodOfLines::MoL_NaN_Check = "yes" -MethodOfLines::initial_data_is_crap = "no" - -############################################################# -# Initial data -############################################################# - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" - -Exact::exact_model = "Minkowski/gauge wave" -Exact::Minkowski_gauge_wave__what_fn = "sin" -Exact::Minkowski_gauge_wave__amplitude = 0.1 -Exact::Minkowski_gauge_wave__omega = 1 -Exact::Minkowski_gauge_wave__lambda = 1 -Exact::Minkowski_gauge_wave__diagonal = "no" - -################################################################################ -# Constraints -################################################################################ - -################################################################################ -# Timers -################################################################################ - -Cactus::cctk_timer_output = full -TimerReport::out_every = 4000 -TimerReport::n_top_timers = 40 -TimerReport::out_filename = "TimerReport" -TimerReport::output_all_timers = yes - -############################################################# -# Output -############################################################# - -IO::out_dir = $parfile -IO::out_fileinfo = "none" - -CarpetIOScalar::outScalar_every = 100 -CarpetIOScalar::outScalar_reductions = "norm_inf norm2" -CarpetIOScalar::outScalar_vars = " - ADMBase::alp - ADMBase::shift - ADMBase::metric - ADMBase::curv - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ML_BSSN::ML_lapse - ML_BSSN::ML_dtlapse - ML_BSSN::ML_dtshift - ML_BSSN::ML_shift - ML_BSSN::ML_Theta - ML_BSSN::ML_Ham - ML_BSSN::ML_mom - ML_BSSN::ML_cons_detg - ML_BSSN::ML_cons_traceA - ML_BSSN::ML_cons_Gamma -" - -CarpetIOBasic::outInfo_every = 1 -CarpetIOBasic::outInfo_vars = " - ADMBase::alp -" - -IOASCII::out1D_every = 4000 -IOASCII::out1D_x = "yes" -IOASCII::out1D_y = "yes" -IOASCII::out1D_z = "yes" -IOASCII::out1D_d = "no" -IOASCII::out1D_vars = " - ADMBase::alp - ADMBase::metric - ADMBase::curv - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ML_BSSN::ML_lapse - ML_BSSN::ML_dtlapse - ML_BSSN::ML_Theta - ML_BSSN::ML_Ham - ML_BSSN::ML_mom - ML_BSSN::ML_cons_detg - ML_BSSN::ML_cons_traceA - ML_BSSN::ML_cons_Gamma -" - -CarpetIOASCII::out_precision = 19 - diff --git a/par/CCZ4/rbst_ccz4.par b/par/CCZ4/rbst_ccz4.par deleted file mode 100644 index 3801ba4..0000000 --- a/par/CCZ4/rbst_ccz4.par +++ /dev/null @@ -1,230 +0,0 @@ -Cactus::cctk_run_title = "CCZ4 Robust Stability Test" - -ActiveThorns = " - ADMBase - ADMCoupling - ADMMacros - Boundary - Carpet - CarpetInterp - CarpetIOASCII - CarpetIOBasic - CarpetIOHDF5 - CarpetIOScalar - CarpetLib - CarpetReduce - CarpetRegrid - CarpetSlab - CartGrid3d - CoordBase - CoordGauge - Exact - GenericFD - IOUtil - LoopControl - ML_BSSN - ML_BSSN_Helper - MoL - NanChecker - Noise - Periodic - Slab - SpaceMask - SphericalSurface - StaticConformal - SymBase - Time - TimerReport - TMuNuBase -" - -ADMBase::evolution_method = "ML_BSSN" -ADMBase::lapse_evolution_method = "ML_BSSN" -ADMBase::shift_evolution_method = "ML_BSSN" -ADMBase::dtlapse_evolution_method = "ML_BSSN" -ADMBase::dtshift_evolution_method = "ML_BSSN" - -ML_BSSN::GammaShift = 0.0 -ML_BSSN::dampk1 = 0.0 -ML_BSSN::harmonicN = 2 -ML_BSSN::harmonicF = 1.0 -ML_BSSN::ShiftGammaCoeff = 0.0 -ML_BSSN::ShiftBCoeff = 0.0 # B to zero -ML_BSSN::BetaDriver = 0.0 -ML_BSSN::LapseACoeff = 0.0 # A to zero -ML_BSSN::LapseAdvectionCoeff = 0.0 -ML_BSSN::ShiftAdvectionCoeff = 0.0 - -ML_BSSN::conformalMethod = 1 # 1 for W -ML_BSSN::formulation = 1 # 1 for CCZ4 - -Boundary::radpower = 2 - -ML_BSSN::ML_log_confac_bound = "none" -ML_BSSN::ML_metric_bound = "none" -ML_BSSN::ML_Gamma_bound = "none" -ML_BSSN::ML_trace_curv_bound = "none" -ML_BSSN::ML_curv_bound = "none" -ML_BSSN::ML_lapse_bound = "none" -ML_BSSN::ML_dtlapse_bound = "none" -ML_BSSN::ML_shift_bound = "none" -ML_BSSN::ML_dtshift_bound = "none" -ML_BSSN::ML_Theta_bound = "none" - -############################################################# -# Grid -############################################################# - -CoordBase::domainsize = minmax - -CoordBase::boundary_size_x_lower = 5 -CoordBase::boundary_size_y_lower = 5 -CoordBase::boundary_size_z_lower = 5 -CoordBase::boundary_shiftout_x_lower = 1 -CoordBase::boundary_shiftout_y_lower = 1 -CoordBase::boundary_shiftout_z_lower = 1 - -CoordBase::boundary_size_x_upper = 5 -CoordBase::boundary_size_y_upper = 5 -CoordBase::boundary_size_z_upper = 5 -CoordBase::boundary_shiftout_x_upper = 0 -CoordBase::boundary_shiftout_y_upper = 0 -CoordBase::boundary_shiftout_z_upper = 0 - -CartGrid3D::type = "coordbase" -CartGrid3D::domain = "full" -CartGrid3D::avoid_origin = "no" - -Periodic::periodic = "yes" - -CoordBase::xmin = -0.5 -CoordBase::ymin = -0.5 -CoordBase::zmin = -0.5 - -CoordBase::xmax = 0.5 -CoordBase::ymax = 0.5 -CoordBase::zmax = 0.5 - -CoordBase::dx = 0.02 -CoordBase::dy = 0.02 -CoordBase::dz = 0.02 - -############################################################# -# CarpetRegrid -############################################################# - -CarpetRegrid::refinement_levels = 1 -CarpetRegrid::refined_regions = "manual-coordinate-list" -CarpetRegrid::smart_outer_boundaries = yes -CarpetRegrid::keep_same_grid_structure = "yes" - -############################################################# -# Carpet -############################################################# - -Carpet::ghost_size = 5 -Carpet::domain_from_coordbase = "yes" -Carpet::max_refinement_levels = 1 -Carpet::init_3_timelevels = "no" -Carpet::prolongation_order_space = 3 -Carpet::prolongation_order_time = 2 -Carpet::use_buffer_zones = yes -Carpet::num_integrator_substeps = 4 - -############################################################# -# CarpetLib -############################################################# - -############################################################# -# Time integration -############################################################# - -Cactus::terminate = "time" -Cactus::cctk_final_time = 1000 -# Use 0.5 for RK4 and 0.25 for ICN -Time::dtfac = 0.5 -MethodOfLines::ode_method = "RK4" -MethodOfLines::MoL_Intermediate_Steps = 4 -MethodOfLines::MoL_Num_Scratch_Levels = 1 - -MethodOfLines::MoL_NaN_Check = "yes" -MethodOfLines::initial_data_is_crap = "no" - -############################################################# -# Initial data -############################################################# - -ADMBase::initial_data = "exact" -ADMBase::initial_lapse = "exact" -ADMBase::initial_shift = "exact" - -Exact::exact_model = "Minkowski" - -Noise::apply_id_noise = yes -Noise::id_vars = " - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ML_BSSN::ML_lapse - ML_BSSN::ML_Theta -" -Noise::Amplitude = 1e-10 - -############################################################# -# Output -############################################################# - -IO::out_dir = $parfile -IO::out_fileinfo = "none" - -CarpetIOScalar::outScalar_every = 100 -CarpetIOScalar::outScalar_reductions = "norm_inf norm2" -CarpetIOScalar::outScalar_vars = " - ADMBase::alp - ADMBase::shift - ADMBase::metric - ADMBase::curv - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ML_BSSN::ML_lapse - ML_BSSN::ML_dtlapse - ML_BSSN::ML_dtshift - ML_BSSN::ML_shift - ML_BSSN::ML_Theta - ML_BSSN::ML_Ham - ML_BSSN::ML_mom - ML_BSSN::ML_cons_detg - ML_BSSN::ML_cons_traceA - ML_BSSN::ML_cons_Gamma -" - -CarpetIOBasic::outInfo_every = 100 -CarpetIOBasic::outInfo_vars = "ADMBase::alp" - -IOASCII::out1D_every = 1000 -IOASCII::out1D_x = "yes" -IOASCII::out1D_y = "yes" -IOASCII::out1D_z = "yes" -IOASCII::out1D_d = "no" -IOASCII::out1D_vars = " - ADMBase::metric - ADMBase::curv - ADMBase::alp - ADMBase::shift - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma - ML_BSSN::ML_lapse - ML_BSSN::ML_dtlapse - ML_BSSN::ML_Theta -" - -CarpetIOASCII::out_precision = 19 - -- cgit v1.2.3