pyvisi/examples/escript_scalarData2D.py

#!/usr/bin/env python


# $Id: escript_scalarData2D.py,v 1.3 2006/01/03 08:46:20 paultcochrane Exp $

from esys.escript import *
from esys import bruce

import numarray

from pyvisi import *
from pyvisi.renderers.vtk import *

tensorDomain = bruce.Rectangle(9,9,10,10)
tensorFunctionSpace=escript.ContinuousFunction(tensorDomain)
domainData = tensorFunctionSpace.getX()

# plotting scalar data in a 2D array
scalarData2D = sin(domainData[0])

scene = Scene()
plot = ContourPlot(scene)
plot.setData(scalarData2D)
scene.render(pause=True)

scene.save(fname="escript_scalarData2D_contourPlot.png", format="png")

scene = Scene()
plot = SurfacePlot(scene)
plot.setData(scalarData2D)
scene.render(pause=True)

scene.save(fname="escript_scalarData2D_surfacePlot.png", format="png")

# add SurfacePlot, MeshPlot etc here

# vim: expandtab shiftwidth=4:
1	#!/usr/bin/env python
2
3
4	# $Id: escript_scalarData2D.py,v 1.3 2006/01/03 08:46:20 paultcochrane Exp $
5
6	from esys.escript import *
7	from esys import bruce
8
9	import numarray
10
11	from pyvisi import *
12	from pyvisi.renderers.vtk import *
13
14	tensorDomain = bruce.Rectangle(9,9,10,10)
15	tensorFunctionSpace=escript.ContinuousFunction(tensorDomain)
16	domainData = tensorFunctionSpace.getX()
17
18	# plotting scalar data in a 2D array
19	scalarData2D = sin(domainData[0])
20
21	scene = Scene()
22	plot = ContourPlot(scene)
23	plot.setData(scalarData2D)
24	scene.render(pause=True)
25
26	scene.save(fname="escript_scalarData2D_contourPlot.png", format="png")
27
28	scene = Scene()
29	plot = SurfacePlot(scene)
30	plot.setData(scalarData2D)
31	scene.render(pause=True)
32
33	scene.save(fname="escript_scalarData2D_surfacePlot.png", format="png")
34
35	# add SurfacePlot, MeshPlot etc here
36
37	# vim: expandtab shiftwidth=4: