User: knarf Date: 2011/01/31 08:56 AM

Added: /examples/bbh/ bbh.par

Log: first version of bbh parameter file

File Changes:

Directory: /examples/bbh/ =========================

File [added]: bbh.par Delta lines: +502 -0 =================================================================== --- examples/bbh/bbh.par (rev 0) +++ examples/bbh/bbh.par 2011-01-31 14:56:58 UTC (rev 23) @@ -0,0 +1,502 @@ +Cactus::cctk_run_title = "BBH" + +Cactus::cctk_full_warnings = yes +Cactus::highlight_warning_messages = no + +Cactus::terminate = "time" +Cactus::cctk_final_time = 300.0 + + + +ActiveThorns = "IOUtil" + +IO::out_dir = $parfile + + + +ActiveThorns = "AEILocalInterp" + +ActiveThorns = "BLAS LAPACK" + +ActiveThorns = "Fortran" + +ActiveThorns = "GSL" + +ActiveThorns = "GenericFD" + +ActiveThorns = "HDF5" + +ActiveThorns = "LocalInterp" + +ActiveThorns = "LoopControl" + +ActiveThorns = "Slab" + + + +ActiveThorns = "SummationByParts" + +SummationByParts::order = 4 + + + +ActiveThorns = "InitBase" + + + +ActiveThorns = "Carpet CarpetLib CarpetInterp CarpetReduce CarpetSlab" + +Carpet::verbose = no +Carpet::veryverbose = no +Carpet::schedule_barriers = no +Carpet::storage_verbose = no +#Carpet::timers_verbose = no +CarpetLib::output_bboxes = no + +Carpet::domain_from_coordbase = yes +Carpet::max_refinement_levels = 7 + +driver::ghost_size = 3 +Carpet::use_buffer_zones = yes + +Carpet::prolongation_order_space = 5 +Carpet::prolongation_order_time = 2 + +Carpet::convergence_level = 0 + +Carpet::init_fill_timelevels = yes +#Carpet::init_3_timelevels = yes + +Carpet::poison_new_timelevels = yes +CarpetLib::poison_new_memory = yes + +Carpet::output_timers_every = 640 +CarpetLib::print_timestats_every = 640 +CarpetLib::print_memstats_every = 640 + + + +ActiveThorns = "NaNChecker" + +NaNChecker::check_every = 1 # 64 +#NaNChecker::verbose = "all" +#NaNChecker::action_if_found = "just warn" +NaNChecker::action_if_found = "terminate" +NaNChecker::check_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 + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift +" + + + +ActiveThorns = "Boundary CartGrid3D CoordBase ReflectionSymmetry RotatingSymmetry180 SymBase" + +CoordBase::domainsize = "minmax" + +CoordBase::xmin = 0.00 +CoordBase::ymin = -120.00 +CoordBase::zmin = 0.00 +CoordBase::xmax = +120.00 +CoordBase::ymax = +120.00 +CoordBase::zmax = +120.00 +CoordBase::dx = 2.00 +CoordBase::dy = 2.00 +CoordBase::dz = 2.00 + +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_z_lower = 1 + +CartGrid3D::type = "coordbase" + +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_z = no + + + +ActiveThorns = "SphericalSurface" + +SphericalSurface::nsurfaces = 5 +SphericalSurface::maxntheta = 39 +SphericalSurface::maxnphi = 76 + +SphericalSurface::ntheta [0] = 39 +SphericalSurface::nphi [0] = 76 +SphericalSurface::nghoststheta[0] = 2 +SphericalSurface::nghostsphi [0] = 2 + +SphericalSurface::ntheta [1] = 39 +SphericalSurface::nphi [1] = 76 +SphericalSurface::nghoststheta[1] = 2 +SphericalSurface::nghostsphi [1] = 2 + +SphericalSurface::ntheta [2] = 39 +SphericalSurface::nphi [2] = 76 +SphericalSurface::nghoststheta[2] = 2 +SphericalSurface::nghostsphi [2] = 2 + +SphericalSurface::ntheta [3] = 39 +SphericalSurface::nphi [3] = 76 +SphericalSurface::nghoststheta[3] = 2 +SphericalSurface::nghostsphi [3] = 2 + +SphericalSurface::ntheta [4] = 39 +SphericalSurface::nphi [4] = 76 +SphericalSurface::nghoststheta[4] = 2 +SphericalSurface::nghostsphi [4] = 2 + + + +ActiveThorns = "CarpetMask" + +CarpetMask::verbose = no + +CarpetMask::excluded_surface [0] = 3 +CarpetMask::excluded_surface_factor[0] = 1.0 + +CarpetMask::excluded_surface [1] = 4 +CarpetMask::excluded_surface_factor[1] = 1.0 + + + +ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal TmunuBase" +ActiveThorns = "CarpetRegrid2 PunctureTracker CarpetTracker" + +CarpetTracker::surface[0] = 0 +CarpetTracker::surface[1] = 1 +PunctureTracker::track [0] = yes +PunctureTracker::initial_x [0] = 3.0 +PunctureTracker::which_surface_to_store_info[0] = 0 +PunctureTracker::track [1] = yes +PunctureTracker::initial_x [1] = -3.0 +PunctureTracker::which_surface_to_store_info[1] = 1 + + +CarpetRegrid2::regrid_every = 16 +CarpetRegrid2::freeze_unaligned_levels = yes +CarpetRegrid2::symmetry_rotating180 = yes +CarpetRegrid2::verbose = yes + +CarpetRegrid2::num_centres = 2 + +CarpetRegrid2::num_levels_1 = 7 +CarpetRegrid2::position_x_1 = +3.0 +CarpetRegrid2::radius_1[ 1] = 64.0 +CarpetRegrid2::radius_1[ 2] = 16.0 +CarpetRegrid2::radius_1[ 3] = 8.0 +CarpetRegrid2::radius_1[ 4] = 4.0 +CarpetRegrid2::radius_1[ 5] = 2.0 +CarpetRegrid2::radius_1[ 6] = 1.0 +CarpetRegrid2::movement_threshold_1 = 0.16 + +CarpetRegrid2::num_levels_2 = 7 +CarpetRegrid2::position_x_2 = -3.0 +CarpetRegrid2::radius_2[ 1] = 64.0 +CarpetRegrid2::radius_2[ 2] = 16.0 +CarpetRegrid2::radius_2[ 3] = 8.0 +CarpetRegrid2::radius_2[ 4] = 4.0 +CarpetRegrid2::radius_2[ 5] = 2.0 +CarpetRegrid2::radius_2[ 6] = 1.0 +CarpetRegrid2::movement_threshold_2 = 0.16 + + +ActiveThorns = "MoL Time" + +MoL::ODE_Method = "RK4" +MoL::MoL_Intermediate_Steps = 4 +MoL::MoL_Num_Scratch_Levels = 1 + +Carpet::time_refinement_factors = "[1, 1, 2, 4, 8, 16, 32, 64, 128, 256]" + +Time::dtfac = 0.25 + + +ADMMacros::spatial_order = 4 + + + +ActiveThorns = "TwoPunctures" + +ADMBase::metric_type = "physical" + +ADMBase::initial_data = "twopunctures" +ADMBase::initial_lapse = "twopunctures-averaged" +ADMBase::initial_shift = "zero" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_dtshift = "zero" + +TwoPunctures::par_b = 3.0 +TwoPunctures::par_m_plus = 0.47656 +TwoPunctures::par_m_minus = 0.47656 +TwoPunctures::par_P_plus [1] = +0.13808 +TwoPunctures::par_P_minus[1] = -0.13808 + +#TODO# TwoPunctures::grid_setup_method = "evaluation" + +TwoPunctures::TP_epsilon = 1.0e-6 +TwoPunctures::TP_Tiny = 1.0e-2 + +TwoPunctures::verbose = yes + + + +ActiveThorns = "ML_BSSN ML_BSSN_Helper NewRad" + +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 = 1 # 1+log +ML_BSSN::harmonicF = 2.0 # 1+log +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::BetaDriver = 1.0 +ML_BSSN::LapseAdvectionCoeff = 1.0 +ML_BSSN::ShiftAdvectionCoeff = 1.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 + +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" + + + +ActiveThorns = "Dissipation" + +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_shift + ML_BSSN::ML_dtlapse + ML_BSSN::ML_dtshift +" + + + +ActiveThorns = "ML_ADMConstraints" + + + +ActiveThorns = "AHFinderDirect" + +AHFinderDirect::find_every = 16 + +AHFinderDirect::run_at_CCTK_ANALYSIS = yes +AHFinderDirect::run_at_CCTK_POSTSTEP = false + +AHFinderDirect::move_origins = yes + +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 = 3 + +AHFinderDirect::origin_x [1] = +3.0 +AHFinderDirect::initial_guess__coord_sphere__x_center[1] = +3.0 +AHFinderDirect::initial_guess__coord_sphere__radius [1] = 0.25 +AHFinderDirect::which_surface_to_store_info [1] = 2 +AHFinderDirect::set_mask_for_individual_horizon [1] = no +AHFinderDirect::reset_horizon_after_not_finding [1] = no +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::max_allowable_horizon_radius [1] = 3 + +AHFinderDirect::origin_x [2] = -3.0 +AHFinderDirect::initial_guess__coord_sphere__x_center[2] = -3.0 +AHFinderDirect::initial_guess__coord_sphere__radius [2] = 0.25 +AHFinderDirect::which_surface_to_store_info [2] = 3 +AHFinderDirect::set_mask_for_individual_horizon [2] = no +AHFinderDirect::reset_horizon_after_not_finding [2] = no +AHFinderDirect::track_origin_from_grid_scalar [2] = yes +AHFinderDirect::track_origin_source_x [2] = "PunctureTracker::pt_loc_x[1]" +AHFinderDirect::track_origin_source_y [2] = "PunctureTracker::pt_loc_y[1]" +AHFinderDirect::track_origin_source_z [2] = "PunctureTracker::pt_loc_z[1]" +AHFinderDirect::max_allowable_horizon_radius [2] = 3 + +AHFinderDirect::origin_x [3] = 0 +AHFinderDirect::find_after_individual_time [3] = 100.0 +AHFinderDirect::initial_guess__coord_sphere__x_center[3] = 0 +AHFinderDirect::initial_guess__coord_sphere__radius [3] = 1.0 +AHFinderDirect::which_surface_to_store_info [3] = 4 +AHFinderDirect::reset_horizon_after_not_finding [3] = no +AHFinderDirect::max_allowable_horizon_radius [3] = 6 + + +ActiveThorns = "WeylScal4 Multipole" + +Multipole::nradii = 4 +Multipole::radius[0] = 30 +Multipole::radius[1] = 40 +Multipole::radius[2] = 50 +Multipole::radius[3] = 60 +Multipole::ntheta = 60 +Multipole::nphi = 120 +Multipole::variables = "WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i' name='psi4'}" +Multipole::out_every = 4 +Multipole::l_max = 8 + + +ActiveThorns = "CarpetIOBasic" + +IOBasic::outInfo_every = 1 +IOBasic::outInfo_reductions = "norm2" +IOBasic::outInfo_vars = " + Carpet::physical_time_per_hour + ML_BSSN::ML_trace_curv +" + + + +ActiveThorns = "CarpetIOScalar" + +IOScalar::one_file_per_group = yes + +IOScalar::outScalar_every = 16 +IOScalar::outScalar_vars = " + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_mom + SphericalSurface::sf_radius +" + + + +ActiveThorns = "CarpetIOASCII" + +IOASCII::one_file_per_group = yes + +IOASCII::output_symmetry_points = no +IOASCII::out3D_ghosts = no + +IOASCII::out0D_every = 16 +IOASCII::out0D_vars = " + Carpet::timing + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_mom + SphericalSurface::sf_active + SphericalSurface::sf_valid + SphericalSurface::sf_info + SphericalSurface::sf_radius + SphericalSurface::sf_origin + SphericalSurface::sf_coordinate_descriptors + PunctureTracker::pt_loc +" + +IOASCII::out1D_every = 16 +IOASCII::out1D_vars = " + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ADMBase::dtlapse + ADMBase::dtshift + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_mom + SphericalSurface::sf_radius +" + +IOASCII::out2D_every = 16 +IOASCII::out2D_vars = " + SphericalSurface::sf_radius +" + + + +Activethorns = "CarpetIOHDF5" + +IOHDF5::out_every = 64 +IOHDF5::one_file_per_group = yes +IOHDF5::output_symmetry_points = no +IOHDF5::out3D_ghosts = no +IOHDF5::compression_level = 1 +IOHDF5::use_checksums = yes +IOHDF5::out_vars = " + CarpetReduce::weight + ADMBase::metric + ADMBase::curv + ADMBase::lapse + ADMBase::shift + ML_ADMConstraints::ML_Ham + ML_ADMConstraints::ML_mom + WeylScal4::Psi4r + WeylScal4::Psi4i +" + +IOHDF5::checkpoint = yes +IO::checkpoint_dir = $parfile +IO::checkpoint_ID = yes +IO::checkpoint_every_walltime_hours = 6.0 +IO::checkpoint_on_terminate = yes + +IO::recover = "autoprobe" +IO::recover_dir = $parfile + + + +ActiveThorns = "Formaline" + + + +ActiveThorns = "TimerReport" + +TimerReport::out_every = 640 +TimerReport::out_filename = "TimerReport" +TimerReport::output_all_timers_together = yes +TimerReport::output_all_timers_readable = yes +TimerReport::n_top_timers = 20