Thanks Jed. <div>My solution to that was to create a 1D vector local to for instance processor zero that holds all the coordinates. Then after I dumped all the parallel data to the *.h5 file, I create another viewer on processor zero or PETSC_COMM_SELF and dump that new vector including tHe coordinates to the end of the existing file. </div>
<div>Do you think that should be possible?</div><div><br></div><div>Thanks, </div><div>Mohamad</div><div><br></div><div><br></div><div><br><br><div class="gmail_quote">On Thu, Jan 19, 2012 at 4:02 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im"><div class="gmail_quote">On Thu, Jan 19, 2012 at 18:00, Mohamad M. Nasr-Azadani <span dir="ltr"><<a href="mailto:mmnasr@gmail.com" target="_blank">mmnasr@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
What I need is just to add 3 1-D arrays of (x[Nx]+y[Ny]+z[Nz]) including the grid coordinates to the end of the *.h5 file and then later on, I can use any visualization software to load the data using those coordinates. I am using orthogonal grid, that's why I don't need all the (x,y,z) coordinates for each cell. </blockquote>
</div><br></div><div>If you want this special case, you have to manage it by hand.</div>
</blockquote></div><br></div>