[cig-commits] r21064 - in seismo/3D/FAULT_SOURCE/branches/new_fault_db: . CUBIT

[surendra at geodynamics.org]

Author: surendra
Date: 2012-11-21 14:32:54 -0800 (Wed, 21 Nov 2012)
New Revision: 21064

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/CUBIT/absorbing_boundary.py
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT
Log:
Added documentation related to Cubit-python scripts. Added back the accidentally deleted function that needs to setup faults

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/CUBIT/absorbing_boundary.py
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/CUBIT/absorbing_boundary.py	2012-11-21 09:01:28 UTC (rev 21063)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/CUBIT/absorbing_boundary.py	2012-11-21 22:32:54 UTC (rev 21064)
@@ -312,6 +312,21 @@
          build_block_side(topo,entity+'_topo',obj=entity)
 
 
+def define_boundaries(entities,xmin,xmax,ymin,ymax,zmin,zmax):
+     bottom=zmin
+     topo=zmax
+     v_list,name_list=define_block()
+     build_block(v_list,name_list)
+     print entities
+     for entity in entities:
+         print "##entity: "+str(entity)
+         build_block_side(xmin,entity+'_abs_xmin',obj=entity)
+         build_block_side(xmax,entity+'_abs_xmax',obj=entity)
+         build_block_side(ymin,entity+'_abs_ymin',obj=entity)
+         build_block_side(ymax,entity+'_abs_ymax',obj=entity)
+         build_block_side(bottom,entity+'_abs_bottom',obj=entity)
+         build_block_side(topo,entity+'_topo',obj=entity)
+
 def define_bc_topo(entities,self):
      # Temporal : Variable zmin should be obtained automatically. 
      xmin = self.xmin

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT	2012-11-21 09:01:28 UTC (rev 21063)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT	2012-11-21 22:32:54 UTC (rev 21064)
@@ -21,6 +21,7 @@
 	Download and updates
 	Installation
 	Mesh generation with split nodes
+	Cubit-python scripts for fault
 	Running a simulation
 	Examples
 	Sign convention for fault quantities
@@ -143,8 +144,19 @@
 
 For more complicated meshes, please refer to the Cubit scripts (*.jou and *.py) in EXAMPLES.
 
+CUBIT-PYTHON SCRIPTS FOR FAULT
+------------------------------
 
+For each fault, one needs to make a list of all surfaces that form each edge of the fault and call the function fault_input with a number corresponding to the fault along with the list of edges of both faults. The number given to fault_input should start from unity going up by one for each additional fault.
 
+For the example is above section, this is how one would call the function fault_input
+Au = [168]   # A_up
+Ad = [160]  # A_down
+faultA = fault_input(1,Au,Ad)
+
+The above lines generate a one file per fault "fault_file_1.dat", in additional to the regular files written from Cubit for bulk mesh.  The file numbers should increase sequentially for meshes with multiple faults. 
+
+
 RUNNING A SIMULATION
 ---------------------