From 342f260edadacf1632d3020fa441da1802a96cce Mon Sep 17 00:00:00 2001 From: Barry Wardell Date: Mon, 30 Apr 2012 17:43:36 +0100 Subject: Add some sample CCZ4 parameter files. --- par/CCZ4/bbh_ccz4.par | 493 +++++++++++++++++++++++++++++++++++++++++++++++++ par/CCZ4/gw1d_ccz4.par | 283 ++++++++++++++++++++++++++++ par/CCZ4/rbst_ccz4.par | 230 +++++++++++++++++++++++ 3 files changed, 1006 insertions(+) create mode 100644 par/CCZ4/bbh_ccz4.par create mode 100644 par/CCZ4/gw1d_ccz4.par create mode 100644 par/CCZ4/rbst_ccz4.par diff --git a/par/CCZ4/bbh_ccz4.par b/par/CCZ4/bbh_ccz4.par new file mode 100644 index 0000000..e009fd5 --- /dev/null +++ b/par/CCZ4/bbh_ccz4.par @@ -0,0 +1,493 @@ +################################################################################ + 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 new file mode 100644 index 0000000..e74f951 --- /dev/null +++ b/par/CCZ4/gw1d_ccz4.par @@ -0,0 +1,283 @@ +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 new file mode 100644 index 0000000..3801ba4 --- /dev/null +++ b/par/CCZ4/rbst_ccz4.par @@ -0,0 +1,230 @@ +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