/[escript]/trunk/dudley/test/python/subduction1_gen.py
ViewVC logotype

Diff of /trunk/dudley/test/python/subduction1_gen.py

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

revision 3773 by jfenwick, Mon Oct 11 01:48:14 2010 UTC revision 3774 by jfenwick, Wed Jan 18 06:29:34 2012 UTC
# Line 37  des=Design(dim=DIM, order=2, element_siz Line 37  des=Design(dim=DIM, order=2, element_siz
37    
38  if DIM ==2 :  if DIM ==2 :
39    if (DEPTH*cos(DIP) > 0.9 * LX):    if (DEPTH*cos(DIP) > 0.9 * LX):
40        raise ValueError,"X edge to short (maybe DIP to large.)"        raise ValueError("X edge to short (maybe DIP to large.)")
41    
42    p0=Point(DEPTH*cos(DIP),-DEPTH)    p0=Point(DEPTH*cos(DIP),-DEPTH)
43    p1=Point(LX,-DEPTH)    p1=Point(LX,-DEPTH)
# Line 54  if DIM ==2 : Line 54  if DIM ==2 :
54  else:  else:
55    
56    if (DEPTH*cos(DIP) > 0.9 * LX):    if (DEPTH*cos(DIP) > 0.9 * LX):
57        raise ValueError,"X edge to short (maybe DIP to large.)"        raise ValueError("X edge to short (maybe DIP to large.)")
58    if (DEPTH*cos(DIP)+LY*sin(STRIKE) > 0.9 * LX):    if (DEPTH*cos(DIP)+LY*sin(STRIKE) > 0.9 * LX):
59        raise ValueError,"X edge to short (maybe DIP to large.)"        raise ValueError("X edge to short (maybe DIP to large.)")
60    if (LY*sin(STRIKE) > 0.9 * LX):    if (LY*sin(STRIKE) > 0.9 * LX):
61        raise ValueError,"X edge to short (maybe DIP to large.)"        raise ValueError("X edge to short (maybe DIP to large.)")
62    
63    
64    p0=Point(DEPTH*cos(DIP),0.,-DEPTH)    p0=Point(DEPTH*cos(DIP),0.,-DEPTH)

Legend:
Removed from v.3773  
changed lines
  Added in v.3774

  ViewVC Help
Powered by ViewVC 1.1.26