/[escript]/trunk/doc/examples/geotutorial/steadystate_variablek.py
ViewVC logotype

Diff of /trunk/doc/examples/geotutorial/steadystate_variablek.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 2156 by gross, Mon Dec 15 05:09:02 2008 UTC revision 2549 by jfenwick, Mon Jul 20 06:43:47 2009 UTC
# Line 1  Line 1 
1  ########################################################  ########################################################
2  #  #
3  # Copyright (c) 2003-2008 by University of Queensland  # Copyright (c) 2003-2009 by University of Queensland
4  # Earth Systems Science Computational Center (ESSCC)  # Earth Systems Science Computational Center (ESSCC)
5  # http://www.uq.edu.au/esscc  # http://www.uq.edu.au/esscc
6  #  #
# Line 10  Line 10 
10  #  #
11  ########################################################  ########################################################
12    
13  __copyright__="""Copyright (c) 2003-2008 by University of Queensland  __copyright__="""Copyright (c) 2003-2009 by University of Queensland
14  Earth Systems Science Computational Center (ESSCC)  Earth Systems Science Computational Center (ESSCC)
15  http://www.uq.edu.au/esscc  http://www.uq.edu.au/esscc
16  Primary Business: Queensland, Australia"""  Primary Business: Queensland, Australia"""
17  __license__="""Licensed under the Open Software License version 3.0  __license__="""Licensed under the Open Software License version 3.0
18  http://www.opensource.org/licenses/osl-3.0.php"""  http://www.opensource.org/licenses/osl-3.0.php"""
19  __url__="http://www.uq.edu.au/esscc/escript-finley"  __url__="https://launchpad.net/escript-finley"
20    
21  # import tools  # import tools
22  from esys.escript import *  from esys.escript import *
# Line 49  mypde.setValue(A=k*kronecker(mydomain),Y Line 49  mypde.setValue(A=k*kronecker(mydomain),Y
49                  q=whereZero(x[1])+whereZero(x[1]-L1))                  q=whereZero(x[1])+whereZero(x[1]-L1))
50  # get temperature:  # get temperature:
51  T=mypde.getSolution()  T=mypde.getSolution()
52  # save to VTK for mayavi:  # save as VTK for visualisation:
53  saveVTK("u.xml",T=T)  saveVTK("u.vtu",T=T)
54    

Legend:
Removed from v.2156  
changed lines
  Added in v.2549

  ViewVC Help
Powered by ViewVC 1.1.26