# Probe integration point results to make surface plot of results on mesh

from vtkpython import *

# create readers for integration points and mesh
ipReader = vtkUnstructuredGridReader()
ipReader.SetFileName('brick_8.ipr')
ipReader.SetScalarsName('data')
meshReader = vtkUnstructuredGridReader()
meshReader.SetFileName('brick_8.vtk')

# create a probefilter
probe = vtkProbeFilter()
probe.SetInput(meshReader.GetOutput())
probe.SetSource(ipReader.GetOutput())

# create actor for mesh
meshMapper = vtkDataSetMapper() 
meshMapper.SetInput(probe.GetOutput())
meshActor = vtkActor()
meshActor.SetMapper(meshMapper)

# create a renderer and add actors to it
ren1 = vtkRenderer()
ren1.AddActor(meshActor)
# create render window and add renderer
renWin = vtkRenderWindow()
renWin.AddRenderer(ren1)

# create the interactor
iren = vtkRenderWindowInteractor() 
iren.SetRenderWindow(renWin)
iren.Initialize()
iren.Start()