[cig-commits] r20933 - in seismo/3D/FAULT_SOURCE/branches/new_fault_db: decompose_mesh_SCOTCH src

[ampuero at geodynamics.org]

Author: ampuero
Date: 2012-10-25 18:05:03 -0700 (Thu, 25 Oct 2012)
New Revision: 20933

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90
Log:
minor fix in rsf_init

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-10-26 00:32:13 UTC (rev 20932)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/decompose_mesh_SCOTCH/fault_scotch.f90	2012-10-26 01:05:03 UTC (rev 20933)
@@ -261,7 +261,7 @@
   subroutine lex_order(xyz_c,loc,nnodes,nspec)
 
   integer, intent(in) :: nnodes,nspec
-  integer :: loc(nspec)
+  integer, intent(out) :: loc(nspec)
   double precision, intent(in) :: xyz_c(3,nspec)
    
   double precision, dimension(nspec) :: work,xp,yp,zp

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-10-26 00:32:13 UTC (rev 20932)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver.f90	2012-10-26 01:05:03 UTC (rev 20933)
@@ -263,9 +263,7 @@
     allocate(bc%MU(bc%nglob))
     if (RATE_AND_STATE) then
       allocate(bc%rsf)
-      call rsf_init(bc%rsf,bc%T0,bc%Fload,bc%coord,IIN_PAR)
-     ! WARNING: the line below is only valid for pure strike-slip faulting
-      bc%V(1,:) = bc%rsf%V_init
+      call rsf_init(bc%rsf,bc%T0,bc%V,bc%Fload,bc%coord,IIN_PAR)
     else 
       allocate(bc%swf)
       call swf_init(bc%swf,bc%MU,bc%coord,IIN_PAR)
@@ -703,10 +701,11 @@
 
 !=====================================================================
 
-subroutine rsf_init(f,T0,nucFload,coord,IIN_PAR)
+subroutine rsf_init(f,T0,V,nucFload,coord,IIN_PAR)
 
   type(rsf_type), intent(out) :: f
   real(kind=CUSTOM_REAL), intent(in) :: T0(:,:)
+  real(kind=CUSTOM_REAL), intent(inout) :: V(:,:)
   real(kind=CUSTOM_REAL), intent(in) :: coord(:,:)
   real(kind=CUSTOM_REAL), pointer :: nucFload(:)
   integer, intent(in) :: IIN_PAR
@@ -857,7 +856,7 @@
                         - f%f0 - f%a*log(f%V_init/f%V0) ) &
                       / f%b )
   else
-     f%theta =  f%a * log(TWO*V0/V_init * sinh(-sqrt(T0(1,:)**2+T0(2,:)**2)/T0(3,:)/f%a))
+     f%theta =  f%a * log(TWO*f%V0/f%V_init * sinh(-sqrt(T0(1,:)**2+T0(2,:)**2)/T0(3,:)/f%a))
   endif
 
  ! WARNING : ad hoc for SCEC benchmark TPV10x
@@ -868,6 +867,9 @@
   nucFload = Fload
   call init_2d_distribution(nucFload,coord,IIN_PAR,nFload)
 
+ ! WARNING: the line below is only valid for pure strike-slip faulting
+  V(1,:) = f%V_init
+
 end subroutine rsf_init
 
 !---------------------------------------------------------------------