sorry for the messages before, hopefully the mail is now small enough, but
i thought that adding a plot might explain better what i mean..
Hello all,
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