[kaya]

I have my own cfd solver in which I was using *.vti for paraview input but decided to do MPI and wanted to use *.pvti format and do node parallelisation.

My *.pvti file is like :

<?xml version="1.0"?>
<VTKFile type="PImageData" version="1.0" byte_order="LittleEndian" header_type="UInt64">
  <PImageData WholeExtent=[/QUOTE][/QUOTE]"       0      31       0      31       0      31" Origin="  0.000000000000E+00  0.000000000000E+00  0.000000000000E+00" Spacing="  0.196349540849E+00  0.196349540849E+00  0.645161290323E-01" GhostLevel="0">
      <PPointData Scalars="u v w">
    <PDataArray type="Float64" Name="u" NumberOfComponents="1" format="appended"/>
<PDataArray type="Float64" Name="v" NumberOfComponents="1" format="appended"/>
<PDataArray type="Float64" Name="w" NumberOfComponents="1" format="appended"/>
</PPointData>
<Piece Extent="0 31 0 31 0 15" Source="data0_0.vti"/>
<Piece Extent="0 31 0 31 16 31" Source="data0_1.vti"/>
</PImageData></VTKFile>

The problem I am facing is I am having a gap in between z-axis 15-16 - which paraview complains.

I know that the solution is to have a ghost (overlapping) layer (can't understand why since I am using point data?), and the issue is that my code is not capable of generating overlapping nodes, unless I code it up and manually transfer ghost layers at each data output stage for the sake of being able to write data in parallel but it feels like an overkill...

Any suggestions?

[kaya]

In case anyone get into the same problem, the solution is to write out cell data instead of point and in the pvti file put the borders of the cells whih you will be duplicating.