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