Hi Ken and Sven,

> Hmm, I'm not sure you have access to the element attribute or block
> properties within ParaView. ParaView is really optimized for reading
> results data.
> 
> Dave, do you know if there is a way to get at this information from an
> Exodus file?


Here is a small patch that will get the VTK reader loading data. A vtkpython 
script is also attached that demonstrates it on the exEdgeFaceElem.exii test 
file in vtkData.

However, I leave it as an exercise to the reader to expose things in 
ParaView... it was hard coming up with all the little icons for different types 
of data arrays and now you need 3 more (edge, face, and element block 
attributes)! Plus you may need to add some convenience methods to the 
vtkExodusIIReader class to expose the different types of attributes in a way 
that lets ParaView map them to properties.

        David

Attachment: 0001-Read-requested-Exodus-block-attributes.patch
Description: Binary data

#!/usr/bin/env vtkpython
from vtk import *
rdr = vtkExodusIIReader()
rdr.SetFileName('/path/to/VTKData/Data/edgeFaceElem.exii')
rdr.UpdateInformation()
rdr.SetObjectAttributeStatus(rdr.ELEM_BLOCK,0,'SPAGHETTI',1)
rdr.SetObjectAttributeStatus(rdr.ELEM_BLOCK,0,'WESTERN',1)
rdr.Update()
na = rdr.GetOutputDataObject(0).GetBlock(0).GetBlock(0).GetCellData().GetNumberOfArrays()
for i in range(na):
  print rdr.GetOutputDataObject(0).GetBlock(0).GetBlock(0).GetCellData().GetArray(i).GetName()
_______________________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Follow this link to subscribe/unsubscribe:
http://www.paraview.org/mailman/listinfo/paraview

Reply via email to