/[escript]/trunk/finley/test/python/seismic_wave.py
ViewVC logotype

Diff of /trunk/finley/test/python/seismic_wave.py

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

revision 928 by gross, Mon Oct 9 03:50:20 2006 UTC revision 929 by gross, Wed Jan 17 07:41:13 2007 UTC
# Line 23  from esys.escript.linearPDEs import Line Line 23  from esys.escript.linearPDEs import Line
23  from esys.finley import Brick  from esys.finley import Brick
24  import time  import time
25    
26    WORKDIR="/raid2/lutz/waves/"
27  output=True  output=True
28  n_end=10000  n_end=10000
29    
30  resolution=1000.  # number of elements per m in the finest region  resolution=1000.  # number of elements per m in the finest region
31    resolution=400.  # number of elements per m in the finest region
32  o=1              # element order  o=1              # element order
33    
34  l=100000.           # width and length m (without obsorber)  l=100000.           # width and length m (without obsorber)
# Line 322  def wavePropagation(dom,rho,mu,lmbd,eta) Line 324  def wavePropagation(dom,rho,mu,lmbd,eta)
324       # ... save current acceleration in units of gravity and displacements       # ... save current acceleration in units of gravity and displacements
325       if output:       if output:
326            if t>=t_write:            if t>=t_write:
327               saveVTK("disp.%i.vtu"%n,displacement=u, amplitude=length(u))               saveVTK(WORKDIR+"disp.%i.vtu"%n_write,displacement=u, amplitude=length(u))
328               t_write+=dt_write               t_write+=dt_write
329               n_write+=1               n_write+=1
330       t+=dt       t+=dt

Legend:
Removed from v.928  
changed lines
  Added in v.929

  ViewVC Help
Powered by ViewVC 1.1.26