[cig-commits] r20914 - in seismo/3D/FAULT_SOURCE/branches/new_fault_db: CUBIT EXAMPLES/splay_faults EXAMPLES/tohoku EXAMPLES/tpv102 EXAMPLES/tpv103 EXAMPLES/tpv15.crack EXAMPLES/tpv16.crack EXAMPLES/tpv5.crack
surendra at geodynamics.org
surendra at geodynamics.org
Wed Oct 24 13:54:19 PDT 2012
Author: surendra
Date: 2012-10-24 13:54:19 -0700 (Wed, 24 Oct 2012)
New Revision: 20914
Modified:
seismo/3D/FAULT_SOURCE/branches/new_fault_db/CUBIT/save_fault_nodes_elements.py
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/splay_faults.py
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tohoku/mesh_japan.py
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv102/TPV102.py
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/TPV103.py
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv15.crack/tpv15.py
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv15.crack/tpv15_coarsening.py
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv15.crack/tpv15_coarsening_base_chunk.py
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv16.crack/TPV16.py
seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/tpv5.py
Log:
factored the Cubit python interface functions for faults into a single one (and replaced in all examples by a single call to it)
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/CUBIT/save_fault_nodes_elements.py
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/CUBIT/save_fault_nodes_elements.py 2012-10-24 20:04:50 UTC (rev 20913)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/CUBIT/save_fault_nodes_elements.py 2012-10-24 20:54:19 UTC (rev 20914)
@@ -10,6 +10,21 @@
self.surface_d = surface_d
self.name = 'MESH/fault_file_'+str(id)+'.dat'
+ quads_Aup,quads_Adp = save_cracks(self.name,self.surface_u,self.surface_d)
+ #Unpacking list.
+ quads_Au=unpack_list(quads_Aup)
+ quads_Ad=unpack_list(quads_Adp)
+
+ print 'len(Au):',len(quads_Au)
+ print 'len(Ad):',len(quads_Ad)
+
+ if not (len(quads_Au)==len(quads_Ad)):
+ print 'Number of elements for each fauld side up and down do not concide'
+ sys.exit('goodbye')
+
+ save_elements_nodes(self.name,quads_Au,quads_Ad)
+
+
def save_cracks(name,list_surface_up,list_surface_down):
quads_fault_up = []
quads_fault_down = []
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/splay_faults.py
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/splay_faults.py 2012-10-24 20:04:50 UTC (rev 20913)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/splay_faults/splay_faults.py 2012-10-24 20:54:19 UTC (rev 20914)
@@ -167,38 +167,12 @@
Au = [8,15] #face_up
Ad = [6,7] #face_down
faultA = fault_input(1,Au,Ad)
-quads_Aup,quads_Adp = save_cracks(faultA.name,faultA.surface_u,faultA.surface_d)
-#Unpacking list.
-quads_Au=unpack_list(quads_Aup)
-quads_Ad=unpack_list(quads_Adp)
-print 'len(Au):',len(quads_Au)
-print 'len(Ad):',len(quads_Ad)
-
-if not (len(quads_Au)==len(quads_Ad)):
- print 'Number of elements for each fauld side up and down do not concide'
- sys.exit('goodbye')
-
-save_elements_nodes(faultA.name,quads_Au,quads_Ad)
-
########## FAULT BC ##############################################################
BCu = [9,10] #face_up
BCd = [14,17] #face_down
faultBC = fault_input(2,BCu,BCd)
-quads_BCup,quads_BCdp = save_cracks(faultBC.name,faultBC.surface_u,faultBC.surface_d)
-#Unpacking list.
-quads_BCu=unpack_list(quads_BCup)
-quads_BCd=unpack_list(quads_BCdp)
-print 'len(BCu):',len(quads_BCu)
-print 'len(BCd):',len(quads_BCd)
-
-if not (len(quads_BCu)==len(quads_BCd)):
- print 'Number of elements for each fauld side up and down do not concide'
- sys.exit('goodbye')
-
-save_elements_nodes(faultBC.name,quads_BCu,quads_BCd)
-
### Exporting the mesh to cubit.
boundary_definition.entities=['face']
boundary_definition.define_bc(boundary_definition.entities,parallel=True)
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tohoku/mesh_japan.py
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tohoku/mesh_japan.py 2012-10-24 20:04:50 UTC (rev 20913)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tohoku/mesh_japan.py 2012-10-24 20:54:19 UTC (rev 20914)
@@ -39,20 +39,7 @@
####################################################################
faultA = fault_input(1,Au,Ad)
-quads_Aup,quads_Adp = save_cracks(faultA.name,faultA.surface_u,faultA.surface_d)
-#Unpacking list.
-quads_Au=unpack_list(quads_Aup)
-quads_Ad=unpack_list(quads_Adp)
-print 'len(Au):',len(quads_Au)
-print 'len(Ad):',len(quads_Ad)
-
-if not (len(quads_Au)==len(quads_Ad)):
- print 'Number of elements for each fauld side up and down do not concide'
- sys.exit('goodbye')
-
-save_elements_nodes(faultA.name,quads_Au,quads_Ad)
-
## FOR THE BULK (Seismic wave propagation part for SESAME)
####### defining absorbing-boundary surface...
xmin = [133,114,124,97,129,106]
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv102/TPV102.py
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv102/TPV102.py 2012-10-24 20:04:50 UTC (rev 20913)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv102/TPV102.py 2012-10-24 20:54:19 UTC (rev 20914)
@@ -53,18 +53,8 @@
Ad = [3] # A_down
faultA = fault_input(1,Au,Ad)
-quads_Aup,quads_Adp = save_cracks(faultA.name,faultA.surface_u,faultA.surface_d)
-#Unpacking list.
-quads_Au=unpack_list(quads_Aup)
-quads_Ad=unpack_list(quads_Adp)
-print 'len(Au):',len(quads_Au)
-print 'len(Ad):',len(quads_Ad)
-if not (len(quads_Au)==len(quads_Ad)):
- print 'Number of elements for each fauld side up and down do not concide'
-
-save_elements_nodes(faultA.name,quads_Au,quads_Ad)
@@ -72,5 +62,3 @@
-
-
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/TPV103.py
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/TPV103.py 2012-10-24 20:04:50 UTC (rev 20913)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv103/TPV103.py 2012-10-24 20:54:19 UTC (rev 20914)
@@ -22,7 +22,7 @@
zmax = [8,15]
zmin = [10,14]
entities=['face']
-define_bc(entities,xmin,xmax,ymin,ymax,zmin,zmax)
+define_boundaries(entities,xmin,xmax,ymin,ymax,zmin,zmax)
cubit.cmd('block 1 name "elastic 1" ')
cubit.cmd('block 1 attribute count 5')
@@ -52,23 +52,11 @@
Ad = [3] # A_down
faultA = fault_input(1,Au,Ad)
-quads_Aup,quads_Adp = save_cracks(faultA.name,faultA.surface_u,faultA.surface_d)
-#Unpacking list.
-quads_Au=unpack_list(quads_Aup)
-quads_Ad=unpack_list(quads_Adp)
-print 'len(Au):',len(quads_Au)
-print 'len(Ad):',len(quads_Ad)
-if not (len(quads_Au)==len(quads_Ad)):
- print 'Number of elements for each fauld side up and down do not concide'
-
-save_elements_nodes(faultA.name,quads_Au,quads_Ad)
-
-
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv15.crack/tpv15.py
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv15.crack/tpv15.py 2012-10-24 20:04:50 UTC (rev 20913)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv15.crack/tpv15.py 2012-10-24 20:54:19 UTC (rev 20914)
@@ -141,39 +141,13 @@
Ad = [5,14] # A_down
faultA = fault_input(1,Au,Ad)
-quads_Aup,quads_Adp = save_cracks(faultA.name,faultA.surface_u,faultA.surface_d)
-#Unpacking list.
-quads_Au=unpack_list(quads_Aup)
-quads_Ad=unpack_list(quads_Adp)
-print 'len(Au):',len(quads_Au)
-print 'len(Ad):',len(quads_Ad)
-
-if not (len(quads_Au)==len(quads_Ad)):
- print 'Number of elements for each fauld side up and down do not concide'
- sys.exit('goodbye')
-
-save_elements_nodes(faultA.name,quads_Au,quads_Ad)
-
####### Branch Fault##############################################################################
Bu = [4] # B_up
Bd = [15] # B_down
faultB = fault_input(2,Bu,Bd)
-quads_Bup,quads_Bdp = save_cracks(faultB.name,faultB.surface_u,faultB.surface_d)
-#Unpacking list.
-quads_Bu=unpack_list(quads_Bup)
-quads_Bd=unpack_list(quads_Bdp)
-print 'len(Au):',len(quads_Bu)
-print 'len(Ad):',len(quads_Bd)
-
-if not (len(quads_Bu)==len(quads_Bd)):
- print 'Number of elements for each fauld side up and down do not concide'
- sys.exit('goodbye')
-
-save_elements_nodes(faultB.name,quads_Bu,quads_Bd)
-
## FOR THE BULK (Seismic wave propagation part for SESAME)
####### This is boundary_definition.py of GEOCUBIT
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv15.crack/tpv15_coarsening.py
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv15.crack/tpv15_coarsening.py 2012-10-24 20:04:50 UTC (rev 20913)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv15.crack/tpv15_coarsening.py 2012-10-24 20:54:19 UTC (rev 20914)
@@ -228,39 +228,13 @@
Ad = [5,14] # A_down
faultA = fault_input(1,Au,Ad)
-quads_Aup,quads_Adp = save_cracks(faultA.name,faultA.surface_u,faultA.surface_d)
-#Unpacking list.
-quads_Au=unpack_list(quads_Aup)
-quads_Ad=unpack_list(quads_Adp)
-print 'len(Au):',len(quads_Au)
-print 'len(Ad):',len(quads_Ad)
-
-if not (len(quads_Au)==len(quads_Ad)):
- print 'Number of elements for each fauld side up and down do not concide'
- sys.exit('goodbye')
-
-save_elements_nodes(faultA.name,quads_Au,quads_Ad)
-
####### Branch Fault##############################################################################
Bu = [4] # B_up
Bd = [15] # B_down
faultB = fault_input(2,Bu,Bd)
-quads_Bup,quads_Bdp = save_cracks(faultB.name,faultB.surface_u,faultB.surface_d)
-#Unpacking list.
-quads_Bu=unpack_list(quads_Bup)
-quads_Bd=unpack_list(quads_Bdp)
-print 'len(Au):',len(quads_Bu)
-print 'len(Ad):',len(quads_Bd)
-
-if not (len(quads_Bu)==len(quads_Bd)):
- print 'Number of elements for each fauld side up and down do not concide'
- sys.exit('goodbye')
-
-save_elements_nodes(faultB.name,quads_Bu,quads_Bd)
-
## FOR THE BULK (Seismic wave propagation part for SESAME)
####### This is boundary_definition.py of GEOCUBIT
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv15.crack/tpv15_coarsening_base_chunk.py
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv15.crack/tpv15_coarsening_base_chunk.py 2012-10-24 20:04:50 UTC (rev 20913)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv15.crack/tpv15_coarsening_base_chunk.py 2012-10-24 20:54:19 UTC (rev 20914)
@@ -235,39 +235,13 @@
Ad = [5,14] # A_down
faultA = fault_input(1,Au,Ad)
-quads_Aup,quads_Adp = save_cracks(faultA.name,faultA.surface_u,faultA.surface_d)
-#Unpacking list.
-quads_Au=unpack_list(quads_Aup)
-quads_Ad=unpack_list(quads_Adp)
-print 'len(Au):',len(quads_Au)
-print 'len(Ad):',len(quads_Ad)
-
-if not (len(quads_Au)==len(quads_Ad)):
- print 'Number of elements for each fauld side up and down do not concide'
- sys.exit('goodbye')
-
-save_elements_nodes(faultA.name,quads_Au,quads_Ad)
-
####### Branch Fault##############################################################################
Bu = [4] # B_up
Bd = [15] # B_down
faultB = fault_input(2,Bu,Bd)
-quads_Bup,quads_Bdp = save_cracks(faultB.name,faultB.surface_u,faultB.surface_d)
-#Unpacking list.
-quads_Bu=unpack_list(quads_Bup)
-quads_Bd=unpack_list(quads_Bdp)
-print 'len(Au):',len(quads_Bu)
-print 'len(Ad):',len(quads_Bd)
-
-if not (len(quads_Bu)==len(quads_Bd)):
- print 'Number of elements for each fauld side up and down do not concide'
- sys.exit('goodbye')
-
-save_elements_nodes(faultB.name,quads_Bu,quads_Bd)
-
## FOR THE BULK (Seismic wave propagation part for SESAME)
####### This is boundary_definition.py of GEOCUBIT
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv16.crack/TPV16.py
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv16.crack/TPV16.py 2012-10-24 20:04:50 UTC (rev 20913)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv16.crack/TPV16.py 2012-10-24 20:54:19 UTC (rev 20914)
@@ -14,7 +14,7 @@
from save_fault_nodes_elements import *
from absorbing_boundary import *
-cubit.cmd('playback "TPV16__tmp1.jou" ')
+cubit.cmd('playback "TPV16.jou" ')
os.system('mkdir -p MESH')
@@ -24,28 +24,16 @@
Ad = [13] # A_down
faultA = fault_input(1,Au,Ad)
-quads_Aup,quads_Adp = save_cracks(faultA.name,faultA.surface_u,faultA.surface_d)
-#Unpacking list.
-quads_Au=unpack_list(quads_Aup)
-quads_Ad=unpack_list(quads_Adp)
-print 'len(Au):',len(quads_Au)
-print 'len(Ad):',len(quads_Ad)
-if not (len(quads_Au)==len(quads_Ad)):
- print 'Number of elements for each fauld side up and down do not concide'
-
-save_elements_nodes(faultA.name,quads_Au,quads_Ad)
-
-
xmin = [4]
xmax = [6]
ymin = [17,10]
ymax = [12,15]
zmax = [9,14]
-zminx = [11,16]
+zmin = [11,16]
entities=['face']
-define_bc(entities,xmin,xmax,ymin,ymax,zmin,zmax)
+define_boundarie(entities,xmin,xmax,ymin,ymax,zmin,zmax)
cubit.cmd('block 1 name "elastic 1" ')
cubit.cmd('block 1 attribute count 5')
Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/tpv5.py
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/tpv5.py 2012-10-24 20:04:50 UTC (rev 20913)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tpv5.crack/tpv5.py 2012-10-24 20:54:19 UTC (rev 20914)
@@ -100,54 +100,11 @@
os.system('mkdir -p MESH')
-fault_u = 2 # fault surface up : surface 2
-fault_d = 3 # fault surface down : surface 3
-txt =''
+Au = 2
+Ad = 3
-fault_file=open('MESH/fault_file_1.dat','w')
-list_hex=cubit.parse_cubit_list('hex','all')
+faultA = fault_input(1,Au,Ad)
-quads_fault_u = cubit.get_surface_quads(fault_u)
-quads_fault_d = cubit.get_surface_quads(fault_d)
-
-# TO DO : stop python properly in case fault nodes at both sides
-# do not match.
-if len(quads_fault_u) != len(quads_fault_d):
- stop
-# Writting number of elements at both sides to make
-# double sure everything is going fine .
-txt='%10i %10i\n' % (len(quads_fault_u),len(quads_fault_d))
-fault_file.write(txt)
-
-dic_quads_fault_u = dict(zip(quads_fault_u,quads_fault_u))
-dic_quads_fault_d = dict(zip(quads_fault_d,quads_fault_d))
-
-
-# FAULT SIDE DOWN
-for h in list_hex:
- faces = cubit.get_sub_elements('hex',h,2)
- for f in faces:
- if dic_quads_fault_d.has_key(f):
- nodes=cubit.get_connectivity('Face',f)
-# print 'h,fault nodes side down :',h,nodes[0],nodes[1],nodes[2],nodes[3]
- txt='%10i %10i %10i %10i %10i\n' % (h,nodes[0],\
- nodes[1],nodes[2],nodes[3])
- fault_file.write(txt)
-
-# FAULT SIDE UP
-for h in list_hex:
- faces = cubit.get_sub_elements('hex',h,2)
- for f in faces:
- if dic_quads_fault_u.has_key(f):
- nodes=cubit.get_connectivity('Face',f)
-# print 'h,fault nodes side up :',h,nodes[0],nodes[1],nodes[2],nodes[3]
- txt='%10i %10i %10i %10i %10i\n' % (h,nodes[0],\
- nodes[1],nodes[2],nodes[3])
- fault_file.write(txt)
-
-# CLOSING FAULT FILE
-fault_file.close()
-
# FOR THE BULK (Seismic wave propagation part for SESAME)
###### This is boundary_definition.py of GEOCUBIT
More information about the CIG-COMMITS
mailing list