[Users] Cell centered and Refluxing
Erik Schnetter
schnetter at cct.lsu.edu
Mon Mar 7 15:41:20 CST 2011
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 at 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 at 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 at einsteintoolkit.org
> http://lists.einsteintoolkit.org/mailman/listinfo/users
>
>
--
Erik Schnetter <schnetter at cct.lsu.edu> http://www.cct.lsu.edu/~eschnett/
More information about the Users
mailing list