Yosef
The thorn LSUDevelopment/Refluxing contains a sample static_tov.par in three versions:
1. static_tov.par the original, vertex centered case 2. static_tov_cc.par with vertex centering 3. static_tov_cc_rf.par with vertex centering and refluxing
If you diff static_tov.par and static_tov_cc.par, you should see what changes are necessary. These are probably related to setting up the coordinates and the extent of the simulation domain, as well as choosing the interpolation orders correspondingly. Note that there are no cell-centred ENO interpolators yet, so you'll have to use linear interpolation for hydrodynamics.
-erik
On Fri, Mar 4, 2011 at 11:36 AM, Yosef Zlochower yosef@astro.rit.edu wrote:
Hi,
I would like to try using cell-centered refinement and Refluxing. I have a parfile (attached) with the grid setup for vertex-centered refinement. What parameters do I need the change to make the parfile work for cell-centered? Based on the error messages when trying cell-centered I made following changes to the parfile
CoordBase::domainsize = "minmax" CoordBase::xmin = -200.00 CoordBase::ymin = -200.00 CoordBase::zmin = 0 CoordBase::xmax = 200 CoordBase::ymax = 200 CoordBase::zmax = 200 CoordBase::spacing = "numcells" CoordBase::ncells_x = 201 CoordBase::ncells_y = 201 CoordBase::ncells_z = 101 CartGrid3D::type = "coordbase" CartGrid3D::domain = "full" CartGrid3D::bitant_plane = "xy" CartGrid3D::avoid_originx = yes CartGrid3D::avoid_originy = yes CartGrid3D::avoid_originz = yes
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 = 0
CoordBase::boundary_size_x_upper = 1 CoordBase::boundary_size_y_upper = 1 CoordBase::boundary_size_z_upper = 1
CoordBase::boundary_staggered_x_lower = "yes" CoordBase::boundary_staggered_x_upper = "yes" CoordBase::boundary_staggered_y_lower = "yes" CoordBase::boundary_staggered_y_upper = "yes" CoordBase::boundary_staggered_z_lower = "yes" CoordBase::boundary_staggered_z_upper = "yes"
ReflectionSymmetry::reflection_x="no" ReflectionSymmetry::reflection_y="no" ReflectionSymmetry::reflection_z="yes" ReflectionSymmetry::avoid_origin_x="yes" ReflectionSymmetry::avoid_origin_y="yes" ReflectionSymmetry::avoid_origin_z="yes" ReflectionSymmetry::avoid_origin_upper_x="yes" ReflectionSymmetry::avoid_origin_upper_y="yes" ReflectionSymmetry::avoid_origin_upper_z="yes"
Carpet::refinement_centering = "cell"
but I then get the error message
(line 136 of prolongate_3d_real8_eno.F90): -> Internal error: array origins are not integer multiples of the strides WARNING level 0 in thorn CarpetLib processor 12 host i182-302.ranger.tacc.utexas.edu
If I further change the HydroBase prolongation operator to Lagrange I get errors of the form
WARNING[L1,P2] (AEILocalInterp): CCTK_InterpLocalUniform(): interpolation point is either outside the grid, or inside but too close to the grid boundary! (this may be caused by a global interpolation with driver::ghost_size too small) 0-origin interpolation point number pt=0 of N_interp_points=81 interpolation point (x,y,z)=(6.09375e-17,0.995185,0.0980171) grid x_min(delta_x)x_max = -200.995(1.99005)-93.5323 grid y_min(delta_y)y_max = -107.463(1.99005)7.9602 grid z_min(delta_z)z_max = -6.93069(1.9802)103.96 which is due to an analysis interpolation (an interpolation over spheres centered on the origin in global mode)
I tried a parfile with no symmetries with similar results.
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"
Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users