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

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

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

revision 2662 by gross, Tue Sep 8 07:11:12 2009 UTC revision 2663 by gross, Tue Sep 15 06:04:44 2009 UTC
# Line 24  from esys.escript import * Line 24  from esys.escript import *
24  from esys.escript.linearPDEs import LinearPDE  from esys.escript.linearPDEs import LinearPDE
25  from esys.escript.models import FaultSystem  from esys.escript.models import FaultSystem
26  from esys.finley import Rectangle  from esys.finley import Rectangle
27    from esys.escript.unitsSI import DEG
28  #... set some parameters ...  #... set some parameters ...
29  lam=1.  lam=1.
30  mu=1  mu=1
31  slip_max=1.  slip_max=1.
 fault_start=[0.5,0.25]  
 fault_end=[0.5,0.75]  
32    
33  mydomain = Rectangle(l0=1.,l1=1.,n0=16, n1=16)  # n1 need to be multiple of 4!!!  mydomain = Rectangle(l0=1.,l1=1.,n0=16, n1=16)  # n1 need to be multiple of 4!!!
34  # .. create the fault system  # .. create the fault system
35  fs=FaultSystem(dim=2)  fs=FaultSystem(dim=2)
36  fs.addFault(top=[fault_start, fault_end], tag=1)  fs.addFault(V0=[0.5,0.25], strikes=90*DEG, ls=0.5, tag=1)
37  # ... create a slip distribution on the fault:  # ... create a slip distribution on the fault:
38  p, m=fs.getParametrization(mydomain.getX(),tag=1)  p, m=fs.getParametrization(mydomain.getX(),tag=1)
39  p0,p1= fs.getW0Range(tag=1)  p0,p1= fs.getW0Range(tag=1)

Legend:
Removed from v.2662  
changed lines
  Added in v.2663

  ViewVC Help
Powered by ViewVC 1.1.26