1 # Copyright (C) 2010-2014 CEA/DEN, EDF R&D
3 # This library is free software; you can redistribute it and/or
4 # modify it under the terms of the GNU Lesser General Public
5 # License as published by the Free Software Foundation; either
6 # version 2.1 of the License, or (at your option) any later version.
8 # This library is distributed in the hope that it will be useful,
9 # but WITHOUT ANY WARRANTY; without even the implied warranty of
10 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
11 # Lesser General Public License for more details.
13 # You should have received a copy of the GNU Lesser General Public
14 # License along with this library; if not, write to the Free Software
15 # Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
17 # See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 if not ('servermanager' in dir()):
22 from pvsimple import *
26 def demo10(fname, impth):
27 """This method demonstrates the Contour filter."""
29 reader = ExodusIIReader(FileName=fname)
31 ##Show(reader, view=v)
32 contourFilter = Contour(reader)
33 contourFilter.ContourBy = 'Temp'
34 contourFilter.Isosurfaces
37 contourFilter.Isosurfaces = [400]
38 camera = GetActiveCamera()
39 camera.Position = [0., 60., -33]
42 dp2 = GetDisplayProperties(contourFilter)
43 contourFilter.PointData[:]
44 ##[Array: Temp, Array: V, Array: Pres, Array: AsH3, Array: GaMe3, Array: CH4, Array: H2]
46 temp = contourFilter.PointData[0]
48 ## (399.99999999999994, 400.00000000000006)
50 dp2.LookupTable = MakeBlueToRedLT(399.999,400)
51 dp2.ColorAttributeType = 'POINT_DATA'
52 dp2.ColorArrayName = 'Temp'
54 WriteImage(filename = (impth + "contourFilter_1.png"), view=v, Magnification=2)
56 pres = contourFilter.PointData[2]
58 ##(0.0074376024368978605, 0.02845284189526591)
60 dp2.LookupTable = MakeBlueToRedLT(0.007437602, 0.02845284)
61 dp2.ColorArrayName = 'Pres'
63 WriteImage(filename = (impth + "contourFilter_2.png"), view=v, Magnification=2)
65 dp1 = GetDisplayProperties(reader)
66 dp1.Representation = 'Wireframe'
67 dp1.LookupTable = dp2.LookupTable
68 dp1.ColorAttributeType = 'POINT_DATA'
69 dp1.ColorArrayName = 'Pres'
74 camera.Elevation(-150)
76 WriteImage(filename = (impth + "contourFilter_3.png"), view=v, Magnification=2)
79 scriptdir = inspect.getframeinfo(inspect.currentframe())[0]
80 testdir = os.path.dirname( os.path.abspath(scriptdir) )
81 pvdata = os.getenv("SAMPLES_SRC_DIR")
83 if __name__ == "__main__":
84 demo10(fname=pvdata+"/Data/disk_out_ref.ex2", impth=testdir+"/Pic/")