<div>You should just create a parallel Vec to hold the 1D data. DAGetLocalInfo() tells you all the local sizes.</div><div><br></div><div>Thanks Matt. But I am not sure how this would help since PETSC_COMM_WORLD includes all the processors and in all three directions whereas the 1-D data is only the size of Nx points. </div>

<div>This vector is not big at all, I could even create is on one processor. The only concern that I have is </div><div>can I create an *.h5 file using PETSC_COMM_WORLD, dump the data, close the file and then, re-open the file (append mode) with a different viewer created via PETSC_COMM_SELF, and dump the coordinates vector to the end of it? Of course, I only call PETSC_COMM_SELF on processor zero. </div>

<div><br></div><div>Thanks, </div><div>M</div><div><br></div><br><div class="gmail_quote">On Thu, Jan 19, 2012 at 4:33 PM, Matthew Knepley <span dir="ltr">&lt;<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>&gt;</span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">On Thu, Jan 19, 2012 at 6:08 PM, Mohamad M. Nasr-Azadani <span dir="ltr">&lt;<a href="mailto:mmnasr@gmail.com" target="_blank">mmnasr@gmail.com</a>&gt;</span> wrote:<br>

</div><div class="gmail_quote"><div class="im"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
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></blockquote><div><br></div></div><div>You should just create a parallel Vec to hold the 1D data. DAGetLocalInfo() tells you all the local sizes.</div><div><br></div><div>   Matt</div>

<div class="im">
<div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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">&lt;<a href="mailto:jedbrown@mcs.anl.gov" target="_blank">jedbrown@mcs.anl.gov</a>&gt;</span> wrote:<br>




<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="gmail_quote">On Thu, Jan 19, 2012 at 18:00, Mohamad M. Nasr-Azadani <span dir="ltr">&lt;<a href="mailto:mmnasr@gmail.com" target="_blank">mmnasr@gmail.com</a>&gt;</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&#39;s why I don&#39;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>
</blockquote></div></div><br><br clear="all"><div class="HOEnZb"><div class="h5"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>


-- Norbert Wiener<br>
</div></div></blockquote></div><br>