i recently read the Reisswig et al. 2012 paper and started trying to switch to cell-centered grid in order to be able to use the LSUThorns/Refluxing thorn. i am trying to simulate BH+torus systems and due to the size of the torus its hard to have it contained within a refinement level, so i thought that using a cell-centered grid in order to use refluxing might improve results..

i however run into problems getting the hdf5 output of 3D grid frundtions to work...

using the static_tov_cc_rf.par file from the Refluxing trunk, i get weird output in the mesh at the reflection boundaries (the par file uses octant mode)..

i am using the following parameters for CarpetIOHDF5:

IOHDF5::out_every = 256

IOHDF5::one_file_per_group = no

IOHDF5::output_symmetry_points = no

IOHDF5::out3D_ghosts = no

IOHDF5::output_ghost_points = no

IOHDF5::out3D_outer_ghosts = no

IOHDF5::output_boundary_points = no

IOHDF5::output_buffer_points = no

IOHDF5::compression_level = 9

IOHDF5::use_checksums = yes

IOHDF5::out_vars = "

HydroBase::rho

HydroBase::press

HydroBase::eps

HydroBase::vel

ADMBase::lapse

ADMBase::metric

ADMBase::curv

GRHydro::scon

the simulations seems to work fine, so is the carpetiohdf5 thorn not suitable for cell-centered output? another effect (that might be output related) is that the different refinement levels are shifted (i'll attach a file showing this in the TOV from the par file (i added extra refinement levels to illustrate the effect) as well as a plot of the density in the torus)...

are the specific parameters for cell-centered hdf5 output or other thorns that could take care of that?

best wishes,

Vassili