ParaView/Automation with Python

1. Extract surface using paraview, the code will look like

from paraview.simple import *

#========== read in data file

print 'Read phasta file'

dataObj = PhastaReader( FileName= 'restart_one.pht' )

SRF = ExtractSurface(Input=dataObj)

MGBLK1 = MergeBlocks(Input=SRF)

vtkWriter = DataSetWriter(FileName = 'Surface.vtk' , Input = MGBLK1, FileType = 1)

vtkWriter.UpdatePipeline()