<div dir="ltr">David,<div><br></div><div>Thanks for the motivation to keep digging, I finally figured it out. My problem was how I was interpreting the indices added to the vtkExtractBlock.</div><div><br></div><div>Just in case anyone else is interested, the following prints out the indices of an Exodus file to use with vtkExtractBlock.</div><div><br></div><div><div># Read Exodus Data</div><div>reader = vtk.vtkExodusIIReader()</div><div>reader.SetFileName('mug_blocks_out.e')</div><div>reader.SetTimeStep(10)</div><div>reader.UpdateInformation()</div><div>reader.Update()</div><div><br></div><div>reader.SetAllArrayStatus(vtk.vtkExodusIIReader.NODAL, 1)</div><div>reader.SetAllArrayStatus(vtk.vtkExodusIIReader.NODAL_TEMPORAL, 1)</div><div><br></div><div>reader.SetObjectStatus(vtk.vtkExodusIIReader.ELEM_BLOCK, 0, 1)</div><div>reader.SetObjectStatus(vtk.vtkExodusIIReader.ELEM_BLOCK, 1, 1)</div><div>reader.SetObjectStatus(vtk.vtkExodusIIReader.SIDE_SET, 0, 1)</div><div>reader.SetObjectStatus(vtk.vtkExodusIIReader.SIDE_SET, 1, 1)</div><div>reader.SetObjectStatus(vtk.vtkExodusIIReader.NODE_SET, 0, 1)</div><div>reader.SetObjectStatus(vtk.vtkExodusIIReader.NODE_SET, 1, 1)</div><div><br></div><div>index = 1<br></div><div>for i in range(reader.GetOutput().GetNumberOfBlocks()):</div><div>    print 'reader.GetOutput().GetBlock(', i, ')', index</div><div>    index += 1</div><div>    for j in range(reader.GetOutput().GetBlock(i).GetNumberOfBlocks()):</div><div>        print '    reader.GetOutput().GetBlock(', i, ')', '.GetBlock(',j, ')', index</div><div>        index += 1</div><div><br></div><div>blocks0 = vtk.vtkExtractBlock()<br></div><div>blocks0.SetInputConnection(reader.GetOutputPort())</div><div>blocks0.AddIndex(8) # Set the indices here to view the various 'blocks'</div><div>blocks0.Update()</div></div><div><br></div><div><br></div><div>-Andrew</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jan 29, 2016 at 9:40 AM, David Thompson <span dir="ltr"><<a href="mailto:david.thompson@kitware.com" target="_blank">david.thompson@kitware.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Andrew,<br>
<br>
Since the Exodus reader generates a 2-level multiblock dataset, you should be able to use vtkExtractBlock to do everything you want; each Exodus block, side set, and node set should be placed in separate leaf-node of the VTK multiblock dataset. The reader also groups all of the leaf nodes that are Exodus blocks together, all of the leaf nodes that are side sets together, and so on -- this is what makes the reader's output a multi-level multiblock dataset.<br>
<br>
        David<br>
