[cig-commits] r18665 - seismo/3D/FAULT_SOURCE/trunk/EXAMPLES/splay_faults

Tue Jun 28 08:01:34 PDT 2011

Author: percygalvez
Date: 2011-06-28 08:01:34 -0700 (Tue, 28 Jun 2011)
New Revision: 18665

Removed:
   seismo/3D/FAULT_SOURCE/trunk/EXAMPLES/splay_faults/tsunami.py
Log:
deleteing tsunami.py

Deleted: seismo/3D/FAULT_SOURCE/trunk/EXAMPLES/splay_faults/tsunami.py
===================================================================

--- seismo/3D/FAULT_SOURCE/trunk/EXAMPLES/splay_faults/tsunami.py	2011-06-28 14:47:58 UTC (rev 18664)
+++ seismo/3D/FAULT_SOURCE/trunk/EXAMPLES/splay_faults/tsunami.py	2011-06-28 15:01:34 UTC (rev 18665)
@@ -1,201 +0,0 @@
-#!/usr/bin/env python
-
-import cubit
-import boundary_definition
-import cubit2specfem3d 
-import math
-import os
-import sys
-import numarray
-
-cubit.cmd('reset')
-
-km = 1000
-y_vert=0*km
-h = 0 #0.1#(h=1)
-h1 = 0 #0.1#(h1=0.1)
-# dipping angle 1
-d = math.pi/180
-dip1=7.5*d
-# dipping angle 2
-dip2=15*d
-
-####  initializing coordinates x,y,z
-x=[]     # fault
-y=[]
-z=[]
-
-xbulk=[] # bulk
-ybulk=[]
-zbulk=[]
-
-xmid=[]
-ymid=[]  # midpoints.
-zmid=[]
-
-# Bulk points ###
-xbulk.append(-10*km)    #x1
-xbulk.append(240*km)    #x2  
-xbulk.append(240*km)    #x3
-xbulk.append(-10*km)    #x4
-
-zbulk.append(0*km)   #y1
-zbulk.append(0*km)   #y2   
-zbulk.append(-50*km) #y3
-zbulk.append(-50*km) #y4
-
-ybulk=[y_vert]*4
-
-### CRACKS ##########
-x.append(0*km)                      #x5  
-x.append(4*km/math.tan(dip1))       #x6  = x11 (triple joint)
-x.append(x[1]+8*km/math.tan(dip2))  #x7
-x.append(x[2])                      #x8  = x7 (kink point)
-x.append(90*km)                     #x9
-x.append(90*km+h1)                  #x10 = x9+h1 (surface point, fault C)
-x.append(x[1])                      #x11 = x6
-x.append(227*km)                    #x12
-x.append(x[7]+h1)                   #x13 = x12+h1 (surface point, fault B)
-
-z.append(-30*km)     #z5
-z.append(-26*km)     #z6  
-z.append(-18*km)     #z7 (kink point)
-z.append(z[2]-h1)    #z8 = z7-h1 (kink point)
-z.append(0)          #z9
-z.append(0)          #z10 = z9 
-z.append(-26*km)  #z11 = z6
-z.append(0)          #z12
-z.append(0)          #z13
-
-
-y=[y_vert]*9
-
-ymid=[y_vert]*6
-
-#####
-
-####################  bulk ###########################################
-for i in range(len(xbulk)): 
-   vert="create vertex x "+str(xbulk[i])+" y "+str(ybulk[i])+" z "+str(zbulk[i]) 
-   cubit.cmd(vert) 
-
-################  Loading fault points profile#############################
-for i in range(len(x)):
-  vert="create vertex x "+str(x[i])+" y "+str(y[i])+" z "+str(z[i])
-  cubit.cmd(vert)
-
-################  Loading fault points profile#############################
-for i in range(len(xmid)):
-  vert="create vertex x "+str(xmid[i])+" y "+str(ymid[i])+" z "+str(zmid[i])
-  cubit.cmd(vert)
-
-
-################ creating fault domains #################################
-
-bulk1="create curve vertex 1 9"   #c1
-bulk2="create curve vertex 10 12" #c2
-bulk3="create curve vertex 13 2"  #c3   
-bulk4="create curve vertex 2 3"   #c4
-bulk5="create curve vertex 3 4"   #c5
-bulk6="create curve vertex 4 1"   #c6
-
-cubit.cmd(bulk1)
-cubit.cmd(bulk2)
-cubit.cmd(bulk3)
-cubit.cmd(bulk4)
-cubit.cmd(bulk5)
-cubit.cmd(bulk6)
-
-fault_up_A1="create curve spline vertex 5 6"       #c7
-fault_up_A2="create curve spline vertex 6 12"      #c8
-
-fault_down_A1="create curve spline vertex 5 11"    #c9
-fault_down_A2="create curve spline vertex 11 13"   #c10 
-
-fault_up_BC1="create curve vertex 9 7"     #c11
-fault_up_BC2="create curve vertex 7 6"     #c12
-fault_down_BC1="create curve vertex 10 8"  #c13
-fault_down_BC2="create curve vertex 8 6"   #c14
- 
-cubit.cmd(fault_up_A1) 
-cubit.cmd(fault_up_A2) 
-cubit.cmd(fault_down_A1) 
-cubit.cmd(fault_down_A2) 
-
-cubit.cmd(fault_up_BC1) 
-cubit.cmd(fault_up_BC2) 
-cubit.cmd(fault_down_BC1) 
-cubit.cmd(fault_down_BC2) 
-
-#surface="create surface curve 1 11 12 7 9 10 3 4 5 6"
-surface="create surface curve 1 11 12 10 3 4 5 6"
-
-cubit.cmd(surface)
-
-surface="create surface curve 2 13 14 8"
-cubit.cmd(surface)
-
-
-cubit.cmd("sweep surface 1 vector 0 1 0 distance "+str(90*km)) 
-cubit.cmd("sweep surface 2 vector 0 1 0 distance "+str(90*km))
-cubit.cmd("sweep curve 7 vector 0 1 0 distance "+str(90*km))
-cubit.cmd("sweep curve 9 vector 0 1 0 distance "+str(90*km))
-# fault_up_A1
-cubit.cmd("sweep curve 7 vector 0 1 0 distance")
-#fault_down_A1
-cubit.cmd("sweep curve 9 vector 0 1 0 distance")
-### Coarsening ########
-xcoorse=[]
-ycoorse=[]
-zcoorse=[]
-
-kx=250*km
-ky=90*km
-kz=50*km
-
-# FAULT_A1_up 
-#cubit.cmd("curve 7 merge off")
-#cubit.cmd("curve 24 merge off")
-cubit.cmd("surface 3 merge off ")
-
-# FAULT_A1_down 
-#cubit.cmd("curve 9 merge off")
-#cubit.cmd("curve 26 merge off")
-cubit.cmd("surface 3 merge off ")
-
-# IMPRINTING 
-cubit.cmd("imprint all")
-# MERGING
-cubit.cmd("merge all")
-# Sticking fault borders ##
-#### Fault A.
-#  UP
-cubit.cmd("merge curve 8 with curve 10 force") # init
-cubit.cmd("merge curve 22 with curve 34 force") # final
-## DOWN.
-#cubit.cmd("merge curve 7 with curve 9 force")
-#cubit.cmd('merge curve 24 with curve 26 force')
-##### Fault BC
-#### Up
-cubit.cmd("merge curve 11 with curve 13 force")
-cubit.cmd("merge curve 30 with curve 42 force")
-#### DOWN
-cubit.cmd("merge curve 12 with curve 14 force")
-cubit.cmd("merge curve 28 with curve 35 force")
-
-# Mesh coinciding fault_A upper and lower boundaries .
-# Meshing the whole domain.
-#cubit.cmd("surface 13 18 size "+str(elementsize))
-#cubit.cmd("volume 1 2 3 size "+str(elementsize))
-#cubit.cmd("surface 13 18 scheme pave")
-#cubit.cmd("mesh surface 13 18 ")
-#cubit.cmd("mesh volume 1 2 3")
-#cubit.cmd("unmerge surface 2 3")
-
-########## Coarse shield #######
-
-
-
-
-
-

