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

_______________________________________________
Users mailing list
Users@einsteintoolkit.org
http://lists.einsteintoolkit.org/mailman/listinfo/users