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

[ampuero at geodynamics.org]

Author: ampuero
Date: 2012-11-26 16:37:53 -0800 (Mon, 26 Nov 2012)
New Revision: 21073

Modified:
   seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_dynamic.f90
Log:
allow slip-strengthening friction

Modified: seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_dynamic.f90
===================================================================
--- seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_dynamic.f90	2012-11-26 23:55:49 UTC (rev 21072)
+++ seismo/3D/FAULT_SOURCE/branches/new_fault_db/src/fault_solver_dynamic.f90	2012-11-27 00:37:53 UTC (rev 21073)
@@ -684,8 +684,7 @@
   type(swf_type), intent(in) :: f
   real(kind=CUSTOM_REAL) :: mu(size(f%theta))
 
-  mu = f%mus -(f%mus-f%mud)/f%dc *f%theta
-  mu = max( mu, f%mud)
+  mu = f%mus -(f%mus-f%mud)* min(f%theta/f%dc, 1e0_CUSTOM_REAL)
 
 end function swf_mu