1 |
|
2 |
######################################################## |
3 |
# |
4 |
# Copyright (c) 2003-2008 by University of Queensland |
5 |
# Earth Systems Science Computational Center (ESSCC) |
6 |
# http://www.uq.edu.au/esscc |
7 |
# |
8 |
# Primary Business: Queensland, Australia |
9 |
# Licensed under the Open Software License version 3.0 |
10 |
# http://www.opensource.org/licenses/osl-3.0.php |
11 |
# |
12 |
######################################################## |
13 |
|
14 |
__copyright__="""Copyright (c) 2003-2008 by University of Queensland |
15 |
Earth Systems Science Computational Center (ESSCC) |
16 |
http://www.uq.edu.au/esscc |
17 |
Primary Business: Queensland, Australia""" |
18 |
__license__="""Licensed under the Open Software License version 3.0 |
19 |
http://www.opensource.org/licenses/osl-3.0.php""" |
20 |
__url__="http://www.uq.edu.au/esscc/escript-finley" |
21 |
|
22 |
""" |
23 |
Author: John Ngui, john.ngui@uq.edu.au |
24 |
""" |
25 |
|
26 |
# Import the necessary modules. |
27 |
from esys.pyvisi import Scene, DataCollector, Contour, Camera |
28 |
from esys.pyvisi.constant import * |
29 |
import os |
30 |
|
31 |
PYVISI_EXAMPLE_MESHES_PATH = "data_meshes" |
32 |
PYVISI_EXAMPLE_IMAGES_PATH = "data_sample_images" |
33 |
X_SIZE = 400 |
34 |
Y_SIZE = 300 |
35 |
|
36 |
SCALAR_FIELD_POINT_DATA_1 = "lava" |
37 |
SCALAR_FIELD_POINT_DATA_2 = "talus" |
38 |
FILE_2D = "phi_talus_lava." |
39 |
|
40 |
IMAGE_NAME = "seriesofreads" |
41 |
JPG_RENDERER = Renderer.ONLINE_JPG |
42 |
|
43 |
# Create a Scene. |
44 |
s = Scene(renderer = JPG_RENDERER, num_viewport = 1, x_size = X_SIZE, |
45 |
y_size = Y_SIZE) |
46 |
|
47 |
# Create a DataCollector reading from a XML file. |
48 |
dc1 = DataCollector(source = Source.XML) |
49 |
dc1.setActiveScalar(scalar = SCALAR_FIELD_POINT_DATA_1) |
50 |
|
51 |
# Create a Contour. |
52 |
mosc1 = Contour(scene = s, data_collector = dc1, |
53 |
viewport = Viewport.SOUTH_WEST, lut = Lut.COLOR, cell_to_point = False, |
54 |
outline = True) |
55 |
mosc1.generateContours(0) |
56 |
|
57 |
# Create a second DataCollector reading from the same XML file |
58 |
# but specifying a different scalar field. |
59 |
dc2 = DataCollector(source = Source.XML) |
60 |
dc2.setActiveScalar(scalar = SCALAR_FIELD_POINT_DATA_2) |
61 |
|
62 |
# Create a second Contour. |
63 |
mosc2 = Contour(scene = s, data_collector = dc2, |
64 |
viewport = Viewport.SOUTH_WEST, lut = Lut.COLOR, cell_to_point = False, |
65 |
outline = True) |
66 |
mosc2.generateContours(0) |
67 |
|
68 |
# Create a Camera. |
69 |
cam1 = Camera(scene = s, viewport = Viewport.SOUTH_WEST) |
70 |
|
71 |
# Read in one file one after another and render the object. |
72 |
for i in range(99, 104): |
73 |
dc1.setFileName(file_name = os.path.join(PYVISI_EXAMPLE_MESHES_PATH, \ |
74 |
FILE_2D + "%04d.vtu") % i) |
75 |
dc2.setFileName(file_name = os.path.join(PYVISI_EXAMPLE_MESHES_PATH, \ |
76 |
FILE_2D + "%04d.vtu") % i) |
77 |
|
78 |
s.render(image_name = os.path.join(PYVISI_EXAMPLE_IMAGES_PATH, \ |
79 |
IMAGE_NAME + "%04d.jpg") % i) |