/[escript]/trunk/pycad/examples/faults.py
ViewVC logotype

Diff of /trunk/pycad/examples/faults.py

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

trunk/pycad/examples/structure.py revision 932 by gross, Fri Jan 19 09:27:15 2007 UTC trunk/pycad/examples/faults.py revision 955 by gross, Mon Feb 12 03:54:54 2007 UTC
# Line 1  Line 1 
1  """  """
2  seismic wave propagation domain  block with some fualts
3    
4  @var __author__: name of author  @var __author__: name of author
5  @var __licence__: licence agreement  @var __licence__: licence agreement
# Line 21  __date__="$Date$" Line 21  __date__="$Date$"
21    
22  from esys.pycad import *  from esys.pycad import *
23  from esys.pycad.gmsh import Design  from esys.pycad.gmsh import Design
24    from esys.finley import MakeDomain
25    
26  l=100000.           # width and length m (without obsorber)  l=100000.           # width and length m (without obsorber)
27  h=30000.            # width and length m (without obsorber)  h=30000.            # width and length m (without obsorber)
# Line 47  m3=Line(q3,q4) Line 48  m3=Line(q3,q4)
48  m4=Line(q4,q1)  m4=Line(q4,q1)
49  c2=CurveLoop(m1,m2,m3,m4)  c2=CurveLoop(m1,m2,m3,m4)
50  c2.setLocalScale(0.1)  c2.setLocalScale(0.1)
51  dsgn=Design(element_size=h/3)  dsgn=Design(element_size=h/5)
52  dsgn.addItems(Volume(SurfaceLoop(*tuple(b.getSurfaces()+[PlaneSurface(c1),PlaneSurface(c2)]))))  dsgn.addItems(Volume(SurfaceLoop(*tuple(b.getSurfaces()+[PlaneSurface(c1),PlaneSurface(c2)]))))
53  dsgn.setScriptFileName("test.geo")  dsgn.setScriptFileName("faults.geo")
54  dsgn.setMeshFileName("test.msh")  dsgn.setMeshFileName("faults.msh")
55  print dsgn.getMeshHandler()  dom=MakeDomain(dsgn,integrationOrder=-1, reducedIntegrationOrder=-1, optimizeLabeling=True)
56    dom.write("faults.fly")
Legend:
Removed from v.932  
changed lines
  Added in v.955
  ViewVC Help
Powered by ViewVC 1.1.26