[Users] Cell centered and Refluxing

-- Dr. Yosef Zlochower Center for Computational Relativity and Gravitation Assistant Professor School of Mathematical Sciences Rochester Institute of Technology 85 Lomb Memorial Drive Rochester, NY 14623 Office:74-2067 Phone: +1 585-475-6103 yosef@astro.rit.edu CONFIDENTIALITY NOTE: The information transmitted, including attachments, is intended only for the person(s) or entity to which it is addressed and may contain confidential and/or privileged material. Any review, retransmission, dissemination or other use of, or taking of any action in reliance upon this information by persons or entities other than the intended recipient is prohibited. If you received this in error, please contact the sender and destroy any copies of this information. #============================================================================== ActiveThorns = "admanalysis admbase admcoupling admmacros coordgauge spacemask StaticConformal boundary time cartgrid3d ioutil iobasic nanchecker CoordBase aeilocalinterp Slab SphericalSurface LocalReduce MoL Carpet CarpetInterp CarpetIOASCII CarpetLib CarpetReduce CarpetSlab CarpetRegrid2 AHFinderDirect CarpetIOHDF5 CarpetIOScalar InitBase SymBase LoopControl GSL ReflectionSymmetry TmunuBase HydroBase" ActiveThorns = "TimerReport" TimerReport::out_every = 512 TimerReport::out_filename = "TimerReport" TimerReport::n_top_timers = 20 TimerReport::output_all_timers = yes Carpet::output_timers_every = 512 CarpetLib::print_memstats_every = 1024 #------------------------------------------------------------------------------ # Run parameters #------------------------------------------------------------------------------ # Grid mol::ode_method = "RK4" mol::mol_intermediate_steps = 4 mol::mol_num_scratch_levels = 1 #mol::initial_data_is_crap = true CoordBase::domainsize = "minmax" CoordBase::xmin = -200.001 CoordBase::ymin = -200.001 CoordBase::zmin = 0 CoordBase::xmax = 200 CoordBase::ymax = 200 CoordBase::zmax = 200 CoordBase::spacing = "numcells" CoordBase::ncells_x = 200 CoordBase::ncells_y = 200 CoordBase::ncells_z = 100 CartGrid3D::type = "coordbase" CartGrid3D::domain = "full" CartGrid3D::bitant_plane = "xy" CartGrid3D::avoid_originx = no CartGrid3D::avoid_originy = no CartGrid3D::avoid_originz = no CoordBase::boundary_size_x_lower = 1 CoordBase::boundary_size_y_lower = 1 CoordBase::boundary_size_z_lower = 4 CoordBase::boundary_shiftout_x_lower = 0 CoordBase::boundary_shiftout_y_lower = 0 CoordBase::boundary_shiftout_z_lower = 1 CoordBase::boundary_size_x_upper = 1 CoordBase::boundary_size_y_upper = 1 CoordBase::boundary_size_z_upper = 1 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" ReflectionSymmetry::avoid_origin_upper_x="no" ReflectionSymmetry::avoid_origin_upper_y="no" ReflectionSymmetry::avoid_origin_upper_z="no" driver::ghost_size = 4 #--------------------------------------------------------------------------- cactus::cctk_itlast = 12800000 time::dtfac = .1 Carpet::domain_from_coordbase = yes Carpet::prolongation_order_space = 5 Carpet::prolongation_order_time = 2 Carpet::max_refinement_levels = 8 Carpet::use_buffer_zones = yes Carpet::num_integrator_substeps = 1 Carpet::additional_buffer_zones = 2 Carpet::verbose = no Carpet::schedule_barriers = no Carpet::veryverbose = no CarpetLib::output_bboxes = no CarpetLib::check_bboxes = no CarpetLib::interleave_communications = yes CarpetLib::combine_sends = yes Carpet::time_refinement_factors = "[1,1,1,2,4,8,16,32,64,128,256]" Carpet::init_3_timelevels = "no" Carpet::init_each_timelevel = "no" Carpet::init_fill_timelevels = "yes" Carpet::enable_all_storage = no Carpet::regrid_during_recovery = no CarpetRegrid2::regrid_every = 32 CarpetRegrid2::snap_to_coarse="yes" CarpetRegrid2::num_centres = 1 CarpetRegrid2::num_levels_1 = 7 CarpetRegrid2::position_x_1 = 0 CarpetRegrid2::position_y_1 = 0 CarpetRegrid2::position_z_1 = 0 CarpetRegrid2::radius_1[ 1] = 90 CarpetRegrid2::radius_1[ 2] = 45 CarpetRegrid2::radius_1[ 3] = 20 CarpetRegrid2::radius_1[ 4] = 8 CarpetRegrid2::radius_1[ 5] = 4 CarpetRegrid2::radius_1[ 6] = 2 # MODEL initbase::initial_data_setup_method = init_some_levels ActiveThorns = "IDAnalyticBH" ADMBase::initial_data="schwarzschild" ADMBase::initial_lapse="psi4_psi2" ADMBase::initial_shift="zero" IDAnalyticBH::mass = 1. IDAnalyticBH::a_Kerr = 0. admbase::metric_type = "Physical" #--- numerics #------------------------------------------------------------------------------ # Matter: #------------------------------------------------------------------------------ ActiveThorns = "Grhydro EOS_Omni Constants" #ActiveThorns = "Refluxing" grhydro::riemann_solver = "Marquina" grhydro::grhydro_eos_type = "general" grhydro::grhydro_eos_table="Ideal_Fluid" grhydro::recon_method = "ppm" grhydro::grhydro_stencil = 3 grhydro::bound = "flat" HydroBase::Prolongation_Type = "WENO" Grhydro::initial_Gamma=1.66666666666666 Grhydro::rho_abs_min=1.0e-13 Grhydro::initial_rho_abs_min=1.0e-13 Grhydro::grhydro_atmo_tolerance=0.0001 eos_omni::poly_gamma = 1.66666666666666 eos_omni::poly_gamma_ini = 1.66666666666666 eos_omni::poly_k = 0.01 EOS_Omni::gl_gamma = 1.66666666666666 EOS_Omni::gl_k = 0.01 grhydro::Grhydro_MaxNumSandRVars = 16 TmunuBase::stress_energy_storage = "yes" TmunuBase::stress_energy_at_RHS = "yes" TmunuBase::timelevels = 1 TmunuBase::prolongation_type = "none" HydroBase::timelevels = 3 HydroBase::evolution_method = "grhydro" SpaceMask::use_mask = "yes"