Cactus::cctk_run_title = "TOV collapse" Cactus::cctk_full_warnings = yes Cactus::terminate = time # never # Cactus::cctk_final_time = 10.0 ActiveThorns = "Fortran" ActiveThorns = "LocalInterp AEILocalInterp" ActiveThorns = "Slab" ActiveThorns = "IOUtil" ActiveThorns = "InitBase Carpet CarpetLib CarpetInterp CarpetReduce LoopControl" ActiveThorns = "NaNChecker" ActiveThorns = "Boundary CartGrid3D CoordBase SymBase " ActiveThorns = "CarpetRegrid2 SphericalSurface" ActiveThorns = "ADMBase ADMCoupling ADMMacros CoordGauge SpaceMask StaticConformal" ActiveThorns = "Time" ActiveThorns = "MoL" ActiveThorns = "WeylScal4 Multipole GenericFD" ActiveThorns = "ML_BSSN ML_BSSN_Helper NewRad ML_ADMConstraints" ActiveThorns = "ADMAnalysis" ActiveThorns = "AHFinderDirect" ActiveThorns = "CarpetIOScalar CarpetIOASCII CarpetIOBasic" # CarpetIOHDF5" ActiveThorns = "TimerReport Formaline" ActiveThorns = "ReflectionSymmetry" # RotatingSymmetry180 ActiveThorns = "HydroBase TmunuBase InitBase" # IsolatedHorizon ActiveThorns = "eos_base eos_polytrope eos_idealfluid" ActiveThorns = "whisky whisky_tovsolverc constants" #Whisky_MeudonData ## Initial data ADMBase::metric_type = physical ADMBase::initial_data = tov ADMBase::initial_lapse = tov ADMBase::initial_shift = tov ADMBase::initial_dtlapse = "zero" ADMBase::initial_dtshift = "zero" whiskytovsolver::TOV_Rho_Central[0] = 8.e-3 whiskytovsolver::TOV_Gamma[0] = 2.0 whiskytovsolver::TOV_K[0] = 110.0 whiskytovsolver::TOV_Num_Radial = 1000000 whiskytovsolver::TOV_dr[0] = 5.e-4 whiskytovsolver::TOV_dr[1] = 5.e-4 InitBase::initial_data_setup_method = "init_all_levels" MoL::initial_data_is_crap = yes ### Grid setup Time::dtfac = 0.35 Carpet::verbose = no Carpet::domain_from_coordbase = yes CartGrid3D::type = coordbase ReflectionSymmetry::reflection_x = yes ReflectionSymmetry::reflection_y = yes ReflectionSymmetry::reflection_z = yes ReflectionSymmetry::avoid_origin_x = no ReflectionSymmetry::avoid_origin_y = no ReflectionSymmetry::avoid_origin_z = no #CarpetRegrid2::symmetry_rotating180 = yes CoordBase::boundary_size_x_lower = 3 CoordBase::boundary_shiftout_x_lower = 1 CoordBase::boundary_size_y_lower = 3 CoordBase::boundary_shiftout_y_lower = 1 CoordBase::boundary_size_z_lower = 3 CoordBase::boundary_shiftout_z_lower = 1 CoordBase::boundary_size_x_upper = 3 CoordBase::boundary_size_y_upper = 3 CoordBase::boundary_size_z_upper = 3 CoordBase::domainsize = minmax CoordBase::xmax = 25.0 CoordBase::ymax = 25.0 CoordBase::zmax = 25.0 CoordBase::xmin = 0.0 #-25.0 # CoordBase::ymin = 0.0 #-25.0 # CoordBase::zmin = 0.0 #-25.0 # CoordBase::dx = 0.5 CoordBase::dy = 0.5 CoordBase::dz = 0.5 Carpet::max_refinement_levels = 1 CarpetRegrid2::num_centres = 1 CarpetRegrid2::num_levels_1 = 1 CarpetRegrid2::position_x_1 = 0.0 CarpetRegrid2::radius_1 [1] = 16.0 Driver::ghost_size = 3 Carpet::use_buffer_zones = yes Carpet::prolongation_order_space = 5 Carpet::prolongation_order_time = 2 Carpet::convergence_level = 0 Carpet::regrid_in_level_mode = yes Carpet::regrid_during_initialisation = no Carpet::init_each_timelevel = no ### Evolution mol::ode_method = rk4 MoL::MoL_Intermediate_Steps = 4 MoL::MoL_Num_Scratch_Levels = 1 Carpet::num_integrator_substeps = 4 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::advectLapse = 1.0 ML_BSSN::advectShift = 1.0 ML_BSSN::MinimumLapse = 1.0e-8 Boundary::radpower = 3 #2 ML_ADMConstraints::verbose = 10 ### Horizon and extraction surfaces # We choose the angular resolution such that it is approx. of same # magnitude as the underlying Cartesian grid. # Note that the angular resolution scales with the radius of the sphere. # -> spheres with smaller R need less angular gridpoints. # (Here I have chosen the same res. as the sphere with largest R) # Also note that IsolatedHorizon is computationally intensive and gets # low res surfaces SphericalSurface::nsurfaces = 2 #5 SphericalSurface::maxntheta = 140 #156 #73 SphericalSurface::maxnphi = 240 #156 #76 SphericalSurface::ntheta [0] = 24 SphericalSurface::nphi [0] = 96 SphericalSurface::symmetric_z [0] = yes SphericalSurface::nghoststheta [0] = 2 SphericalSurface::nghostsphi [0] = 2 SphericalSurface::ntheta [1] = 24 SphericalSurface::nphi [1] = 96 SphericalSurface::symmetric_z [1] = yes SphericalSurface::nghoststheta [1] = 2 SphericalSurface::nghostsphi [1] = 2 #### Apaprent Horizons AHFinderDirect::N_horizons = 2 AHFinderDirect::find_every = 32 AHFinderDirect::output_h_every = 320 AHFinderDirect::max_Newton_iterations__initial = 50 AHFinderDirect::max_Newton_iterations__subsequent = 50 AHFinderDirect::max_allowable_Theta_growth_iterations = 10 AHFinderDirect::max_allowable_Theta_nonshrink_iterations = 10 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::verbose_level = "physics details" AHFinderDirect::move_origins = yes AHFinderDirect::origin_x [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::initial_guess__coord_sphere__radius [1] = 0.25 AHFinderDirect::which_surface_to_store_info [1] = 0 AHFinderDirect::set_mask_for_individual_horizon [1] = no AHFinderDirect::reset_horizon_after_not_finding [1] = no AHFinderDirect::find_after_individual_time [1] = 3000 AHFinderDirect::max_allowable_horizon_radius [1] = 15.0 AHFinderDirect::origin_x [2] = 0.0 AHFinderDirect::initial_guess__coord_sphere__x_center [2] = 0.0 AHFinderDirect::initial_guess__coord_sphere__radius [2] = 1.50 AHFinderDirect::which_surface_to_store_info [2] = 1 AHFinderDirect::set_mask_for_individual_horizon [2] = no AHFinderDirect::reset_horizon_after_not_finding [2] = no AHFinderDirect::find_after_individual_time [2] = 3000 AHFinderDirect::max_allowable_horizon_radius [2] = 15.0 AHFinderDirect::reshape_while_moving = yes AHFinderDirect::predict_origin_movement = yes ### NaNCheck NaNChecker::check_every = 6400 NanChecker::check_after = 0 NaNChecker::action_if_found = terminate NaNChecker::check_vars = "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" #all # NaNChecker::out_NaNmask = yes ### Output IO::out_dir = "TOV_k100" IOBasic::outInfo_every = 1 IOBasic::outInfo_reductions = "norm2" IOBasic::outInfo_vars = "ML_ADMConstraints::ML_Ham" IOScalar::one_file_per_group = no IOScalar::outScalar_reductions = "minimum maximum norm1 norm2" IOScalar::outScalar_criterion = divisor IOScalar::outScalar_every = 2 #64 IOScalar::outScalar_vars = " hydrobase::rho{ reductions = 'minimum maximum'} hydrobase::eps{ reductions = 'minimum maximum'} admbase::lapse{ reductions = 'minimum'} ML_ADMConstraints::ML_Ham{ reductions = 'norm1 norm2 norm_inf'} ML_ADMConstraints::ML_mom{ reductions = 'norm1 norm2 norm_inf'} " IOASCII::out0D_criterion = divisor IOASCII::out0D_every = 2 #64 IOASCII::out0D_vars = " Carpet::timing " CarpetIOASCII::one_file_per_group = no IOASCII::out1D_criterion = divisor IOASCII::out1D_d = no IOASCII::out1D_every = 10 IOASCII::out1D_vars = " ML_ADMConstraints::ML_Ham ML_ADMConstraints::ML_mom hydrobase::rho hydrobase::press whisky::dens hydrobase::eps hydrobase::vel ADMBase::metric " CarpetLib::print_timestats_every = 12800 CarpetLib::print_memstats_every = 12800 #====================================== # Whisky with PPM parameters #====================================== TmunuBase::stress_energy_storage = yes TmunuBase::stress_energy_at_RHS = yes TmunuBase::timelevels = 1 TmunuBase::prolongation_type = none TmunuBase::support_old_CalcTmunu_mechanism = yes hydrobase::timelevels = 3 hydrobase::evolution_method = "whisky" hydrobase::prolongation_type = ENO #default whisky::riemann_solver = "Marquina" whisky::whisky_eos_type = "Polytype" #General # whisky::whisky_eos_table = "2D_Polytrope" #Ideal_Fluid # whisky::recon_method = "ppm" whisky::ppm_detect = "yes" whisky::ppm_omega2 = 0.0 whisky::ppm_flatten = "stencil_3" whisky::whisky_stencil = 3 whisky::bound = "none" spacemask::use_mask = "true" whisky::whisky_c2p_warnlevel = 0 whisky::whisky_c2p_warn_from_reflevel = 1 whisky::c2p_reset_pressure = yes #====================================== # Whisky atmosphere parameters #====================================== #whisky::rho_abs_min = 1.e-13 whisky::rho_rel_min = 1.e-7 whisky::initial_atmosphere_factor = 0.1 whisky::whisky_atmo_tolerance = 0.0001 eos_polytrope::eos_gamma = 2.0 eos_polytrope::eos_k = 100.0 #123.6