[cig-commits] [commit] master: export stress in the dip direction (rake=90 degrees). (f642e54)

cig_noreply at geodynamics.org cig_noreply at geodynamics.org
Sat Apr 5 08:57:01 PDT 2014


Repository : ssh://shell.geodynamics.org/relax

On branch  : master
Link       : https://github.com/geodynamics/relax/compare/924da9babd73caf2a86a6e5fcaa2eb85d620c27b...f642e5469b9fd9a355d4162fa424d663d2407f59

>---------------------------------------------------------------

commit f642e5469b9fd9a355d4162fa424d663d2407f59
Author: Sylvain Barbot <sbarbot at ntu.edu.sg>
Date:   Sat Apr 5 23:56:36 2014 +0800

    export stress in the dip direction (rake=90 degrees).


>---------------------------------------------------------------

f642e5469b9fd9a355d4162fa424d663d2407f59
 src/export.f90 | 24 ++++++++++++------------
 1 file changed, 12 insertions(+), 12 deletions(-)

diff --git a/src/export.f90 b/src/export.f90
index 618dbbf..5ca56a0 100644
--- a/src/export.f90
+++ b/src/export.f90
@@ -1781,9 +1781,9 @@ END SUBROUTINE exportcreep_vtk
        sdip=sin(dip)
  
        ! surface normal vector components
-       n(1)=+cdip*cstrike
-       n(2)=-cdip*sstrike
-       n(3)=-sdip
+       n(1)=-cdip*cstrike
+       n(2)=+cdip*sstrike
+       n(3)=+sdip
 
        ! strike-slip unit direction
        s(1)=sstrike
@@ -1791,9 +1791,9 @@ END SUBROUTINE exportcreep_vtk
        s(3)=0._8
 
        ! dip-slip unit direction
-       d(1)=+cstrike*sdip
-       d(2)=-sstrike*sdip
-       d(3)=+cdip
+       d(1)=-cstrike*sdip
+       d(2)=+sstrike*sdip
+       d(3)=-cdip
 
        ! traction vector
        t=lsig .tdot. n
@@ -1922,9 +1922,9 @@ END SUBROUTINE exportcreep_vtk
        sdip=sin(dip)
  
        ! surface normal vector components
-       n(1)=+cdip*cstrike
-       n(2)=-cdip*sstrike
-       n(3)=-sdip
+       n(1)=-cdip*cstrike
+       n(2)=+cdip*sstrike
+       n(3)=+sdip
 
        ! strike-slip unit direction
        s(1)=sstrike
@@ -1932,9 +1932,9 @@ END SUBROUTINE exportcreep_vtk
        s(3)=0._8
 
        ! dip-slip unit direction
-       d(1)=+cstrike*sdip
-       d(2)=-sstrike*sdip
-       d(3)=+cdip
+       d(1)=-cstrike*sdip
+       d(2)=+sstrike*sdip
+       d(3)=-cdip
 
        ! traction vector
        t=lsig .tdot. n



More information about the CIG-COMMITS mailing list