[Users] Cell centered and Refluxing

[Yosef Zlochower]

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.
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: grid.par
Url: http://lists.einsteintoolkit.org/pipermail/users/attachments/20110304/da745f64/attachment.pl