[cig-commits] r21064 - in seismo/3D/FAULT_SOURCE/branches/new_fault_db: . CUBIT
surendra at geodynamics.org
surendra at geodynamics.org
Wed Nov 21 14:32:54 PST 2012
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
---------------------
More information about the CIG-COMMITS
mailing list