/[escript]/trunk/doc/examples/usersguide/wave.py
ViewVC logotype

Diff of /trunk/doc/examples/usersguide/wave.py

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

trunk/doc/user/examples/wave.py revision 582 by lkettle, Wed Mar 8 05:54:06 2006 UTC trunk/doc/examples/wave.py revision 2344 by jfenwick, Mon Mar 30 02:13:58 2009 UTC
# Line 1  Line 1 
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__="https://launchpad.net/escript-finley"
21    
22    # You can shorten the execution time by reducing variable tend from 60 to 0.5
23    
24  from esys.escript import *  from esys.escript import *
25  from esys.escript.pdetools import Locator  from esys.escript.pdetools import Locator
26  from esys.escript.linearPDEs import LinearPDE  from esys.escript.linearPDEs import LinearPDE
27  from esys.finley import Brick  from esys.finley import Brick
28  from numarray import identity,zeros,ones  from numarray import identity,zeros,ones
29    
30    if not os.path.isdir("data"):
31       print "\nCreating subdirectory 'data'\n"
32       os.mkdir("data")
33    
34  ne=32          # number of cells in x_0 and x_1 directions  ne=32          # number of cells in x_0 and x_1 directions
35  width=10000.  # length in x_0 and x_1 directions  width=10000.  # length in x_0 and x_1 directions
36  lam=3.462e9  lam=3.462e9
# Line 42  def wavePropagation(domain,h,tend,lam,mu Line 69  def wavePropagation(domain,h,tend,lam,mu
69     t=0     t=0
70    
71     # define the location of the point source     # define the location of the point source
72     L=Locator(domain,xc)     L=Locator(domain,numarray.array(xc))
73     # find potential at point source     # find potential at point source
74     u_pc=L.getValue(u)     u_pc=L.getValue(u)
75     print "u at point charge=",u_pc     print "u at point charge=",u_pc
# Line 70  def wavePropagation(domain,h,tend,lam,mu Line 97  def wavePropagation(domain,h,tend,lam,mu
97       t+=h       t+=h
98       n+=1       n+=1
99       print n,"-th time step t ",t       print n,"-th time step t ",t
      L=Locator(domain,xc)  
100       u_pc=L.getValue(u)       u_pc=L.getValue(u)
101       print "u at point charge=",u_pc       print "u at point charge=",u_pc
102            
# Line 83  def wavePropagation(domain,h,tend,lam,mu Line 109  def wavePropagation(domain,h,tend,lam,mu
109    
110       # ... save current acceleration in units of gravity and displacements       # ... save current acceleration in units of gravity and displacements
111       if n==1 or n%10==0: saveVTK("./data/usoln.%i.vtu"%(n/10),acceleration=length(a)/9.81,       if n==1 or n%10==0: saveVTK("./data/usoln.%i.vtu"%(n/10),acceleration=length(a)/9.81,
112       displacement = length(u), Ux = u[0] )       displacement = length(u), tensor = stress, Ux = u[0] )
113    
114     u_pc_data.close()     u_pc_data.close()
115        

Legend:
Removed from v.582  
changed lines
  Added in v.2344

  ViewVC Help
Powered by ViewVC 1.1.26