[cig-commits] r20774 - seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tohoku

[percygalvez at geodynamics.org]

Author: percygalvez
Date: 2012-09-25 06:48:07 -0700 (Tue, 25 Sep 2012)
New Revision: 20774

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tohoku/mesh_japan.py
Log:
opening tohoku crack with cubit

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-09-24 19:01:14 UTC (rev 20773)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/EXAMPLES/tohoku/mesh_japan.py	2012-09-25 13:48:07 UTC (rev 20774)
@@ -7,7 +7,7 @@
 import sys
 from save_fault_nodes_elements import *
 from absorbing_boundary import *
-from funcions import m2km
+from functions import m2km
 
 cubit.cmd('reset')
 
@@ -26,13 +26,18 @@
                                 # avoid fault segmentation fault errors
 
 ########### Slab elements and nodes ###############
-
 ########  SLAB  ################################################
 os.system('mkdir -p MESH') 
 
 Au = [42]  # A_up
 Ad = [36]  # A_down
 
+####### Fault opening #############################################
+cubit.cmd('set node constraint off')
+cubit.cmd('node in surf 42 move X 0 Y 0 Z 0.001')   # delta = 2e-3 km. in Z direction.
+cubit.cmd('node in surf 36 move X 0 Y 0 Z -0.001')  # In general the shift should be normal to the fault.
+
+####################################################################
 faultA = fault_input(1,Au,Ad)
 quads_Aup,quads_Adp = save_cracks(faultA.name,faultA.surface_u,faultA.surface_d)
 #Unpacking list.