[cig-commits] r18678 - in seismo/3D/FAULT_SOURCE/branches/new_fault_db: . src

Wed Jun 29 07:03:50 PDT 2011

Author: ampuero
Date: 2011-06-29 07:03:50 -0700 (Wed, 29 Jun 2011)
New Revision: 18678

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_object.f90
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90
Log:
cleaning up

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT
===================================================================

--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT	2011-06-29 12:54:24 UTC (rev 18677)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/README_SPECFEM3D_FAULT	2011-06-29 14:03:50 UTC (rev 18678)
@@ -1,14 +1,14 @@
 This package contains a version of the 3D spectral element code SPECFEM3D-SESAME (Komatitsch et al)
 modified by Galvez, Ampuero and Nissen-Meyer to model dynamic earthquake on non-planar faults
-with slip-weakening friction. The main modifications are two new modules (fault_object.f90 and
-fault_solver.f90). We also include some examples and Matlab function for post-processing and 
-visualization. 
+with slip-weakening friction. The main modifications are encapsulated in three new modules: 
+src/fault_object.f90, src/fault_solver.f90) and decompose_mesh_SCOTCH/fault_scotch.f90. 
+We also include some examples and Matlab functions for post-processing and visualization. 
 
 This is a preliminary version, still under development and testing. The format of inputs, outputs
 and features are subject to change.
 
 Details about the original package are described in its manual (manual_SPECFEM3D.pdf).
-The current README file describes how to install and run our modified version the code.
+The current README file describes how to install and run our modified version of the code.
 It is not a replacement for the original manual, which you should read first.
 
 This document assumes that the root directory containing the code is ~/SPECFEM3D.
@@ -136,8 +136,10 @@
   in directory DATABASES_MPI.
 
 
-6. If this is the first time you run a simulation with this mesh, edit DATA/Par_file to make sure that the variable DT (the time step)
-   is smaller than the "Maximum suggested time step" in OUTPUT_FILES/output_mesher.txt. ALso edit variables NSTEP and NTSTEP_* accordingly.
+6. If this is the first time you run a simulation with this mesh, set the parameter
+   DT (the time step) DATA/Par_file to a value smaller than the 
+   "Maximum suggested time step" in OUTPUT_FILES/output_mesher.txt. 
+   Also edit the variables NSTEP and NTSTEP_* in DATA/Par_file accordingly.
 
 
 7. Run the solver:
@@ -166,7 +168,7 @@
 
 	>> cp -r EXAMPLES/tpv5/DATA ./
 
-Then follow all the steps in section II above.
+3. follow all the steps in section II above.
 
 
 
@@ -193,14 +195,14 @@
 	Line NF+4: Slip velocity threshold below which frictional healing is set (friction coefficient
                is reset to its static value). If this value is negative healing is disabled.
 	Line NF+5: Slip velocity threshold to define the rupture front. Only used for outputs.
-       Line NF+6: Fault parameters id.  
-       Line NF+7: Initial fault stress for fault #1: 
+        Line NF+6: Fault parameters id.  
+        Line NF+7: Initial fault stress for fault #1: 
                	S1 = along-strike shear 
 			S2 = along-dip shear
 			S3 = normal stress (negative in compresion)
                   See figure 1 below. 
-       Line NF+8: Static and dynamic friction coefficients.
-       Repeat NF+6 to NF+8 for each fault.
+        Line NF+8: Static and dynamic friction coefficients.
+        Repeat NF+6 to NF+8 for each fault.
 
 DATA/FAULT/FAULT_STATIONS.in	Stations in the fault plane.
 	Line 1: number of stations.
@@ -253,6 +255,9 @@
 
 4. The face FAULT plane reference used in the computatinos takes fault_down as reference. 
 
+
+
+
 V. POST-PROCESSING AND VISUALIZATION
 -------------------------------------
 

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_object.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_object.f90	2011-06-29 12:54:24 UTC (rev 18677)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_object.f90	2011-06-29 14:03:50 UTC (rev 18678)
@@ -177,6 +177,8 @@
   do iflt=1,size(fault_db)
 ! TO DO : Percy, iface is computed with fault open , we have to change this and close  *_dummy
 !                 at this point.
+!jpa: we don't need to change it, the small fault opening does not affect this subroutine 
+!     (see get_element_face_id)
     call setup_iface(fault_db(iflt),nnodes_ext_mesh,nodes_coords_ext_mesh,nspec,nglob,ibool)
 
     ! saving gll indices for each fault face, needed for ibulks
@@ -515,7 +517,7 @@
     igll = 0
     do j=1,NGLLY
       do i=1,NGLLX
-        ! adds all gll points on thas face
+        ! adds all gll points on that face
         igll = igll + 1  
         ! stores weighted jacobian and normals
         fdb%jacobian2Dw(igll,ispec_flt) = jacobian2Dw_face(i,j)

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90	2011-06-29 12:54:24 UTC (rev 18677)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90	2011-06-29 14:03:50 UTC (rev 18678)
@@ -421,7 +421,7 @@
         stop 'bc_dynflt_3d::init_2d_distribution:: unknown shape'
     end select
    
-    where (b /= 0) a = b
+    where (b /= 0e0_CUSTOM_REAL) a = b
   enddo
     
   end subroutine init_2d_distribution

