diff options
Diffstat (limited to 'par/CCZ4/rbst_ccz4.par')
| -rw-r--r-- | par/CCZ4/rbst_ccz4.par | 230 |
1 files changed, 230 insertions, 0 deletions
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 + |