[Paraview] Internal volume vs. external surface via python interface?

[Utkarsh Ayachit]

Try this:

> Clip(Input=servermanager.OutputPort(ll9cellf523_ncdf,1))

I've added a bug for the tracing issue:
http://www.paraview.org/Bug/view.php?id=10651

Utkarsh


On Fri, Apr 30, 2010 at 1:09 PM, Greg Schussman
<schussman at slac.stanford.edu> wrote:
> Hi, Everyone.
>
> I am trying to display a slice from a volume dataset, but when I use the
> python interface, all I get is a slice of the external surface.  Within
> paraview's gui, there is no problem selecting "External Surface" vs.
> "Internal Volume".
>
> I have tried comparing python traces of slicing entire volume vs. exterior
> surface, but the traces are, for all practical purposes, identical.  I also
> tried looking at a trace of changing the slice's input, but the trace seems
> to have no relevant information.
>
> I tried comparing trace state, with no important difference found.   I tried
> comparing saved state (.pvsm) files, but saw no difference there except for
> id kinds of numbers.
>
> I think that ll9cellf523_ncdf (see below) has two output ports, where 0 is
> the surface and 1 is the entire volume, but I'm not sure of that.   From
> "Paraview Scripting with Python" (August 2007), page 15, it looks like I
> should be able to:
>
>      Slice1.Input = OutputPort(ll9cellf523_ncdf, 1)
>
> But that gives
>
>     NameError: name 'OutputPort' is not defined
>
> Trying either
>
>     Slice1.Input = ll9cellf523_ncdf.Port(0)
>     Slice1.Input = ll9cellf523_ncdf.Port(1)
>
> runs, but the slice is still only slicing the external surface.
>
> My test script is below.
>
> Thanks for any help.
>
> Greg
>
> #!/usr/bin/python
>
> try: paraview.simple
> except: from paraview.simple import *
>
> ll9cellf523_ncdf = SLACDataReader(MeshFileName='./ll-9cell-f523.ncdf' )
> ll9cellf523_ncdf.ModeFileName = ['./mode0.l0.R2.457087E+09I5.096430E+04.m2']
> ll9cellf523_ncdf.ReadInternalVolume = 1
>
> a3_efield_PVLookupTable = \
>    GetLookupTableForArray( "efield", 3,
>                            RGBPoints=[0.0,   0.0, 0.0, 1.0,
>                                       19.0,  1.0, 0.0, 0.0],
>                            VectorMode='Magnitude',
>                            ColorSpace='HSV',
>                            ScalarRangeInitialized=1.0 )
>
> RenderView1 = GetRenderView()
> RenderView1.CameraPosition = [-1.1, 0.0, 0.7]
> RenderView1.CameraFocalPoint = [0.0, 0.0, 0.7]
> RenderView1.CameraClippingRange = [2.68, 3.33]
>
> a3_efield_PVLookupTable.RGBPoints = [0.0,   0.0, 0.0, 1.0,
>                                      19.0,  1.0, 0.0, 0.0]
>
> Slice1 = Slice( SliceType="Plane" )
> Slice1.SliceType.Origin = [-0.023, 0.0, 0.7]
> Slice1.SliceOffsetValues = [0.0]
> Slice1.SliceType = "Plane"
>
> DataRepresentation1 = Show()
> DataRepresentation1.ColorArrayName = 'efield'
> DataRepresentation1.LookupTable = a3_efield_PVLookupTable
> DataRepresentation1.ColorAttributeType = 'POINT_DATA'
>
> ###########################################################
>
> UpdatePipeline()
>
> view            = GetActiveView()
> view.Background = [1, 1, 1]
> view.ViewSize   = [800, 300]
>
> WriteImage("test.png")
>
> p.s., I actually want to do something else (use a calculator to get efield
> magnitude, and then plot isocontour surfaces from that), but this slice
> example seemed simplest, and I'm assuming that once that works, the other
> stuff will work in the same way.
>
>
>
>
>
> _______________________________________________
> 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
>
>