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

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

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

revision 3135 by ahallam, Thu Sep 2 03:28:39 2010 UTC revision 3195 by ahallam, Wed Sep 22 00:28:04 2010 UTC
# Line 25  from esys.pycad.gmsh import Design Line 25  from esys.pycad.gmsh import Design
25  from esys.finley import MakeDomain  from esys.finley import MakeDomain
26  from esys.escript.unitsSI import *  from esys.escript.unitsSI import *
27  from math import tan  from math import tan
 import os  
28    
29  #  #
30  #   input data:  #   input data:
# Line 37  fault_w_front=200*m Line 36  fault_w_front=200*m
36  fault_w_back=200*m  fault_w_back=200*m
37  fault_mid_front=2*km  fault_mid_front=2*km
38  fault_mid_back=4*km  fault_mid_back=4*km
39  fault_dip_front=45*DEG  fault_dip_front=30*DEG
40  fault_dip_back=40*DEG  fault_dip_back=30*DEG
41    
42  layers_left_at_front=[ [ 'limestone' ,  1100*m ] , [  'xx' , 150*m ] , ['shale',250*m] , [ 'limestone' ] ]  layers_left_at_front=[ [ 'limestone' ,  1100*m ] , [  'xx' , 150*m ] ,['shale',300*m],  [ 'limestone' ] ]
43  layers_right_at_front=[ [ 'limestone' , 500*m ], [ 'xx' , 150*m ], ['shale',250*m], ['limestone' ] ]  layers_right_at_front=[ [ 'limestone' , 500*m ], [ 'xx' , 150*m ],['shale',300*m], ['limestone' ] ]
44  slop=2*DEG  slop=2*DEG
45  #  #
46  # ====================================================================================  # ====================================================================================
# Line 176  FF.append(-PlaneSurface(CurveLoop(*tuple Line 175  FF.append(-PlaneSurface(CurveLoop(*tuple
175    
176    
177  # war 120  # war 120
178  des=Design(dim=3, order=1, element_size = 200*m, keep_files=True)  des=Design(dim=3, order=1, element_size = 40*m, keep_files=True)
179  des.addItems(*tuple(PS.values()))  des.addItems(*tuple(PS.values()))
180  des.addItems(PropertySet("fault",Volume(SurfaceLoop( *tuple(FF)))))  des.addItems(PropertySet("fault",Volume(SurfaceLoop( *tuple(FF)))))
181  des.setScriptFileName(os.path.join("data","example10","fault.geo"))  des.setMeshFileName("fault.msh")
 des.setMeshFileName(os.path.join("data","example10","fault.msh"))  
182  dom=MakeDomain(des)  dom=MakeDomain(des)
183  dom.write(os.path.join("data","example10","fault.fly"))  dom.write("fault.fly")
184    

Legend:
Removed from v.3135  
changed lines
  Added in v.3195

  ViewVC Help
Powered by ViewVC 1.1.26