Thanks, I was aware of such functions. I think I may have found where my confusion is. If I have a regular Vec, HDF5 SetTimestep and IncrementTimestep work as expected. Instead the following does not work for me:
[...] DMCreateGlobalVector(da, &solution); [...] // fill the vector solution and restore HDF5Open HDF5SetTimestep(H5_viewer, 0); ierr = VecView(solution,H5_viewer); ERRCHKQ(ierr); HDF5IncrementTimestep(H5_viewer); ierr = VecView(solution,H5_viewer); ERRCHKQ(ierr); [...] I get the following external library error from HDF at the line of the second VecView: HDF5-DIAG: Error detected in HDF5 (1.8.6) MPI-process 0: #000: H5D.c line 170 in H5Dcreate2(): unable to create dataset major: Dataset minor: Unable to initialize object #001: H5Dint.c line 431 in H5D_create_named(): unable to create and link to dataset major: Dataset minor: Unable to initialize object #002: H5L.c line 1640 in H5L_link_object(): unable to create new link to object major: Links minor: Unable to initialize object #003: H5L.c line 1866 in H5L_create_real(): can't insert link major: Symbol table minor: Unable to insert object #004: H5Gtraverse.c line 929 in H5G_traverse(): internal path traversal failed major: Symbol table minor: Object not found #005: H5Gtraverse.c line 718 in H5G_traverse_real(): traversal operator failed major: Symbol table minor: Callback failed #006: H5L.c line 1675 in H5L_link_cb(): name already exists major: Symbol table minor: Object already exists
