/[escript]/trunk/doc/examples/pyvisi/velocityonplanecut.py
ViewVC logotype

Annotation of /trunk/doc/examples/pyvisi/velocityonplanecut.py

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1153 - (hide annotations)
Thu May 17 06:33:54 2007 UTC (13 years, 9 months ago) by jongui
File MIME type: text/x-python
File size: 1342 byte(s)
Made some minor modifications and fixed minor bugs. Most tests appear to pass on the Altix although there are two tests which occasionally throw an error related to memory corruption.
1 ksteube 1147 # Import the necessary modules
2     from esys.pyvisi import Scene, DataCollector, VelocityOnPlaneCut, Camera
3     from esys.pyvisi.constant import *
4 jongui 1148 import os
5 ksteube 1147
6 jongui 1148 PYVISI_EXAMPLE_MESHES_PATH = "data_meshes"
7     PYVISI_EXAMPLE_IMAGES_PATH = "data_sample_images"
8 ksteube 1147 X_SIZE = 400
9     Y_SIZE = 400
10    
11     VECTOR_FIELD_CELL_DATA = "velocity"
12     FILE_3D = "interior_3D.xml"
13     IMAGE_NAME = "velocity.jpg"
14     JPG_RENDERER = Renderer.ONLINE_JPG
15    
16     # Create a Scene with four viewports
17     s = Scene(renderer = JPG_RENDERER, num_viewport = 1, x_size = X_SIZE,
18     y_size = Y_SIZE)
19    
20     # Create a DataCollector reading from a XML file.
21     dc1 = DataCollector(source = Source.XML)
22 jongui 1148 dc1.setFileName(file_name = os.path.join(PYVISI_EXAMPLE_MESHES_PATH, FILE_3D))
23 ksteube 1147 dc1.setActiveVector(vector = VECTOR_FIELD_CELL_DATA)
24    
25     # Create VelocityOnPlaneCut.
26     vopc1 = VelocityOnPlaneCut(scene = s, data_collector = dc1,
27     viewport = Viewport.SOUTH_WEST, color_mode = ColorMode.VECTOR,
28     arrow = Arrow.THREE_D, lut = Lut.COLOR, cell_to_point = False,
29     outline = True)
30     vopc1.setScaleFactor(scale_factor = 0.5)
31     vopc1.setPlaneToXY(offset = 0.5)
32     vopc1.setRatio(2)
33     vopc1.randomOn()
34    
35     # Create a Camera.
36 jongui 1148 c1 = Camera(scene = s, viewport = Viewport.SOUTH_WEST)
37 ksteube 1147 c1.isometricView()
38     c1.elevation(angle = -20)
39    
40     # Render the object.
41 jongui 1148 s.render(image_name = os.path.join(PYVISI_EXAMPLE_IMAGES_PATH, IMAGE_NAME))
42 ksteube 1147

  ViewVC Help
Powered by ViewVC 1.1.26