<div><div class="h5"><br>
> On Jan 29, 2016, at 11:12 AM, Slaughter, Andrew E <<a href="mailto:andrew.slaughter@inl.gov">andrew.slaughter@inl.gov</a>> wrote:<br>
><br>
> I have been working on this for sometime and still haven't come up with a solution, is there anyone who could help me find a solution, if it is possible?<br>
><br>
> I will provide another description of what I am trying to accomplish:<br>
><br>
> (1) I would like to read an ExodusII file that contains BLOCKS, SIDE_SETS, and NODE_SETS. When reading the file I want to enable everything (SetAllArrayStatus enabled for each type).<br>
><br>
> (2) I would like to extract the various parts (each BLOCK, SIDE_SET, NODE_SET) so that I can create a unique geometry->mapper->actor objects. This will enable me to show any combination of the parts, each with a different colormap, etc., but with only a single reader open.<br>
><br>
> I was able to get number (2) working using vtkExtractBlock for BLOCKS, but I can't figure out how to also extract out the SIDE_SET and NODE_SET. One idea I have was to use vtkExtractCell, but was unable to find the correct methodology for that idea.<br>
><br>
> Is this type extraction possible? If so, any hint in the correct direction would be greatly appreciated.<br>
><br>
> Thanks,<br>
> Andrew<br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
> On Tue, Dec 8, 2015 at 9:02 PM, Slaughter, Andrew E <<a href="mailto:andrew.slaughter@inl.gov">andrew.slaughter@inl.gov</a>> wrote:<br>
> I am creating a visualization tool for Exodus data capable of displaying any combination of nodeset, sideset, and subdomain information from an ExodusII file. I would like to be able to have a single reader object and multiple renderers for showing the nodesets, sidesets, and/or subdomains.<br>
><br>
> However, I am unable to figure out how to enable all of the above in vtkExodusIIReader but limit what is shown in the renderer.<br>
><br>
> I have attached a script (OSX with VTK6.3) that I would like to have display the nodesets on the left-side of the screen and the subdomains on the right, but I am unable to get the desired result.<br>
><br>
> I would appreciate any help that could be provided.<br>
><br>
> Thanks,<br>
> Andrew<br>
><br>
</div></div>> _______________________________________________<br>
> Powered by <a href="http://www.kitware.com" rel="noreferrer" target="_blank">www.kitware.com</a><br>
><br>
> Visit other Kitware open-source projects at <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.kitware.com_opensource_opensource.html&d=BQIFAg&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=h7heP8xwI1i_HikChvhFbEBurKirgfOCdwgBxB9lM8c&m=IzvSV8J8_UY6qt_rnM6naTRnGNk_stEmtgGBVwg5eCA&s=kl_8I3anInVvggyEaod6WrX7KeURy8h_C7mbTsZyOVI&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.kitware.com_opensource_opensource.html&d=BQIFAg&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=h7heP8xwI1i_HikChvhFbEBurKirgfOCdwgBxB9lM8c&m=IzvSV8J8_UY6qt_rnM6naTRnGNk_stEmtgGBVwg5eCA&s=kl_8I3anInVvggyEaod6WrX7KeURy8h_C7mbTsZyOVI&e=</a><br>
><br>
> Please keep messages on-topic and check the VTK FAQ at: <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__www.vtk.org_Wiki_VTK-5FFAQ&d=BQIFAg&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=h7heP8xwI1i_HikChvhFbEBurKirgfOCdwgBxB9lM8c&m=IzvSV8J8_UY6qt_rnM6naTRnGNk_stEmtgGBVwg5eCA&s=aiN7SFhWFQJe3Wj4LjzkP_bvTS-kyIr411JaPQb5Z6Y&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__www.vtk.org_Wiki_VTK-5FFAQ&d=BQIFAg&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=h7heP8xwI1i_HikChvhFbEBurKirgfOCdwgBxB9lM8c&m=IzvSV8J8_UY6qt_rnM6naTRnGNk_stEmtgGBVwg5eCA&s=aiN7SFhWFQJe3Wj4LjzkP_bvTS-kyIr411JaPQb5Z6Y&e=</a><br>
><br>
> Search the list archives at: <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__markmail.org_search_-3Fq-3Dvtkusers&d=BQIFAg&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=h7heP8xwI1i_HikChvhFbEBurKirgfOCdwgBxB9lM8c&m=IzvSV8J8_UY6qt_rnM6naTRnGNk_stEmtgGBVwg5eCA&s=_RvOk4cwnd3I4lZci8ZiXALDQ5S68i2l8goeAGzkXNs&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__markmail.org_search_-3Fq-3Dvtkusers&d=BQIFAg&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=h7heP8xwI1i_HikChvhFbEBurKirgfOCdwgBxB9lM8c&m=IzvSV8J8_UY6qt_rnM6naTRnGNk_stEmtgGBVwg5eCA&s=_RvOk4cwnd3I4lZci8ZiXALDQ5S68i2l8goeAGzkXNs&e=</a><br>
><br>
> Follow this link to subscribe/unsubscribe:<br>
> <a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__public.kitware.com_mailman_listinfo_vtkusers&d=BQIFAg&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=h7heP8xwI1i_HikChvhFbEBurKirgfOCdwgBxB9lM8c&m=IzvSV8J8_UY6qt_rnM6naTRnGNk_stEmtgGBVwg5eCA&s=wzPIk1KYN2NI_Inwd0yJDYEjbO7-t_4M7aAcPeej758&e=" rel="noreferrer" target="_blank">https://urldefense.proofpoint.com/v2/url?u=http-3A__public.kitware.com_mailman_listinfo_vtkusers&d=BQIFAg&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=h7heP8xwI1i_HikChvhFbEBurKirgfOCdwgBxB9lM8c&m=IzvSV8J8_UY6qt_rnM6naTRnGNk_stEmtgGBVwg5eCA&s=wzPIk1KYN2NI_Inwd0yJDYEjbO7-t_4M7aAcPeej758&e=</a><br>
<br>
</blockquote></div><br></div>