[cig-commits] r20781 - seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH

[ampuero at geodynamics.org]

Author: ampuero
Date: 2012-09-26 14:13:47 -0700 (Wed, 26 Sep 2012)
New Revision: 20781

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90
Log:
fixed bug with SMALLVALTOL in fault_scotch:lex_order

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90	2012-09-26 19:06:04 UTC (rev 20780)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90	2012-09-26 21:13:47 UTC (rev 20781)
@@ -264,7 +264,6 @@
   double precision, dimension(nspec) :: work,xp,yp,zp
   integer, dimension(nspec) :: ind,ninseg,iwork
   logical :: ifseg(nspec)
-  double precision :: UTM_X_MIN,UTM_X_MAX
   integer :: ispec,i,j
   integer :: ieoff,ilocnum,nseg,ioff,iseg
   double precision :: SMALLVALTOL
@@ -274,9 +273,7 @@
   zp=xyz_c(3,:)
 
   ! define geometrical tolerance based upon typical size of the model
-  UTM_X_MAX = maxval(xp)
-  UTM_X_MIN = minval(xp)
-  SMALLVALTOL = 1.d-10 * dabs(UTM_X_MAX - UTM_X_MIN)
+  SMALLVALTOL = 1.d-10 * maxval( maxval(xyz_c,2) - minval(xyz_c,2) )
 
   ! establish initial pointers
   do ispec=1,nspec