XMDF reader error

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