[cig-commits] r18329 - seismo/3D/FAULT_SOURCE/branches/new_fault_db/CUBIT

[percygalvez at geodynamics.org]

Author: percygalvez
Date: 2011-05-06 09:38:04 -0700 (Fri, 06 May 2011)
New Revision: 18329

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/CUBIT/crack_fault.py
Log:
python script

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/CUBIT/crack_fault.py
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/CUBIT/crack_fault.py	2011-05-06 15:27:54 UTC (rev 18328)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/CUBIT/crack_fault.py	2011-05-06 16:38:04 UTC (rev 18329)
@@ -1,8 +1,6 @@
 #!/usr/bin/env python
 
 import cubit
-import boundary_definition
-import cubit2specfem3d 
 import math
 import os
 import sys
@@ -84,7 +82,7 @@
 cubit.cmd("merge all")
 cubit.cmd("mesh surface 1")
 cubit.cmd("mesh volume 1")
-cubit.cmd("unmerge surface 2 3")
+#cubit.cmd("unmerge surface 2 3")
 
 ########### Elements nodes    ###############
 fault_A_elements_up="sideset 200 surface 2"
@@ -92,27 +90,60 @@
 cubit.cmd(fault_A_elements_up)
 cubit.cmd(fault_A_elements_down)
 
+#### SIDESETS
 
-#nsets=cubit.get_sideset_id_list()
-### SIDESETS ARE ELEMENTSSSSSSSSSSSS #####
-# To get nodes on fault elments used : cubit.get_connectivity(entity_type,entity_id)
-# To get faces on fault elments used : cubit.get_sub_elements(entity_type,entity_id,dimension)
+sidesetlist = cubit.get_sideset_id_list()
 
+for sidesetid in sidesetlist:
+    sidesetname=cubit.get_exodus_entity_name('sideset',sidesetid)
+    print 'sideset '+sidesetname+' ids'
+    print 'sideset '+str(sidesetid)
 
-for nset in nsets:
-    el=cubit.get_sub_elements("hex",nset,2)
-    print el
-    
+
+for sidesetid in sidesetlist:
+    sidesetquads = cubit.get_sideset_quads(sidesetid)
+    sidesetsurfaces= cubit.get_sideset_surfaces(sidesetid)
+    print 'sideset quads :' , sidesetquads  # empty no quads ?
+    print 'sideset surface :' +str(sidesetsurfaces)
+
+
+#    sidesetdescription = cubit.get_exodus_entity_description('sideset',sidesetid)
+#    sidesetel = cubit.get_exodus_id('sideset',sidesetid)
+#    face = cubit.get_sub_elements('sideset',sidesetid,2)
+#    print 'sideset description :' + sidesetdescription # ....
+#    print 'sideset el:'      +str(sidesetel)         # -1
+#    print 'face : '+str(face)
+
+#    sideset_list = cubit.parse_cubit_list('node',sidesetid)
 #    name=cubit.get_exodus_entity_name("sideset",nset)
 #    id_sideset=cubit.get_exodus_id("hex",nset)
 #    print name
 #    print id_sideset
 
+### INTERFACE BLOCK
+fault_u = 300
+fault_d = 301
+side_u = "block "+str(fault_u)+" surface 2"
+side_d = "block "+str(fault_d)+" surface 3"
+cubit.cmd(side_u)
+cubit.cmd(side_d)
+el_u=cubit.get_block_faces(fault_u)
+el_d=cubit.get_block_faces(fault_d)
 
+for el in el_u:
+    print 'el_u :', el
 
+for el in el_d:
+    print 'el_d :', el
 
+#cubit.get_block_hexes(block_id)
 
 
 
 
 
+
+
+
+
+