Hi all,
I am trying to define a vector field on a grid using the data from a
pair of hdf5 arrays, the hdf5 data layout is:
phil@work $ h5ls data.h5
map_x Dataset {11, 11}
map_y Dataset {11, 11}
registered Dataset {100, 100}
and I have produced an xdmf file with the following:
<?xml version="1.0" ?>
<!DOCTYPE Xdmf SYSTEM "Xdmf.dtd">
<Xdmf Version="3.0">
<Domain>
<Grid Name="MapGrid" GridType="Uniform">
<Topology TopologyType="2DCoRectMesh" Dimensions="11 11"/>
<Geometry GeometryType="Origin_DxDy">
<DataItem Name="Map_Origin" Dimensions="2" NumberType="Float" Format="XML">
0 0
</DataItem>
<DataItem Name="Map_Spacing" Dimensions="2" NumberType="Float" Format="XML">
10 10
</DataItem>
</Geometry>
<Attribute Name="Displacements" AttributeType="Vector" Center="Node">
<DataItem ItemType="Function" Function="join($0, $1)" Dimensions="11 11 2">
<DataItem Format="HDF" Dimensions="11 11">
data.h5:/map_x
</DataItem>
<DataItem Format="HDF" Dimensions="11 11">
data.h5:/map_y
</DataItem>
</DataItem>
</Attribute>
</Grid>
</Domain>
</Xdmf>
Trying to plot in paraview using the xdmf2 reader this gives an hdf5 error:
HDF5-DIAG: Error detected in HDF5 (1.10.1) MPI-process 0:
#000:
/var/tmp/portage/sci-libs/hdf5-1.10.1/work/hdf5-1.10.1/src/H5Dio.c line
256 in H5Dwrite(): memory selection+offset not within extent
major: Dataspace
minor: Out of range
XDMF Error in
/var/tmp/portage/sci-visualization/paraview-5.5.2/work/ParaView-v5.5.2/VTK/ThirdParty/xdmf2/vtkxdmf2/libsrc/XdmfHDF.cxx
line 897 (Can't Write Temp Dataset)
Is there something wrong with my xdmf?
Thanks
Phil
–
Phil Tooley
Research Software Engineering
University of Sheffield