diff options
Diffstat (limited to 'test/schw-0050.par')
| -rw-r--r-- | test/schw-0050.par | 73 |
1 files changed, 50 insertions, 23 deletions
diff --git a/test/schw-0050.par b/test/schw-0050.par index bf28d18..fb3e6ad 100644 --- a/test/schw-0050.par +++ b/test/schw-0050.par @@ -2,20 +2,20 @@ # cartoon, using BSSN_MoL #------------------------------------------------------------------------------ -ActiveThorns = "admbase admcoupling admmacros coordgauge spacemask staticconformal idanalyticbh time pugh pughreduce pughslab cartgrid3d ioascii ioutil iobasic cartoon2d bssn_mol mol admconstraints boundary nanchecker coordbase symbase localreduce" +ActiveThorns = "admbase admcoupling admmacros coordgauge spacemask staticconformal idanalyticbh time pugh pughreduce pughslab cartgrid3d ioascii ioutil iobasic cartoon2d ml_bssn ml_bssn_helper genericfd newrad loopcontrol tmunubase mol admconstraints boundary nanchecker coordbase symbase localreduce reflectionsymmetry" -IOBasic::outScalar_vars = " admbase::metric admbase::curv admbase::lapse admbase::shift staticconformal::psi adm_bssn::adm_bssn_metric adm_bssn::adm_bssn_curv adm_bssn::adm_bssn_k adm_bssn::adm_bssn_gamma adm_bssn::adm_bssn_phi admconstraints::momy admconstraints::ham " +IOBasic::outScalar_vars = " admbase::metric admbase::curv admbase::lapse admbase::shift ml_bssn::ml_metric ml_bssn::ml_curv ml_bssn::ml_trace_curv ml_bssn::ml_Gamma ml_bssn::phi admconstraints::momy admconstraints::ham " IOBasic::outInfo_vars = " admconstraints::ham admbase::alp " -IOASCII::out1D_vars = " admbase::metric admbase::curv admbase::lapse admbase::shift staticconformal::psi adm_bssn::adm_bssn_metric adm_bssn::adm_bssn_curv adm_bssn::adm_bssn_k adm_bssn::adm_bssn_gamma adm_bssn::adm_bssn_phi admconstraints::momy admconstraints::ham " +IOASCII::out1D_vars = " admbase::metric admbase::curv admbase::lapse admbase::shift ml_bssn::ml_metric ml_bssn::ml_curv ml_bssn::ml_trace_curv ml_bssn::ml_Gamma ml_bssn::phi admconstraints::momy admconstraints::ham " #============================================================================== # schwarzschild initial data #============================================================================== admbase::initial_data = "bl_bh" -admbase::metric_type = "static conformal" +admbase::metric_type = "physical" idanalyticbh::bl_M_1 = 1.0 #============================================================================== @@ -40,13 +40,30 @@ grid::xmax = 4.85 grid::ymin = -0.2 grid::ymax = 0.2 -grid::zmin = -4.85 +grid::zmin = -0.15 grid::zmax = 4.85 time::dtfac = 0.25 -grid::domain = "bitant" -grid::bitant_plane = "xy" +grid::domain = "full" + +coordbase::xmin = -0.15 +coordbase::xmax = 4.85 + +coordbase::ymin = -0.2 +coordbase::ymax = 0.2 + +coordbase::zmin = -0.15 +coordbase::zmax = 4.85 + +coordbase::boundary_size_x_upper = 2 +coordbase::boundary_size_y_upper = 2 +coordbase::boundary_size_z_upper = 2 +coordbase::boundary_size_x_lower = 2 +coordbase::boundary_size_y_lower = 2 +coordbase::boundary_size_z_lower = 2 + +reflectionsymmetry::reflection_z = yes driver::ghost_size_x = 2 driver::ghost_size_y = 2 @@ -59,7 +76,6 @@ grid::avoid_originy = "no" cartoon2d::stencil = "yes" cartoon2d::order = 3 -adm_bssn::cartoon = "yes" cartoon2d::cartoon_active = "yes" admconstraints::cartoon = "yes" @@ -67,14 +83,23 @@ admconstraints::cartoon = "yes" # Standard ADM_BSSN evolution #============================================================================== -admbase::evolution_method = "adm_bssn" +admbase::evolution_method = "ml_bssn" methodoflines::ode_method = "icn" -adm_bssn::ICN_itnum = 3 -adm_bssn::stencil_size = 2 -adm_bssn::advection = "upwind2" -adm_bssn::bound = "newrad" +ml_bssn::fdOrder = 2 +ml_bssn::my_initial_boundary_condition = "extrapolate-gammas" +ml_bssn::my_rhs_boundary_condition = "newrad" +# these are necessary +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" boundary::radpower = 3 #============================================================================== @@ -82,21 +107,23 @@ boundary::radpower = 3 #============================================================================== admbase::initial_lapse = "one" -admbase::lapse_evolution_method = "1+log" -adm_bssn::lapsesource = "modified" -adm_bssn::harmonic_f = 2.0 +admbase::lapse_evolution_method = "ML_BSSN" +admbase::dtlapse_evolution_method = "ML_BSSN" +ml_bssn::alphadriver = 1 +ml_bssn::lapseadvectioncoeff = 0 +ml_bssn::harmonicf = 2.0 #============================================================================== -# Gamma2 shift +# Gamma0 shift #============================================================================== admbase::initial_shift = "zero" -admbase::shift_evolution_method = "gamma2" - -adm_bssn::ShiftGammaCoeff = 0.75 -adm_bssn::BetaDriver = 4.0 -adm_bssn::ShiftAlpPower = 1 -adm_bssn::ShiftPsiPower = 4 +admbase::initial_dtshift = "zero" +admbase::shift_evolution_method = "ML_BSSN" +admbase::dtshift_evolution_method = "ML_BSSN" +ml_bssn::shiftgammacoeff = 0.75 +ml_bssn::betadriver = 4.0 +ml_bssn::shiftadvectioncoeff = 0 #============================================================================== # 1D output |