Hi, all, I have implemented&nbsp; an MPI code and I want to use VTK's parallel file format(pvtu). The computational domain is already partition using METIS. What&nbsp; I want to do is: for i=1:1:N A) Each MPI has its own data to file   &nbsp;&nbsp; (data_i_0.vtu, data_i_1.vtu, ... data_i_MPI_Size-1.vtu). B) Then combine all the files to create data_i.pvtu file to visualoze the  &nbsp;&nbsp;&nbsp; whole domain.  end So far&nbsp; I have, vtkUnstructuredGrid* mesh = vtkUnstructuredGrid::New();  &nbsp; addNodes(mesh, nodeNum, nodeArray); &nbsp; addCells(mesh, tetNum, tetArray); &nbsp; addCellMaterial(mesh, tetNum, tetArray); &nbsp; addPointScalars(mesh, nodeNum, eField, "E"); &nbsp; addPointScalars(mesh, nodeNum, hField, "H");  &nbsp;&nbsp; // Parallel  &nbsp; string outputFile(outFilePrefix); &nbsp; outputFile.append(".pvtu"); &nbsp; vtkXMLPUnstructuredGridWriter *pwriter = vtkXMLPUnstructuredGridWriter::New(); &nbsp;&nbsp;&nbsp; pwriter-&gt;SetInput( mesh );  &nbsp;&nbsp;&nbsp; pwriter-&gt;SetFileName(outputFile.c_str()); //&nbsp;&nbsp; pwriter-&gt;SetDataModeToBinary(); &nbsp; pwriter-&gt;SetDataModeToAscii(); &nbsp; if (MPI_RankID == 0) &nbsp; { &nbsp;&nbsp;&nbsp; &nbsp; pwriter-&gt;SetNumberOfPieces( MPI_Size ); &nbsp;&nbsp;&nbsp; pwriter-&gt;Write();  &nbsp; } &nbsp; else &nbsp; { &nbsp;&nbsp;&nbsp; pwriter-&gt;SetNumberOfPieces(1); &nbsp;&nbsp;&nbsp; pwriter-&gt;SetStartPiece(MPI_RankID); &nbsp;&nbsp;&nbsp; pwriter-&gt;SetEndPiece(MPI_RankID); &nbsp;&nbsp;&nbsp; pwriter-&gt;Write(); &nbsp; } &nbsp; mesh-&gt;Delete(); &nbsp; pwriter-&gt;Delete();   I get the files&nbsp; that I want but when tryig to visualize the&nbsp; pvtu ifiles in VisIt  I get an error If anyone has an idea on how to do this&nbsp; I would appreciate it.