hdf5 functionality in io/petsc.py
rebased