I am trying to read the can.e.4.* Parallel Exodus II files in the VTKData tar ball using the vtkPExodusIIReader. It seems that the reader recognises which files needs to be read and returns without error. However, the reader does not recognise the blocks and arrays that should be part of that file. The returned vtkMultiBlockDataset appears empty, and note that I am able to use the vtkExodusIIReader on each of the 4 files. I am not sure what I have missed. I have attached my effort below, noting that it includes some pyvista
to streamline the visualisation.
import vtk
from pyvista import MultiBlock
def main():
controller = vtk.vtkDummyController()
vtk.vtkDummyController.SetGlobalController(controller)
reader = vtk.vtkPExodusIIReader()
reader.SetFilePattern(r"%s.e.4.%d")
reader.SetFilePrefix("can")
reader.SetFileRange(0,3)
reader.UpdateInformation()
print(reader.GetNumberOfElementResultArrays())
print(reader.GetNumberOfFiles())
print(reader.GetTotalNumberOfElements())
reader.SetPointResultArrayStatus("temp", 1)
tstep = reader.GetNumberOfTimeSteps() - 1
reader.SetTimeStep(tstep)
reader.Update()
print(reader)
data = MultiBlock(reader.GetOutputDataObject(0))
print(data)
print(data)
data.plot()
if __name__ == '__main__':
main()
Thanks in advance!
Matthew