[cig-commits] r19877 - seismo/3D/FAULT_SOURCE/branches/new_fault_db/src

[surendra at geodynamics.org]

Author: surendra
Date: 2012-03-26 13:39:06 -0700 (Mon, 26 Mar 2012)
New Revision: 19877

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90
Log:
Generalized TPV16 coordinate locations for input file

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	2012-03-26 20:02:18 UTC (rev 19876)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90	2012-03-26 20:39:06 UTC (rev 19877)
@@ -249,6 +249,7 @@
        dyn_fc,swcd,cohes,tim_forcedRup
 
   integer :: iLoc
+  real(kind=CUSTOM_REAL) :: minX
 
 
   NAMELIST / INIT_STRESS / S1,S2,S3,n1,n2,n3
@@ -378,10 +379,12 @@
     allocate( bc%swf%theta(bc%nglob) )
     allocate( bc%swf%T(bc%nglob) )
     allocate( bc%swf%C(bc%nglob) )
-  
+
+    minX = minval(bc%coord(1,:))
+
     do iLoc=1,bc%nglob
   
-      ipar = minloc( (-24000.0+loc_str(:)-bc%coord(1,iLoc))**2 + (-loc_dip(:)-bc%coord(3,iLoc))**2 , 1) !iloc_dip is negative of Z-coord
+      ipar = minloc( (minX+loc_str(:)-bc%coord(1,iLoc))**2 + (-loc_dip(:)-bc%coord(3,iLoc))**2 , 1) !iloc_dip is negative of Z-coord
       bc%T0(3,iLoc) = -sigma0(ipar)
       bc%T0(1,iLoc) = tau0_str(ipar)
       bc%T0(2,iLoc) = tau0_dip(ipar)