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

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)