<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><br class=""><div> I looked at the code VecView_MPI_HDF5 and VecView_MPI_HDF5_DA they both seem to use the timestep variable to allow putting multiple datasets into the fire. Perhaps you forgot call PetscViewerHDF5SetTimestep for the DMDA case? Otherwise it looks like there must be a bug in VecView_MPI_HDF5_DA where it does not properly save multiple fields. You would need to step through that routine in the debugger.</div><div><br class=""></div><div><br class=""></div><div> Barry</div><div><br class=""><blockquote type="cite" class=""><div class="">On Jun 23, 2020, at 7:55 AM, Lucas Banting <<a href="mailto:bantingl@myumanitoba.ca" class="">bantingl@myumanitoba.ca</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">Hello,</div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class=""><br class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">I am trying to use the HDF5 viewer to store the solution of my system at every iteration.<span class="Apple-converted-space"> </span><br class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">I am storing the solutions in order to compare two different results, one from an old sequential program and one from a new parallel program that is parallelized using the DMDA structure.</div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">I am using matlab to compare the results.<br class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class=""><br class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">In the sequential program, the vector for the solution is just a simple one dimensional vector made from VecCreate. In the parallel program, the vector is made from DMCreateGlobalVector.</div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class=""><br class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class="">In both versions of the program, I create the viewer as:</div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class=""><br class=""></div><div style="font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; font-family: Calibri, Arial, Helvetica, sans-serif; font-size: 12pt;" class=""><div style="background-color: rgb(255, 255, 255); font-family: Consolas, "Courier New", monospace; font-weight: normal; font-size: 14px; line-height: 19px;" class=""><div class=""><span style="color: rgb(175, 0, 219);" class="">call</span><span class=""> </span><span style="color: rgb(121, 94, 38);" class="">PetscViewerHDF5Open</span><span class="">(PETSC_COMM_WORLD,</span><span style="color: rgb(163, 21, 21);" class="">"filename.hdf5"</span><span class="">,&</span></div><div class=""><span class=""> FILE_MODE_WRITE,hdf5viewer,ierr)</span></div><div class=""><span class=""><br class=""></span></div><div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; line-height: normal; background-color: rgba(0, 0, 0, 0);" class="">When I load the files into MATLAB after words, the sequential program creates a new dataset each iteration, while the parallel version using DMCreateGlobalVector only stores a single dataset.</span><span class=""><br class=""></span></div><div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; line-height: normal; background-color: rgba(0, 0, 0, 0);" class="">Is there a way I can force both viewers to create a new dataset in the same way each time vecview is called? Or is there someway I can manually create a new dataset?</span><span class=""><br class=""></span></div><div class=""><span class=""><br class=""></span></div><div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; line-height: normal; background-color: rgba(0, 0, 0, 0);" class="">I am using Open MPI 3.1.4, Petsc 3.12.4, Intel Fortran 2019.5, and hdf5/1.10.5.</span></div><div class=""><span class=""><br class=""></span></div><div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; line-height: normal; background-color: rgba(0, 0, 0, 0);" class="">Thanks,</span></div><div class=""><span class=""><br class=""></span></div><div class=""><span style="font-family: calibri, arial, helvetica, sans-serif; font-size: 12pt; line-height: normal; background-color: rgba(0, 0, 0, 0);" class="">Lucas Banting</span></div></div></div></div></blockquote></div><br class=""></body></html>