[Paraview] XDMF Reader

[D Haley]

Dear Paraview List,

I'm trying to read some HDF files, which contain dense voxel data. In my
calculations, I write out the HDF and wrapper XDMF data to store
timestep and other metadata.

This all works fine, until I want to have 2 different voxel sizes at 2
different timesteps (different number of voxels in each dimension).

When loading up the data as a series of time-steps, everything works
until I reach the first set of data from the new size. When loading the
data, the values in the voxels seem to become invalid, and the bounding
box around the data starts moving in the direction of (1,1,1).

I can load each timestep individually correctly, it is just when I try
to load the entire timeseries that it fails.

Does anyone know a solution for this? I have raw voxel grids in the
program that generates the data, so I need to wrap them in something to
provide a timestep. I'd rather avoid having to interface with VTK, as
the code I am working on does not have a VTK dependency - I'm just
post-processing with paraview.

I'm opening the series via:
File->Open, then select the tree item data-..xmdf

each file is written as:

data-1.hdf.xmdf
data-2.hdf.xmdf
....

The XMDF files contain something like this

Before problem timestep:
==========
<?xml version="1.0" ?>
<!DOCTYPE Xdmf SYSTEM "Xdmf.dtd" []>
<Xdmf xmlns:xi="http://www.w3.org/2003/XInclude" Version="2.2">
<Domain>
<Grid Name="Data" GridType="Uniform">
<Topology TopologyType="3DCORECTMesh" Dimensions="39 39 39"/>
<Geometry GeometryType="ORIGIN_DXDYDZ">
<DataItem Name="Origin" Dimensions="3" NumberType="Float" Precision="4"
Format="XML">-1.26667 -1.26667 -1.26667</DataItem>
<DataItem Name="Spacing" Dimensions="3" NumberType="Float" Precision="4"
Format="XML">0.0666667 0.0666667 0.0666667</DataItem>
</Geometry>
<Attribute Name="second 0" AttributeType="Scalar" Center="Node">
<DataItem Format="HDF" NumberType="Float" Precision="4" Dimensions="39
39 39">data-59.hdf:/voxels</DataItem></Attribute>
</Grid>
</Domain>
</Xdmf>
==========

After problem timestep:
==========
<?xml version="1.0" ?>
<!DOCTYPE Xdmf SYSTEM "Xdmf.dtd" []>
<Xdmf xmlns:xi="http://www.w3.org/2003/XInclude" Version="2.2">
<Domain>
<Grid Name="Data" GridType="Uniform">
<Topology TopologyType="3DCORECTMesh" Dimensions="31 31 31"/>
<Geometry GeometryType="ORIGIN_DXDYDZ">
<DataItem Name="Origin" Dimensions="3" NumberType="Float" Precision="4"
Format="XML">-1 -1 -1</DataItem>
<DataItem Name="Spacing" Dimensions="3" NumberType="Float" Precision="4"
Format="XML">0.0666667 0.0666667 0.0666667</DataItem>
</Geometry>
<Attribute Name="second 0" AttributeType="Scalar" Center="Node">
<DataItem Format="HDF" NumberType="Float" Precision="4" Dimensions="31
31 31">data-60.hdf:/voxels</DataItem></Attribute>
</Grid>
</Domain>
</Xdmf>
==========


Thanks!