diff options
Diffstat (limited to 'par/MagTOVstar.par')
| -rw-r--r-- | par/MagTOVstar.par | 671 |
1 files changed, 671 insertions, 0 deletions
diff --git a/par/MagTOVstar.par b/par/MagTOVstar.par new file mode 100644 index 0000000..5621c00 --- /dev/null +++ b/par/MagTOVstar.par @@ -0,0 +1,671 @@ +#============================================================================== +# Poloidally Magnetized TOV star: mclachlan TOVSolver GRHydro_InitData +#============================================================================== + +#------------------------------------------------------------------------------ +# Cactus parameters: +#------------------------------------------------------------------------------ +Cactus::cctk_run_title = "Poloidally Magnetized TOV star" +Cactus::cctk_full_warnings = "yes" + +#Cactus::terminate = "never" +#Cactus::terminate = "time" +#Cactus::cctk_final_time = 1500.0 +Cactus::terminate = "iteration" +#Cactus::cctk_itlast = 0 +#Cactus::cctk_itlast = 256 +Cactus::cctk_itlast = 1024 +#Cactus::terminate="runtime" +#Cactus::max_runtime=28800 # 46 days +#Cactus::max_runtime=1410 # 23.5 hours +#Cactus::max_runtime=1380 # 23 hours +#Cactus::max_runtime=690 # 11.5 hours +#Cactus::max_runtime=450 # 7.5 hours +#Cactus::max_runtime=210 # 3.5 hours + +#------------------------------------------------------------------------------ +# Activate all necessary thorns: +#------------------------------------------------------------------------------ + +ActiveThorns = "Boundary CartGrid3D CoordBase Fortran InitBase IOUtil LocalReduce SymBase Time" +ActiveThorns = "AEILocalInterp LocalInterp" +ActiveThorns = "MoL ReflectionSymmetry RotatingSymmetry90 Slab SpaceMask SphericalSurface" +#ActiveThorns = "MoL ReflectionSymmetry RotatingSymmetry180 Slab SpaceMask SphericalSurface" +#ActiveThorns = "MoL ReflectionSymmetry Slab SpaceMask SphericalSurface" +ActiveThorns = "Carpet CarpetInterp CarpetIOASCII CarpetIOHDF5 CarpetIOScalar CarpetLib CarpetIOBasic CarpetReduce CarpetRegrid2 CarpetSlab CarpetTracker CarpetMask LoopControl" +#ActiveThorns = "Formaline NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "NaNChecker TerminationTrigger TimerReport" +ActiveThorns = "ADMbase ADMcoupling ADMmacros CoordGauge StaticConformal" +#ActiveThorns = "PunctureTracker" +#ActiveThorns = "ADManalysis Constants TmunuBase HydroBase " +ActiveThorns = "Constants TmunuBase HydroBase " +ActiveThorns = "QuasiLocalMeasures" +#ActiveThorns = "ADMConstraints" +#ActiveThorns = "EOS_Base EOS_IdealFluid EOS_Polytrope" +ActiveThorns = "EOS_Omni" +ActiveThorns = "GRHydro GRHydro_InitData" +ActiveThorns = "BLAS LAPACK GSL HDF5" +ActiveThorns = "SummationByParts" +ActiveThorns = "TOVSolver" +ActiveThorns = "GenericFD NewRad" +ActiveThorns = "ML_BSSN ML_BSSN_Helper ML_ADMConstraints" +ActiveThorns = "Hydro_Analysis" +ActiveThorns = "Dissipation" +ActiveThorns = "AHFinderDirect" +ActiveThorns = "WeylScal4 Multipole" +ActiveThorns = "SetMask_SphericalSurface" +ActiveThorns = "SystemStatistics" + +#------------------------------------------------------------------------------ +# Diagnostic parameters: +#------------------------------------------------------------------------------ +AHFinderDirect::verbose_level = "physics details" + +Carpet::output_timers_every = 1024 +Carpet::output_timer_tree_every = 256 +Carpet::output_initialise_timer_tree = "yes" +Carpet::storage_verbose = "no" +#Carpet::timers_verbose = "no" +#Carpet::verbose = "yes" # only for Initial Data +Carpet::verbose = "no" +Carpet::veryverbose = "no" +Carpet::grid_structure_filename = "carpet-grid-structure.asc" +Carpet::grid_coordinates_filename = "carpet-grid-coordinates.asc" + +CarpetLib::output_bboxes = "no" +CarpetLib::print_memstats_every = 1024 +CarpetLib::memstat_file = "carpetlib-memory-statistics.asc" +#CarpetLib::print_timestats_every = 1024 +#CarpetLib::timestat_file = "carpetlib-timing-statistics.asc" + +CarpetMask::verbose = "yes" +CarpetReduce::verbose = "no" +CarpetRegrid2::verbose = "yes" +#CarpetRegrid2::veryverbose = "yes" + +CarpetTracker::verbose = "yes" + +Dissipation::verbose = "no" + +#NaNChecker::verbose = "all" + +#PunctureTracker::verbose = "yes" + +TimerReport::out_every = 1024 +TimerReport::out_filename = "TimerReport" +#TimerReport::output_all_timers = "yes" +#TimerReport::output_all_timers_together = "yes" +#TimerReport::output_all_timers_readable = "yes" +#TimerReport::before_checkpoint = "yes" +TimerReport::n_top_timers = 40 + +QuasiLocalMeasures::verbose = "yes" +SphericalSurface::verbose = "yes" + +#------------------------------------------------------------------------------ +# Utility parameters: +#------------------------------------------------------------------------------ + +#NaNChecker::check_every = 256 +NaNChecker::check_every = 1 +NaNChecker::check_vars = " + ADMBase::curv + ADMBase::metric + ADMBase::lapse + ADMBase::shift + HydroBase::Bvec + HydroBase::rho + HydroBase::eps + HydroBase::press + HydroBase::vel + HydroBase::w_lorentz +" +NaNChecker::action_if_found = "terminate" +#NaNChecker::action_if_found = "abort" +#NaNChecker::action_if_found = "just warn" #"terminate", "just warn", "abort" +#NaNChecker::check_after=0 + +#TerminationTrigger::max_walltime = 0.0 # hours +TerminationTrigger::max_walltime = 24.0 # hours +TerminationTrigger::on_remaining_walltime = 10.0 # minutes +TerminationTrigger::termination_from_file = "yes" +TerminationTrigger::create_termination_file = "yes" +TerminationTrigger::termination_file = "../TERMINATE" + + +#------------------------------------------------------------------------------ +# Run parameters: +#------------------------------------------------------------------------------ + +#------ +# Grid: +#------ + +Time::dtfac = 0.25 +MoL::ODE_Method = "rk4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +CartGrid3D::type = "coordbase" +CartGrid3D::domain = "full" +CartGrid3D::avoid_origin = "no" + +CoordBase::domainsize = "minmax" +CoordBase::spacing = "gridspacing" # "gridspacing" or "numcells" + +CoordBase::xmin = 0.0 +CoordBase::ymin = 0.0 +CoordBase::zmin = 0.0 +CoordBase::xmax = 256.0 +CoordBase::ymax = 256.0 +CoordBase::zmax = 256.0 +CoordBase::dx = 8.0 # dx or ncells_x +CoordBase::dy = 8.0 # dy or ncells_y +CoordBase::dz = 8.0 # dz or ncells_z +#CoordBase::dx = 6.4 # dx or ncells_x +#CoordBase::dy = 6.4 # dy or ncells_y +#CoordBase::dz = 6.4 # dz or ncells_z +#CoordBase::dx = 5.12 # dx or ncells_x +#CoordBase::dy = 5.12 # dy or ncells_y +#CoordBase::dz = 5.12 # dz or ncells_z + +CoordBase::boundary_size_x_lower = 3 +CoordBase::boundary_size_y_lower = 3 +CoordBase::boundary_size_z_lower = 3 +CoordBase::boundary_size_x_upper = 3 +CoordBase::boundary_size_y_upper = 3 +CoordBase::boundary_size_z_upper = 3 +CoordBase::boundary_shiftout_x_lower = 1 +CoordBase::boundary_shiftout_y_lower = 1 +CoordBase::boundary_shiftout_z_lower = 1 +CoordBase::boundary_shiftout_x_upper = 0 +CoordBase::boundary_shiftout_y_upper = 0 +CoordBase::boundary_shiftout_z_upper = 0 + +Driver::ghost_size = 3 +GRHydro::GRHydro_stencil = 3 + +ReflectionSymmetry::reflection_x = "no" +ReflectionSymmetry::reflection_y = "no" +ReflectionSymmetry::reflection_z = "yes" +ReflectionSymmetry::avoid_origin_x = "no" +ReflectionSymmetry::avoid_origin_y = "no" +ReflectionSymmetry::avoid_origin_z = "no" + + +# General Carpet parameters: +Carpet::domain_from_coordbase = "yes" +Carpet::enable_all_storage = "no" +Carpet::use_buffer_zones = "yes" +Carpet::schedule_barriers = "no" + +Carpet::poison_new_timelevels = "yes" +Carpet::check_for_poison = "no" +#Carpet::poison_value = 113 # git version only! + +Carpet::init_3_timelevels = "no" +Carpet::init_fill_timelevels = "yes" + +CarpetLib::poison_new_memory = "yes" +CarpetLib::poison_value = 114 +CarpetLib::check_bboxes = "no" +CarpetLib::interleave_communications = "yes" +#CarpetLib::vary_tags = "yes" # git version only! +CarpetLib::combine_sends = "yes" +#CarpetLib::reduce_mpi_waitall = "yes" # git version only! + +CarpetInterp::tree_search = "yes" +CarpetInterp::check_tree_search = "no" + +CarpetRegrid2::freeze_unaligned_levels = "yes" +CarpetRegrid2::snap_to_coarse = "yes" +CarpetRegrid2::symmetry_rotating90 = "yes" + +# System specific Carpet parameters: +Carpet::max_refinement_levels = 9 +Carpet::prolongation_order_space = 3 +Carpet::prolongation_order_time = 2 + +CarpetRegrid2::regrid_every = 0 +CarpetRegrid2::num_centres = 1 + +CarpetRegrid2::num_levels_1 = 5 +CarpetRegrid2::position_x_1 = 0.0 +CarpetRegrid2::position_y_1 = 0.0 +CarpetRegrid2::position_z_1 = 0.0 +CarpetRegrid2::radius_1[1] =128.0 # 4 +CarpetRegrid2::radius_1[2] = 64.0 # 2 +CarpetRegrid2::radius_1[3] = 32.0 # 1 +CarpetRegrid2::radius_1[4] = 16.0 # 0.5 +CarpetRegrid2::radius_1[5] = 8.0 # 0.25 +CarpetRegrid2::radius_1[6] = 4.0 # 0.125 +CarpetRegrid2::radius_1[7] = 2.0 # 0.0625 +CarpetRegrid2::radius_1[8] = 1.0 # 0.03125 + +#CarpetMask::excluded_surface [0] = 0 +#CarpetMask::excluded_surface_factor[0] = 1.0 + +#CarpetTracker::surface[0] = 0 + +#------ +# MODEL: +#------ + +HydroBase::initial_hydro = "tov" +GRHydro_InitData::initial_Bvec = "poloidalmagfield" + +GRHydro_InitData::poloidal_A_b = 200.0 +GRHydro_InitData::poloidal_n_p = 3 +GRHydro_InitData::poloidal_P_cut = 1.0e-8 +GRHydro_InitData::poloidal_rho_max = 1.221209344680609e-03 + +ADMBase::metric_type = "physical" +ADMBase::initial_data = "tov" +ADMBase::initial_lapse = "tov" +ADMBase::initial_shift = "tov" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" +# needed for AHFinderDirect: +ADMBase::metric_timelevels = 3 + +TOVSolver::TOV_Enforce_Interpolation = "yes" +TOVSolver::TOV_Num_Radial = 1000000 +TOVSolver::TOV_dr[0] = 1.e-5 +#TOVSolver::TOV_Position_x[0] = 0.0 +TOVSolver::TOV_Rho_Central[0] = 1.221209344680609e-03 +TOVSolver::TOV_Gamma[0] = 2.0 +TOVSolver::TOV_K[0] = 1.019575500488281e+02 +#TOVSolver::TOV_Mass[0] = 1.4 +#TOVSolver::TOV_C[0] = 0.1438356164383561644 +TOVSolver::TOV_Momentum_Psi_Power = 6 + +EOS_Omni::poly_K = 1.019575500488281e+02 +EOS_Omni::poly_gamma = 2.0 +EOS_Omni::gl_K = 1.019575500488281e+02 +EOS_Omni::gl_gamma = 2.0 +GRHydro::initial_Gamma = 2.0 + +#---------- +# Numerics: +#---------- + +InitBase::initial_data_setup_method = "init_all_levels" + +TmunuBase::stress_energy_storage = "yes" +TmunuBase::stress_energy_at_RHS = "yes" +TmunuBase::timelevels = 1 +TmunuBase::prolongation_type = "none" +TmunuBase::support_old_CalcTmunu_mechanism = "no" + +HydroBase::timelevels = 3 + +ADMMacros::spatial_order = 4 +SummationByParts::order = 4 + +SpaceMask::use_mask = "yes" + +#Dissipation::ah_slope = 0.2 +#Dissipation::epsdis = 0.01 +#Dissipation::extra_dissipation_in_horizons = "yes" +#Dissipation::ah_max_epsdis=0.1 +#Dissipation::update_ah_every = 16 +#Dissipation::surface_number[0] = 0 +#Dissipation::horizon_number[0] = 1 +Dissipation::vars = " + ML_BSSN::ML_log_confac + ML_BSSN::ML_metric + ML_BSSN::ML_curv + ML_BSSN::ML_trace_curv + ML_BSSN::ML_Gamma + ML_BSSN::ML_lapse + ML_BSSN::ML_shift +" + +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 +SphericalSurface::radius [0] = 0.5 + +SetMask_SphericalSurface::SetMask_SurfaceIndex[0] = 0 +SetMask_SphericalSurface::SetMask_RadiusFactor[0] = 0.75 + + +#----------- +# Evolution: +#----------- + +HydroBase::evolution_method = "GRHydro" +HydroBase::Bvec_evolution_method = "GRHydro" + +GRHydro::method_type = "RSA FV" +#GRHydro::recon_method = "tvd" +GRHydro::recon_method = "ppm" +#GRHydro::tvd_limiter = "vanleerMC" +#GRHydro::tvd_limiter = "minmod" +GRHydro::recon_vars = "primitive" +GRHydro::riemann_solver = "HLLE" +#GRHydro::riemann_solver = "LLF" +GRHydro::bound = "none" +GRHydro::GRHydro_eos_type = "General" +GRHydro::GRHydro_eos_table = "Ideal_Fluid" + +GRHydro::Grhydro_MaxNumConstrainedVars = 13 +GRHydro::GRHydro_MaxNumEvolvedVars = 10 + +GRHydro::clean_divergence = "yes" +GRHydro::kap_dc = 1.0 +GRHydro::track_divB = "yes" +GRHydro::calculate_bcom = "yes" + +GRHydro::rho_abs_min = 1.0e-10 +#Grhydro::initial_rho_abs_min = 1.0e-13 +#Grhydro::grhydro_atmo_tolerance = 0.0001 + +GRHydro::c2p_reset_pressure = "yes" +#GRHydro::GRHydro_enable_internal_excision = "false" + +grhydro::grhydro_c2p_warnlevel = 0 +grhydro::grhydro_c2p_warn_from_reflevel = 2 +grhydro::con2prim_oct_hack = yes + + + +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::timelevels = 3 + +ML_BSSN::harmonicN = 1 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::LapseACoeff = 0.0 # don't evolve A in time (default = 0) +ML_BSSN::ShiftBCoeff = 1.0 # evolve B^i in time (default = 1) +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::AlphaDriver = 0.0 # (default = 0) +ML_BSSN::BetaDriver = 0.2 # \equiv eta; must have dt < 1/eta always. +#ML_BSSN::BetaDriver = 1.0 # Usual for BHB once dt < 1/eta satisfied. +ML_BSSN::LapseAdvectionCoeff = 1.0 +ML_BSSN::ShiftAdvectionCoeff = 1.0 + +ML_BSSN::UseSpatialBetaDriver = "no" +ML_BSSN::SpatialBetaDriverRadius = 100.0e0 + +ML_BSSN::MinimumLapse = 1.0e-8 + +ML_BSSN::my_initial_boundary_condition = "extrapolate-gammas" +ML_BSSN::my_rhs_boundary_condition = "NewRad" +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" + + +#------------------------------------------------------------------------------ +# Output: +#------------------------------------------------------------------------------ + +IO::out_dir = $parfile +IO::out_fileinfo = "none" + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_reductions = "maximum" +IOBasic::outInfo_vars = " + Carpet::physical_time_per_hour + GRHydro::bcom_sq + HydroBase::rho + HydroBase::press + ML_ADMConstraints::ML_Ham + SystemStatistics::maxrss_mb +" + +IOScalar::outScalar_every = 1024 +IOScalar::one_file_per_group = "yes" +IOScalar::outScalar_reductions = "minimum maximum average norm1 norm2" +IOScalar::outScalar_vars = " + ADMBase::curv + ADMBase::lapse + ADMBase::metric + ADMBase::shift + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GrHydro::divB + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_mom + ML_BSSN::ML_Ham + ML_BSSN::ML_mom + SystemStatistics::process_memory_mb + SphericalSurface::sf_radius +" + +IOASCII::one_file_per_group = "yes" +IOASCII::output_symmetry_points = "no" +IOASCII::output_ghost_points = "no" #hg version only! +IOASCII::output_boundary_points = "no" #hg version only! +#IOASCII::out3D_ghosts = "no" #hg (DEPRECATED) +#IOASCII::out3D_outer_ghosts = "no" #hg (DEPRECATED) + +IOASCII::out0D_every = 256 +IOASCII::out0D_vars = " + Carpet::timing + Hydro_Analysis::Hydro_Analysis_rho_max_loc +# PunctureTracker::pt_loc + QuasiLocalMeasures::qlm_scalars + SphericalSurface::sf_active + SphericalSurface::sf_coordinate_descriptors + SphericalSurface::sf_info + SphericalSurface::sf_origin + SphericalSurface::sf_radius + SphericalSurface::sf_valid +" + +#Set these IOASCII options for initial data only: +IOASCII::out1D_every = 1 +#IOASCII::out1D_dt = 10 +IOASCII::out1D_d = "yes" +IOASCII::out1D_vars = " + ADMBase::curv + ADMBase::lapse + ADMBase::metric + ADMBase::shift + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GrHydro::divB + GRHydro::GRHydro_C2P_failed + GRHydro::scon + GRHydro::tau + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz + ML_BSSN::ML_Ham + ML_BSSN::ML_mom + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_mom +" + +IOASCII::out2D_every = 256000000 +IOASCII::out2D_vars = " + HydroBase::rho + HydroBase::vel + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_mom + SphericalSurface::sf_radius{out_every=256} +" + +CarpetIOHDF5::one_file_per_group = "yes" +CarpetIOHDF5::open_one_input_file_at_a_time = "yes" + +#CarpetIOHDF5::output_symmetry_points = "no" +#CarpetIOHDF5::output_ghost_points = "no" +#CarpetIOHDF5::output_boundary_points = "no" +#CarpetIOHDF5::output_buffer_points = "no" + +#CarpetIOHDF5::out2D_every = 256 +CarpetIOHDF5::out2D_every = 1 +CarpetIOHDF5::out2D_xy = "yes" +CarpetIOHDF5::out2D_xz = "no" +CarpetIOHDF5::out2D_yz = "no" +CarpetIOHDF5::out2D_xyplane_z = 0.0 +CarpetIOHDF5::out2D_vars = " + ADMBase::curv + ADMBase::lapse + ADMBase::metric + ADMBase::shift + CarpetReduce::weight + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::Bcons + GRHydro::dens + GRHydro::divB + GRHydro::GRHydro_C2P_failed + GRHydro::scon + GRHydro::tau + Grid::coordinates + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_Mom + WeylScal4::Psi4i + WeylScal4::Psi4r +" + +#CarpetIOHDF5::out3D_every = 4096 +CarpetIOHDF5::out3D_every = 1 +CarpetIOHDF5::out3D_vars = " + ADMBase::lapse + CarpetReduce::weight + GRHydro::bcom + GRHydro::bcom0 + GRHydro::bcom_sq + GRHydro::divB + GRHydro::GRHydro_C2P_failed + HydroBase::Bvec + HydroBase::eps + HydroBase::press + HydroBase::rho + HydroBase::vel + HydroBase::w_lorentz + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_Mom +" + +#------------------------------------------------------------------------------ +# Analysis: +#------------------------------------------------------------------------------ +AHFinderDirect::find_every = 512 + +#AHFinderDirect::run_at_CCTK_ANALYSIS = "yes" +#AHFinderDirect::run_at_CCTK_POSTSTEP = "no" +AHFinderDirect::run_at_CCTK_POST_RECOVER_VARIABLES = "no" + +AHFinderDirect::move_origins = "yes" +#AHFinderDirect::reshape_while_moving = "yes" +#AHFinderDirect::predict_origin_movement = "yes" + +# Hermite to order 3 to avoid discontinuities in the metric spatial derivatives: +#AHFinderDirect::geometry_interpolator_name = "Hermite polynomial interpolation" +#AHFinderDirect::geometry_interpolator_pars = "order=3" +#AHFinderDirect::surface_interpolator_name = "Hermite polynomial interpolation" +#AHFinderDirect::surface_interpolator_pars = "order=3" +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::output_h_every = 0 + +AHFinderDirect::N_horizons = 1 + +AHFinderDirect::reset_horizon_after_not_finding [1] = "no" +AHFinderDirect::initial_guess__coord_sphere__radius [1] = 0.2 +AHFinderDirect::origin_x [1] = 0.0 +AHFinderDirect::origin_y [1] = 0.0 +AHFinderDirect::origin_z [1] = 0.0 +AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0.0 +AHFinderDirect::initial_guess__coord_sphere__y_center[1] = 0.0 +AHFinderDirect::initial_guess__coord_sphere__z_center[1] = 0.0 +#AHFinderDirect::track_origin_from_grid_scalar [1] = "yes" +#AHFinderDirect::track_origin_source_x [1] = "PunctureTracker::pt_loc_x[0]" +#AHFinderDirect::track_origin_source_y [1] = "PunctureTracker::pt_loc_y[0]" +#AHFinderDirect::track_origin_source_z [1] = "PunctureTracker::pt_loc_z[0]" +AHFinderDirect::which_surface_to_store_info [1] = 0 + + +Hydro_Analysis::Hydro_Analysis_comp_rho_max = "true" + +#PunctureTracker::track [0] = "yes" +#PunctureTracker::initial_x [0] = 10.64502257 +#PunctureTracker::initial_y [0] = 0 +#PunctureTracker::initial_z [0] = 0 +#PunctureTracker::which_surface_to_store_info[0] = 0 + +QuasiLocalMeasures::num_surfaces = 1 +QuasiLocalMeasures::spatial_order = 4 +QuasiLocalMeasures::interpolator = "Lagrange polynomial interpolation" +QuasiLocalMeasures::interpolator_options = "order=4" +QuasiLocalMeasures::surface_index [0] = 0 + +Multipole::nradii = 4 +Multipole::radius[0] = 80 +Multipole::radius[1] = 100 +Multipole::radius[2] = 120 +Multipole::radius[3] = 140 +Multipole::ntheta = 120 +Multipole::nphi = 240 +Multipole::variables = "WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i' name='psi4'}" +Multipole::out_every = 4 +Multipole::l_max = 4 +Multipole::integration_method = "Simpson" + +WeylScal4::fd_order = "4th" + +#------------------------------------------------------------------------------ +# Checkpoint/Recovery: +#------------------------------------------------------------------------------ +#IOHDF5::checkpoint = "yes" +#IO::checkpoint_dir = $parfile +#IO::checkpoint_ID = "yes" +#IO::checkpoint_every = 6144 +##IO::checkpoint_every_walltime_hours = 6.0 +#IO::checkpoint_keep=3 +#IO::checkpoint_on_terminate = "yes" +# +#IO::recover = "autoprobe" +#IO::recover_dir = $parfile + + |