<div dir="ltr"><div><p dir="ltr">Hi Luca,</p>Try this -> open up both data sets, select both in the Pipeline Browser and apply the Append Attributes filter.<br><br></div>hope that help<br><div><br><br>On Mar 30, 2016 5:11 AM, "Luca Pavarino" <<a href="mailto:luca.pavarino@unimi.it" target="_blank">luca.pavarino@unimi.it</a>> wrote:<br type="attribution"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
I need to visualize in Paraview at every time step a scalar field (in my PETSc code Vec cc_1 below) over a moving structured mesh (Vec um). Both the field and the mesh are solutions of some a parallel PETSc time-dependent PDE code. Do you have a simple PETSc example to output these vectors at each time step and see them in Paraview (either each snapshot separately or with a movie)?<br>
<br>
I tried with the simple PETSc code (Vec cc_1 is the scalar field and Vec uum_1 is the vector of 3D coordinates):<br>
<br>
 write (filename,"(A, I4.4, A)") "uuvtk", int(k*dt),".vts"<br>
 call PetscViewerVTKOpen(PETSC_COMM_WORLD,filename,<br>
 & FILE_MODE_WRITE,view_out_vtk,ierr);<br>
 call VecView(cc_1,view_out_vtk,ierr)<br>
 call PetscViewerDestroy(view_out_vtk,ierr)<br>
<br>
 write (filename,"(A, I4.4, A)") "uumvtk", int(k*dt),".vts"<br>
 call PetscViewerVTKOpen(PETSC_COMM_WORLD,filename,<br>
 & FILE_MODE_WRITE,view_out_vtk,ierr);<br>
 call VecView(uum_1,view_out_vtk,ierr)<br>
 call PetscViewerDestroy(view_out_vtk,ierr)<br>
<br>
but I get 2 separate vectors and I do not know how to make Paraview see that the scalar field is defined over the coordinate mesh. Thanks for any help,<br>
<br>
 Luca<br>
_______________________________________________<br>
Powered by <a href="http://www.kitware.com" rel="noreferrer" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" rel="noreferrer" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the ParaView Wiki at: <a href="http://paraview.org/Wiki/ParaView" rel="noreferrer" target="_blank">http://paraview.org/Wiki/ParaView</a><br>
<br>
Search the list archives at: <a href="http://markmail.org/search/?q=ParaView" rel="noreferrer" target="_blank">http://markmail.org/search/?q=ParaView</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/paraview" rel="noreferrer" target="_blank">http://public.kitware.com/mailman/listinfo/paraview</a><br>
</blockquote></div>
</div></div>