1#!/usr/bin/env python 2import vtk 3from vtk.util.misc import vtkGetDataRoot 4VTK_DATA_ROOT = vtkGetDataRoot() 5 6# Debugging parameters 7sze = 300 8 9reader = vtk.vtkExodusIIReader() 10reader.SetFileName(VTK_DATA_ROOT + "/Data/disk_out_ref.ex2") 11reader.UpdateInformation() 12reader.SetPointResultArrayStatus("CH4", 1) 13reader.Update() 14 15tree = vtk.vtkSpanSpace() 16 17contour = vtk.vtkSMPContourGrid() 18#contour = vtk.vtkContourGrid() 19contour.SetInputConnection(reader.GetOutputPort()) 20contour.SetInputArrayToProcess(0, 0, 0, vtk.vtkAssignAttribute.POINT_DATA, "CH4") 21contour.SetValue(0, 0.000718448) 22contour.MergePiecesOff() 23contour.UseScalarTreeOn() 24contour.SetScalarTree(tree) 25contour.Update() 26 27mapper = vtk.vtkCompositePolyDataMapper2() 28mapper.SetInputConnection(contour.GetOutputPort()) 29mapper.ScalarVisibilityOff() 30 31actor = vtk.vtkActor() 32actor.SetMapper(mapper) 33 34# Create graphics stuff 35# 36ren1 = vtk.vtkRenderer() 37ren1.AddActor(actor) 38ren1.ResetCamera() 39renWin = vtk.vtkRenderWindow() 40renWin.AddRenderer(ren1) 41renWin.SetSize(sze, sze) 42 43iren = vtk.vtkRenderWindowInteractor() 44iren.SetRenderWindow(renWin) 45 46renWin.Render() 47iren.Start() 48