/[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

revision 1058 by jongui, Fri Mar 23 04:06:56 2007 UTC revision 2420 by gross, Thu May 14 02:28: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 52  def wavePropagation(domain,h,tend,lam,mu Line 79  def wavePropagation(domain,h,tend,lam,mu
79     u_pc_z = u_pc[2]     u_pc_z = u_pc[2]
80    
81     # open file to save displacement at point source     # open file to save displacement at point source
82     u_pc_data=open('./data/U_pc.out','w')     u_pc_data=FileWriter('./data/U_pc.out')
83     u_pc_data.write("%f %f %f %f\n"%(t,u_pc_x,u_pc_y,u_pc_z))     u_pc_data.write("%f %f %f %f\n"%(t,u_pc_x,u_pc_y,u_pc_z))
84    
85     while t<tend:     while t<tend:
# Line 82  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.1058  
changed lines
  Added in v.2420

  ViewVC Help
Powered by ViewVC 1.1.26