From 5121d3bcf501be5532e62b8c06c9a025e9d532c0 Mon Sep 17 00:00:00 2001 From: Anton Khirnov Date: Wed, 23 Oct 2013 09:37:32 +0200 Subject: Initial commit. --- contrib/skel.par | 350 +++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 350 insertions(+) create mode 100644 contrib/skel.par (limited to 'contrib/skel.par') diff --git a/contrib/skel.par b/contrib/skel.par new file mode 100644 index 0000000..59793b5 --- /dev/null +++ b/contrib/skel.par @@ -0,0 +1,350 @@ +############################################################################## +# GRID SETUP # +############################################################################## +ActiveThorns = "CoordBase" + +CoordBase::domainsize = "minmax" +CoordBase::xmin = -256 +CoordBase::xmax = 256 +CoordBase::ymin = 0 +CoordBase::ymax = 256 +CoordBase::zmin = 0 +CoordBase::zmax = 256 + +CoordBase::spacing = "numcells" +CoordBase::ncells_x = 64 +CoordBase::ncells_y = 32 +CoordBase::ncells_z = 32 + +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 + +ActiveThorns = "CartGrid3d" + +CartGrid3d::type = "coordbase" + +ActiveThorns = "Carpet CarpetLib CarpetRegrid2 CarpetTracker SphericalSurface IOUtil" + +Carpet::verbose = 0 +Carpet::domain_from_coordbase = 1 +Carpet::init_fill_timelevels = 1 +Carpet::poison_new_timelevels = 1 +Carpet::max_refinement_levels = 10 +Carpet::refinement_factor = 2 # anything else not implemented +Carpet::refinement_centering = "vertex" +Carpet::ghost_size = 3 +Carpet::use_buffer_zones = 1 +Carpet::prolongation_order_space = 5 +Carpet::prolongation_order_time = 2 +Carpet::num_threads = 4 +CarpetLib::poison_new_memory = 1 + +CarpetRegrid2::verbose = 1 +CarpetRegrid2::regrid_every = 256 +CarpetRegrid2::num_centres = 1 +CarpetRegrid2::movement_threshold_1 = 0.1 +CarpetRegrid2::position_x_1 = 0 +CarpetRegrid2::num_levels_1 = 8 +CarpetRegrid2::radius_1[1] = 128 +CarpetRegrid2::radius_1[2] = 64 +CarpetRegrid2::radius_1[3] = 32 +CarpetRegrid2::radius_1[4] = 16 +CarpetRegrid2::radius_1[5] = 8 +CarpetRegrid2::radius_1[6] = 4 +CarpetRegrid2::radius_1[7] = 2 +CarpetRegrid2::radius_1[8] = 1 +CarpetRegrid2::radius_1[9] = 0.5 + +CarpetTracker::verbose = 0 +CarpetTracker::surface[0] = 1 + +ActiveThorns = "ReflectionSymmetry" + +ReflectionSymmetry::reflection_z = 1 +ReflectionSymmetry::avoid_origin_z = 0 +ReflectionSymmetry::reflection_y = 1 +ReflectionSymmetry::avoid_origin_y = 0 +#ReflectionSymmetry::reflection_x = 1 +#ReflectionSymmetry::avoid_origin_x = 0 + + +############################################################################## +# INITIAL DATA # +############################################################################## + +ActiveThorns = "ADMBase ADMCoupling ADMMacros StaticConformal Trumpet" + +ADMBase::metric_timelevels = 3 +ADMBase::initial_data = "trumpet" +ADMBase::initial_lapse = "trumpet" +ADMBase::initial_dtlapse = "zero" +ADMBase::initial_shift = "trumpet" +ADMBase::initial_dtshift = "zero" + +ADMMacros::spatial_order = 4 + +Trumpet::boost_velocity = 0.5 + +############################################################################## +# EVOLUTION # +############################################################################## + +Cactus::terminate = "time" +Cactus::cctk_final_time = 100 + +ActiveThorns = "MoL Time" + +MoL::initial_data_is_crap = 0 +MoL::ODE_Method = "RK3" +MoL::MoL_Intermediate_Steps = 3 +MoL::MoL_Num_Scratch_Levels = 1 +MoL::MoL_NaN_Check = 1 + +Time::dtfac = 0.125 + +ActiveThorns = "ML_BSSN ML_BSSN_Helper CoordGauge GenericFD NewRad Boundary StaticConformal TmunuBase + SymBase LoopControl" + +ADMBase::evolution_method = "ML_BSSN" +ADMBase::lapse_evolution_method = "ML_BSSN" +ADMBase::dtlapse_evolution_method = "ML_BSSN" +ADMBase::shift_evolution_method = "ML_BSSN" +ADMBase::dtshift_evolution_method = "ML_BSSN" + +ML_BSSN::verbose = 1 +ML_BSSN::timelevels = 3 +ML_BSSN::fdOrder = 4 +ML_BSSN::harmonicF = 2 # \dot{α} = - f α^n K +ML_BSSN::harmonicN = 1 # 1+log +ML_BSSN::conformalMethod = 1 # χ method +ML_BSSN::ShiftGammaCoeff = 0.75 +ML_BSSN::HarmonicShift = 0 +ML_BSSN::LapseACoeff = 0.0 +ML_BSSN::ShiftBCoeff = 1.0 +ML_BSSN::AlphaDriver = 0.0 +ML_BSSN::BetaDriver = 2.0 +ML_BSSN::LapseAdvectionCoeff = 1.0 +ML_BSSN::ShiftAdvectionCoeff = 1.0 +ML_BSSN::MinimumLapse = 1e-08 +ML_BSSN::dt_lapse_shift_method = "correct" +ML_BSSN::EpsDiss = 0 +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 SpaceMask" + +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 = "PunctureTracker" + +PunctureTracker::track[0] = 1 +PunctureTracker::initial_x[0] = 0 +PunctureTracker::which_surface_to_store_info[0] = 1 + +############################################################################### +# ANALYSIS # +############################################################################### + +ActiveThorns = "NaNChecker CarpetReduce" + +NaNChecker::check_every = 1024 +NaNChecker::action_if_found = "just warn" +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 +" + +SphericalSurface::nsurfaces = 6 +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 +SphericalSurface::ntheta [5] = 39 +SphericalSurface::nphi [5] = 76 +SphericalSurface::nghoststheta[5] = 2 +SphericalSurface::nghostsphi [5] = 2 + +ActiveThorns = "CarpetInterp AEILocalInterp" +ActiveThorns = "AHFinderDirect CarpetInterp AEILocalInterp" + +AHFinderDirect::find_every = 256 +AHFinderDirect::move_origins = 1 +AHFinderDirect::predict_origin_movement = 0 +AHFinderDirect::N_horizons = 5 + +AHFinderDirect::origin_x[1] = 0 +AHFinderDirect::initial_guess__coord_sphere__x_center[1] = 0 +AHFinderDirect::initial_guess__coord_sphere__radius[1] = 1 +AHFinderDirect::which_surface_to_store_info[1] = 0 + +AHFinderDirect::surface_selection[2] = "areal radius" +AHFinderDirect::origin_x[2] = 0 +AHFinderDirect::initial_guess__coord_sphere__x_center[2] = 0 +AHFinderDirect::initial_guess__coord_sphere__radius[2] = 50 +AHFinderDirect::desired_value[2] = 50 +AHFinderDirect::which_surface_to_store_info[2] = 2 + +AHFinderDirect::surface_selection[3] = "areal radius" +AHFinderDirect::origin_x[3] = 0 +AHFinderDirect::initial_guess__coord_sphere__x_center[3] = 0 +AHFinderDirect::initial_guess__coord_sphere__radius[3] = 100 +AHFinderDirect::desired_value[3] = 100 +AHFinderDirect::which_surface_to_store_info[3] = 3 + +AHFinderDirect::surface_selection[4] = "areal radius" +AHFinderDirect::origin_x[4] = 0 +AHFinderDirect::initial_guess__coord_sphere__x_center[4] = 0 +AHFinderDirect::initial_guess__coord_sphere__radius[4] = 150 +AHFinderDirect::desired_value[4] = 150 +AHFinderDirect::which_surface_to_store_info[4] = 4 + +AHFinderDirect::surface_selection[5] = "areal radius" +AHFinderDirect::origin_x[5] = 0 +AHFinderDirect::initial_guess__coord_sphere__x_center[5] = 0 +AHFinderDirect::initial_guess__coord_sphere__radius[5] = 250 +AHFinderDirect::desired_value[5] = 250 +AHFinderDirect::which_surface_to_store_info[5] = 5 + +ActiveThorns = "Multipole WeylScal4" + +WeylScal4::verbose = 1 +Multipole::integration_method = "simpson" +Multipole::nradii = 4 +Multipole::radius[0] = 50 +Multipole::radius[1] = 100 +Multipole::radius[2] = 150 +Multipole::radius[3] = 250 +Multipole::variables = "WeylScal4::Psi4r{sw=-2 cmplx='WeylScal4::Psi4i'}" +Multipole::l_max = 4 +Multipole::interpolator_name = "Lagrange polynomial interpolation" + +ActiveThorns = "ML_Kretschmann ML_Kretschmann_Helper" +ML_Kretschmann::conformalMethod = 1 +ML_Kretschmann::ML_Kretschmann_kretschmann_calc_every = 128 + +ActiveThorns = "QuasiLocalMeasures" +QuasiLocalMeasures::num_surfaces = 5 +QuasiLocalMeasures::surface_index[0] = 0 +QuasiLocalMeasures::surface_index[1] = 2 +QuasiLocalMeasures::surface_index[2] = 3 +QuasiLocalMeasures::surface_index[3] = 4 +QuasiLocalMeasures::surface_index[4] = 5 + +############################################################################### +# OUTPUT # +############################################################################### + +ActiveThorns = "CarpetIOBasic CarpetIOHDF5" + +IO::out_dir = $parfile +IO::out_criterion = "time" +IO::out_dt = 0.25 +IO::out_mode = "onefile" +IO::out_unchunked = 1 +IOHDF5::checkpoint = 1 +IO::checkpoint_every = 4096 +IO::checkpoint_on_terminate = 1 +IO::checkpoint_dir = $parfile +IO::recover = "autoprobe" +IO::recover_dir = $parfile + +IOBasic::outInfo_criterion = "iteration" +IOBasic::outInfo_every = 128 +IOBasic::outInfo_reductions = "minimum maximum average norm1 norm2" +IOBasic::outInfo_vars = " + ML_BSSN::H + PunctureTracker::pt_loc_x[0] + PunctureTracker::pt_vel_x[0] + QuasiLocalMeasures::qlm_adm_energy[0] + QuasiLocalMeasures::qlm_adm_energy[1] + QuasiLocalMeasures::qlm_adm_energy[2] + QuasiLocalMeasures::qlm_adm_energy[3] + QuasiLocalMeasures::qlm_adm_energy[4] +" + +IOHDF5::compression_level = 5 +IOHDF5::output_symmetry_points = 0 +IOHDF5::out3D_ghosts = 0 +IOHDF5::out3D_outer_ghosts = 0 +IOHDF5::out_vars = " + ML_BSSN::ML_metric + ML_BSSN::ML_curv + ML_BSSN::ML_trace_curv + ML_BSSN::ML_log_confac + ML_BSSN::ML_lapse + ML_BSSN::ML_shift + ML_BSSN::ML_dtshift + ML_BSSN::ML_Ham + ML_BSSN::ML_mom + ML_BSSN::ML_cons_detg + ML_BSSN::ML_cons_traceA + ML_BSSN::ML_Gamma + AHFinderDirect::ah_radius + AHFinderDirect::ah_origin + ML_Kretschmann::Kretsch + PunctureTracker::pt_loc_x[0] + PunctureTracker::pt_vel_x[0] + SphericalSurface::sf_info + SphericalSurface::sf_radius + QuasiLocalMeasures::qlm_scalars + " -- cgit v1.2.3