/[escript]/trunk/doc/examples/cookbook/example09b.py
ViewVC logotype

Diff of /trunk/doc/examples/cookbook/example09b.py

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

revision 3099 by ahallam, Sun Aug 22 04:35:15 2010 UTC revision 3387 by caltinay, Thu Nov 25 07:09:23 2010 UTC
# Line 28  __url__="https://launchpad.net/escript-f Line 28  __url__="https://launchpad.net/escript-f
28  #######################################################EXTERNAL MODULES  #######################################################EXTERNAL MODULES
29  from esys.escript import *  from esys.escript import *
30  from esys.finley import Rectangle  from esys.finley import Rectangle
31    from esys.weipa import saveVTK
32  import os  import os
33  # smoothing operator  # smoothing operator
34  from esys.escript.pdetools import Projector, Locator  from esys.escript.pdetools import Projector, Locator
35  from esys.escript.unitsSI import *  from esys.escript.unitsSI import *
36  import numpy as np  import numpy as np
37    import matplotlib
38    matplotlib.use('agg') #It's just here for automated testing
39    
40  import pylab as pl  import pylab as pl
41  import matplotlib.cm as cm  import matplotlib.cm as cm
42  from esys.escript.linearPDEs import LinearPDE  from esys.escript.linearPDEs import LinearPDE
# Line 64  mu2=vel2**2.*rho2/4.;  mu1=vel1**2.*rho1 Line 68  mu2=vel2**2.*rho2/4.;  mu1=vel1**2.*rho1
68  lam2=vel2**2.*rho2/2.; lam1=vel1**2.*rho1/2.  #lames constant  lam2=vel2**2.*rho2/2.; lam1=vel1**2.*rho1/2.  #lames constant
69    
70  # Time related variables.  # Time related variables.
71  tend=0.1    # end time  testing=True
72    if testing:
73        print 'The testing end time is currently selected. This severely limits the number of time iterations.'
74        print "Try changing testing to False for more iterations."
75        tend=0.001
76    else:
77        tend=0.1    # end time
78    
79  h=0.00005    # time step  h=0.00005    # time step
80  # data recording times  # data recording times
81  rtime=0.0 # first time to record  rtime=0.0 # first time to record

Legend:
Removed from v.3099  
changed lines
  Added in v.3387

  ViewVC Help
Powered by ViewVC 1.1.26