[cig-commits] commit: make the loose rake constraint of afterslip optional and fixes the rake angle for dip-slip faults.

[Mercurial]

changeset:   7:ab4a669ee6d7
user:        Sylvain Barbot <sylbar.vainbot at gmail.com>
date:        Mon Apr 11 08:53:51 2011 -0700
files:       friction3d.f90 relax.f90
description:
make the loose rake constraint of afterslip optional and fixes the rake angle for dip-slip faults.


diff -r d75b50af47e5 -r ab4a669ee6d7 friction3d.f90
--- a/friction3d.f90	Mon Apr 11 08:33:10 2011 -0700
+++ b/friction3d.f90	Mon Apr 11 08:53:51 2011 -0700
@@ -335,8 +335,8 @@ CONTAINS
              ! effective shear stress on fault plane
              tau=taus+friction*taun-cohesion
 
-             ! rake direction test
-             IF (SUM(ts*r).LT.0.d0) CYCLE
+             ! rake direction test only if | rake | < 3*Pi
+             IF (SUM(ts*r).LT.0.d0 .AND. ABS(rake).LT.pi2*1.5d0) CYCLE
 
              ! warning for wrong input
              IF ((tau/tauc) .gt. 20) THEN
diff -r d75b50af47e5 -r ab4a669ee6d7 relax.f90
--- a/relax.f90	Mon Apr 11 08:33:10 2011 -0700
+++ b/relax.f90	Mon Apr 11 08:53:51 2011 -0700
@@ -156,6 +156,9 @@ PROGRAM relax
   !                  (02-28-11) - add constraints on the broad direction of 
   !                               afterslip, export faults to GMT xy format
   !                               and allow scaling of computed time steps.
+  !                  (04-11-11) - make the constraint on the broad direction
+  !                               of afterslip optional and fixes the case
+  !                               of dip-slip afterslip.
   !-----------------------------------------------------------------------
 
   USE green
@@ -1372,7 +1375,14 @@ CONTAINS
                 WRITE (0,'("error in input file: plane index misfit")')
                 STOP 1
              END IF
-             
+
+             ! modify rake for consistency with slip model
+             IF (n(k)%rake .GE. 0.d0) THEN
+                n(k)%rake=n(k)%rake-180.d0
+             ELSE             
+                n(k)%rake=n(k)%rake+180.d0
+             END IF
+
              ! comply to Wang's convention
              CALL wangconvention(dummy,n(k)%x,n(k)%y,n(k)%z,&
                   n(k)%length,n(k)%width,n(k)%strike,n(k)%dip,n(k)%rake,rot)