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

revision 3099 by ahallam, Sun Aug 22 04:35:15 2010 UTC revision 3387 by caltinay, Thu Nov 25 07:09:23 2010 UTC
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