So you understand what I'm doing, below is the full script. I added those lines and now its giving one above the number of points. I have 800 points but it pints out 801
import numpy as np import pylab as pl try: paraview.simple except: from paraview.simple import * paraview.simple._DisableFirstRenderCameraReset() num_plots = 10 time = ['00010','00050','00090','00130','00184','00250','00300','00350','00400','00500'] lims = np.zeros(num_plots) pl.gca().set_color_cycle(['black', 'green', 'blue', 'yellow','purple','chocolate','gray','brown','cyan','gold']) nameoffile = 'TEST_nodal' for i in range(0,num_plots): pfile = time[i] nameof_file = nameoffile+'.'+pfile+'.vts' print nameof_file paraview.simple._DisableFirstRenderCameraReset() TEST_nodal_00001_vts = XMLStructuredGridReader(FileName=nameof_file ) TEST_nodal_00001_vts.UpdatePipeline() TEST_nodal_00001_vts.PointArrayStatus = ['Pressure'] ExtractSubset0 = ExtractSubset(guiName="ExtractSubset0",IncludeBoundary=0,VOI=[399, 399, 0, 0, 0, 799], SampleRateK=1, SampleRateJ=1, SampleRateI=1) ExtractSubset0.UpdatePipeline() data = servermanager.Fetch(ExtractSubset0) print ExtractSubset0.GetDataInformation().GetNumberOfPoints() numPoints=data.GetNumberOfPoints() print numPoints p = np.zeros(numPoints) coord = np.zeros((numPoints,3)) for j in range(0,numPoints): coord[j,:] = np.array(data.GetPoint(j)) p[j]=data.GetPointData().GetArray("Pressure").GetValue(j) pl.plot(coord[:,2], p,linewidth=2.0) pl.xlabel('Position Along Length (m)',size = 20) pl.ylabel('Pore pressure (MPa)', size = 20) pl.ticklabel_format(axis='y', style='sci', scilimits=(-2,2)) pl.legend(['t = 10 s','t = 50 s','t = 90 s', 't = 130 s', 't = 184 s', 't = 250 s', 't = 300 s', 't = 350', 't = 400 s', 't = 500 s'],'best',loc='upper left', numpoints=1, prop={'size':16}) pl.savefig('pressures.pdf') On Mon, Nov 3, 2014 at 9:00 AM, Utkarsh Ayachit <utkarsh.ayac...@kitware.com > wrote: > Hi, > > I wonder if the problem is with Fetch(). What does the following yeild? > > ExtractSubset0.UpdatePipeline() > print ExtractSubset0.GetDataInformation().GetNumberOfPoints() > > Utkarsh > > On Thu, Oct 30, 2014 at 5:29 PM, Chukwudi Chukwudozie > <cchu...@tigers.lsu.edu> wrote: > > Hi, > > > > I have a structured grid computational domain of size 800,2,800. I want > to > > extract the information along a line and so I do the following in my > python > > script. > > > > ExtractSubset0 = > > ExtractSubset(guiName="ExtractSubset0",IncludeBoundary=0,VOI=[399, 399, > 0, > > 0, 0, 799], SampleRateK=1, SampleRateJ=1, SampleRateI=1) > > data = servermanager.Fetch(ExtractSubset0) > > numPoints=data.GetNumberOfPoints() > > > > However, when I plot the information on that line, I see it has only > > extracted information between 0-399 in the z-coordinate direction. I even > > printed numPoints to confirm how much information is extracted and I get > > "numPoints=400". I am using PV 3.98.0. > > > > Any help on this will be appreciated. > > > > > > > > > > -- > > Chukwudi Chukwudozie > > cchu...@tigers.lsu.edu > > Craft and Hawkins Department of Petroleum Engineering > > Louisiana State University > > Patrick F. Taylor Hall (CEBA) > > Room 3429 > > Baton Rouge LA, 70803 > > Cell:3373546822 > > Office: 2255786064 > > > > _______________________________________________ > > 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://public.kitware.com/mailman/listinfo/paraview > > > -- Chukwudi Chukwudozie cchu...@tigers.lsu.edu Craft and Hawkins Department of Petroleum Engineering Louisiana State University Patrick F. Taylor Hall (CEBA) Room 3429 Baton Rouge LA, 70803 Cell:3373546822 Office: 2255786064
_______________________________________________ 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://public.kitware.com/mailman/listinfo/paraview