diff options
Diffstat (limited to 'ML_BSSN_Test/test/gw3d_McLachlan_ord4_15/gw3d_McLachlan_ord4_15.par')
-rw-r--r-- | ML_BSSN_Test/test/gw3d_McLachlan_ord4_15/gw3d_McLachlan_ord4_15.par | 164 |
1 files changed, 164 insertions, 0 deletions
diff --git a/ML_BSSN_Test/test/gw3d_McLachlan_ord4_15/gw3d_McLachlan_ord4_15.par b/ML_BSSN_Test/test/gw3d_McLachlan_ord4_15/gw3d_McLachlan_ord4_15.par new file mode 100644 index 0000000..b1d9d88 --- /dev/null +++ b/ML_BSSN_Test/test/gw3d_McLachlan_ord4_15/gw3d_McLachlan_ord4_15.par @@ -0,0 +1,164 @@ + +ActiveThorns = "CoordBase SymBase Boundary NanChecker CartGrid3d Time MoL CarpetIOBasic CarpetIOScalar IOUtil Carpet CarpetLib CarpetReduce CarpetInterp CarpetSlab CarpetRegrid CarpetIOASCII ADMBase StaticConformal SpaceMask Slab Periodic Exact CoordGauge ADMCoupling LoopControl ML_BSSN ML_BSSN_Helper Dissipation GenericFD TMuNuBase NewRad SphericalSurface ADMMacros" + + +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::harmonicN = 2 +ML_BSSN::harmonicF = 1.0 +ML_BSSN::ShiftGammaCoeff = 0.0 +ML_BSSN::BetaDriver = 1.0 +ML_BSSN::LapseAdvectionCoeff = 1.0 +ML_BSSN::ShiftAdvectionCoeff = 0.0 + +ML_BSSN::MinimumLapse = 1.0e-8 + +ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" +#ML_BSSN::my_rhs_boundary_condition = "NewRad" +Boundary::radpower = 2 + +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_dtlapse + ML_BSSN::ML_dtshift +" +Dissipation::epsdis = 0.02 + + + +############################################################# +# 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 +CoordBase::ymin = 0 +CoordBase::zmin = 0 + +CoordBase::xmax = 1 +CoordBase::ymax = 1 +CoordBase::zmax = 1 + +CoordBase::dx = 0.0666666666666667 +CoordBase::dy = 0.0666666666666667 +CoordBase::dz = 0.0666666666666667 + +############################################################# +# 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::refinement_factor = 2 +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 = 1 +# Use 0.5 for RK4 and 0.25 for ICN +Time::dtfac = 0.5 + +MethodOfLines::ode_method = "RK4" +MethodOfLines::MoL_Intermediate_Steps = 4 + +#MethodOfLines::ode_method = "generic" +#MethodOfLines::generic_type = "RK" +#MethodOfLines::MoL_Intermediate_Steps = 1 + +MethodOfLines::MoL_Num_Scratch_Levels = 1 +MethodOfLines::MoL_NaN_Check = "no" +MethodOfLines::initial_data_is_crap = "no" + +############################################################# +# Initial data +############################################################# + +ADMBase::initial_data = "exact" +ADMBase::initial_lapse = "exact" +ADMBase::initial_shift = "exact" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +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 = 0.577350269189626 +Exact::Minkowski_gauge_wave__diagonal = "no" + +Exact::rotation_euler_psi = -1.9216757376671543544 +Exact::rotation_euler_theta = 0.66214523564555227398 +Exact::rotation_euler_phi = 1.2199169159226388841 + +Exact::exact_order = 4 +Exact::exact_eps = 2e-4 + +############################################################# +# Output +############################################################# + +IO::out_dir = $parfile +IO::out_fileinfo = "none" + +CarpetIOBasic::outInfo_every = 1 +CarpetIOBasic::outInfo_vars = "ADMBase::alp" + +IOASCII::out1D_every = 30 +IOASCII::out_precision = 19 +IOASCII::out1D_x = "yes" +IOASCII::out1D_y = "yes" +IOASCII::out1D_z = "yes" +IOASCII::out1D_d = "no" +IOASCII::out1D_vars = "ADMBase::gxx ADMBase::kxx ML_BSSN::ML_cons_detg ML_BSSN::ML_cons_Gamma ML_BSSN::ML_cons_traceA ML_BSSN::ML_curv ML_BSSN::ML_dtlapse ML_BSSN::ML_dtshift ML_BSSN::ML_Gamma ML_BSSN::ML_lapse ML_BSSN::ML_log_confac ML_BSSN::ML_metric ML_BSSN::ML_shift ML_BSSN::ML_trace_curv grid::x grid::y grid::z" |