[cig-commits] [commit] master: use a right-hand coordinate system for stress change calculations. (64bc932)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/relax

On branch  : master
Link       : https://github.com/geodynamics/relax/compare/f58c3e4a3b537ce548f0dc9b6301d68585ba7587...64bc932cf611ccf4683972590830ed18b974c6e9

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

commit 64bc932cf611ccf4683972590830ed18b974c6e9
Author: Sylvain Barbot <sbarbot at ntu.edu.sg>
Date:   Sun Apr 6 02:11:51 2014 +0800

    use a right-hand coordinate system for stress change calculations.


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

64bc932cf611ccf4683972590830ed18b974c6e9
 src/export.f90 | 12 ++++++------
 1 file changed, 6 insertions(+), 6 deletions(-)

diff --git a/src/export.f90 b/src/export.f90
index 5ca56a0..96c5264 100644
--- a/src/export.f90
+++ b/src/export.f90
@@ -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
@@ -2146,9 +2146,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