From 95c73d61458dc832692f575d8967a0f9884068c2 Mon Sep 17 00:00:00 2001 From: Erik Schnetter Date: Wed, 5 Dec 2007 20:31:26 -0600 Subject: Make single black hole evolution work --- par/ks_bssn_evol.par | 186 ++++++++++++++++++++++++++++++++++----------------- 1 file changed, 125 insertions(+), 61 deletions(-) (limited to 'par') diff --git a/par/ks_bssn_evol.par b/par/ks_bssn_evol.par index b7b1810..a9118ca 100644 --- a/par/ks_bssn_evol.par +++ b/par/ks_bssn_evol.par @@ -5,12 +5,10 @@ ActiveThorns = " IOUtil InitBase MoL - time SymBase + Time - ADMAnalysis ADMBase - ADMConstraints ADMCoupling ADMMacros CoordGauge @@ -23,27 +21,29 @@ ActiveThorns = " CarpetIOASCII CarpetIOBasic CarpetIOScalar + CarpetInterp CarpetLib CarpetReduce CarpetRegrid LoopControl + IsolatedHorizon ReflectionSymmetry + AEILocalInterp + AHFinderDirect + Dissipation Exact Formaline GSL + SphericalSurface Slab GenericFD ML_BSSN - SphericalSurface TmunuBase - Dissipation - - nanchecker " # output @@ -57,91 +57,130 @@ IOBasic::outInfo_vars = " " IOScalar::one_file_per_group = yes -IOScalar::outScalar_every = 1 +IOScalar::outScalar_every = 32 IOScalar::outScalar_vars = " ADMBase::metric ADMBase::curv ADMBase::lapse ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma ML_BSSN::Ham ML_BSSN::mom ML_BSSN::cons_detg ML_BSSN::cons_Gamma ML_BSSN::cons_traceA " +# 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_shift +# ML_BSSN::ML_dtshift +# ML_BSSN::ML_log_confacrhs +# ML_BSSN::ML_metricrhs +# ML_BSSN::ML_Gammarhs +# ML_BSSN::ML_trace_curvrhs +# ML_BSSN::ML_curvrhs +# ML_BSSN::ML_curvrhs +# ML_BSSN::ML_lapserhs +# ML_BSSN::ML_dtlapserhs +# ML_BSSN::ML_shiftrhs +# ML_BSSN::ML_dtshiftrhs IOASCII::one_file_per_group = yes +IOASCII::out0D_every = 32 +IOASCII::out0D_vars = " + IsolatedHorizon::ih_shapes + IsolatedHorizon::ih_state + IsolatedHorizon::ih_grid_int + IsolatedHorizon::ih_grid_real + IsolatedHorizon::ih_shapes + IsolatedHorizon::ih_coordinates + IsolatedHorizon::ih_tetrad_l + IsolatedHorizon::ih_tetrad_n + IsolatedHorizon::ih_tetrad_m + IsolatedHorizon::ih_newman_penrose + IsolatedHorizon::ih_weyl_scalars + IsolatedHorizon::ih_twometric + IsolatedHorizon::ih_killing_vector + IsolatedHorizon::ih_killed_twometric + IsolatedHorizon::ih_scalars + IsolatedHorizon::ih_invariant_coordinates + IsolatedHorizon::ih_multipole_moments + IsolatedHorizon::ih_3determinant +" + IOASCII::out1D_every = 32 IOASCII::out1D_vars = " ADMBase::metric ADMBase::curv ADMBase::lapse ADMBase::shift - ADMAnalysis::trace_of_K - ADMAnalysis::ricci_scalar - ADMConstraints::Hamiltonian - ADMConstraints::momentum - ML_BSSN::ML_log_confac - ML_BSSN::ML_metric - ML_BSSN::ML_trace_curv - ML_BSSN::ML_curv - ML_BSSN::ML_Gamma ML_BSSN::Ham ML_BSSN::mom ML_BSSN::cons_detg ML_BSSN::cons_Gamma ML_BSSN::cons_traceA - ML_BSSN::ML_log_confacrhs - ML_BSSN::ML_metricrhs - ML_BSSN::ML_Gammarhs - ML_BSSN::ML_trace_curvrhs - ML_BSSN::ML_curvrhs + IsolatedHorizon::ih_shapes + IsolatedHorizon::ih_newman_penrose + IsolatedHorizon::ih_weyl_scalars + IsolatedHorizon::ih_killing_vector + IsolatedHorizon::ih_killed_twometric + IsolatedHorizon::ih_3determinant " +# 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_shift +# ML_BSSN::ML_dtshift +# ML_BSSN::ML_log_confacrhs +# ML_BSSN::ML_metricrhs +# ML_BSSN::ML_Gammarhs +# ML_BSSN::ML_trace_curvrhs +# ML_BSSN::ML_curvrhs +# ML_BSSN::ML_curvrhs +# ML_BSSN::ML_lapserhs +# ML_BSSN::ML_dtlapserhs +# ML_BSSN::ML_shiftrhs +# ML_BSSN::ML_dtshiftrhs #--- driver +Carpet::init_fill_timelevels = yes + Cactus::cctk_itlast = 1024 Carpet::use_buffer_zones = yes +Carpet::use_tapered_grids = yes -CarpetLib::interleave_communications = yes -CarpetLib::vary_tags = yes -CarpetLib::combine_sends = yes -CarpetLib::reduce_mpi_waitall = yes CarpetLib::print_memstats_every = 1024 -Carpet::regrid_in_level_mode = no -Carpet::output_timers_every = 1024 CarpetLib::print_timestats_every = 1024 -Carpet::print_timestats_every = 1 +Carpet::output_timers_every = 1024 +Carpet::print_timestats_every = 1024 -carpet::verbose = no +carpet::verbose = yes carpet::veryverbose = no #--- Carpet prolongation order settings -Carpet::prolongation_order_space = 3 -Carpet::prolongation_order_time = 2 +Carpet::prolongation_order_space = 5 +Carpet::prolongation_order_time = 0 + # MoL time integration -MoL::ODE_Method = RK4 -Carpet::num_integrator_substeps = 4 +MoL::ODE_Method = "RK4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 -#mol::ode_method = "generic" -#mol::generic_type = "RK" -MoL::MoL_Intermediate_Steps = 4 -MoL::MoL_Num_Scratch_Levels = 1 -time::dtfac = 0.25 +Time::dtfac = 0.25 #--- grid and symmetries @@ -196,11 +235,8 @@ carpetregrid::refinement_levels = 7 ADMBase::initial_data = "exact" ADMBase::initial_lapse = "exact" ADMBase::initial_shift = "exact" -ADMBase::initial_dtlapse = "zero" -ADMBase::initial_dtshift = "zero" -ADMBase::lapse_timelevels = 3 -ADMBase::shift_timelevels = 3 -ADMBase::metric_timelevels = 3 +ADMBase::initial_dtlapse = "exact" +ADMBase::initial_dtshift = "exact" Exact::exact_model = "Kerr/Kerr-Schild" Exact::Kerr_KerrSchild__mass = 1.0 @@ -213,7 +249,7 @@ Exact::Kerr_KerrSchild__epsilon = 1.0 #--- ML_BSSN -ML_BSSN::harmonicN = 0 # 1+log +ML_BSSN::harmonicN = 1 # 1+log ML_BSSN::harmonicF = 2.0 # 1+log ML_BSSN::ShiftGammaCoeff = 0.75 ML_BSSN::BetaDriver = 0.5 @@ -228,11 +264,11 @@ ML_BSSN::ML_dtlapse_bound = "radiative" ML_BSSN::ML_shift_bound = "radiative" ML_BSSN::ML_dtshift_bound = "radiative" -ADMMacros::spatial_order = 4 - #--- Dissipation -Dissipation::vars = " +Dissipation::order = 5 +Dissipation::epsdis = 0.3 +Dissipation::vars = " ML_BSSN::ML_Gamma ML_BSSN::ML_lapse ML_BSSN::ML_shift @@ -244,11 +280,39 @@ Dissipation::vars = " #--- Analysis -admconstraints::constraints_timelevels = 3 -admconstraints::constraints_persist = yes -admconstraints::bound = "static" - nanchecker::check_every = 1 nanchecker::check_vars = "admbase::gxx" nanchecker::action_if_found = "terminate" +#--- Horizons + +SphericalSurface::nsurfaces = 1 +SphericalSurface::maxntheta = 39 +SphericalSurface::maxnphi = 76 + +SphericalSurface::ntheta [0] = 39 +SphericalSurface::nphi [0] = 76 +SphericalSurface::nghoststheta[0] = 2 +SphericalSurface::nghostsphi [0] = 2 + +AHFinderDirect::find_every = 32 +AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = no + +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::N_horizons = 1 + +AHFinderDirect::initial_guess__coord_sphere__radius [1] = 2.0 +AHFinderDirect::which_surface_to_store_info [1] = 0 +AHFinderDirect::reset_horizon_after_not_finding [1] = no + +IsolatedHorizon::verbose = yes +IsolatedHorizon::interpolator = "Lagrange polynomial interpolation" +IsolatedHorizon::interpolator_options = "order=4" +IsolatedHorizon::spatial_order = 4 + +IsolatedHorizon::num_horizons = 1 +IsolatedHorizon::surface_index[0] = 0 -- cgit v1.2.3