[cig-commits] [commit] master: fix the polarity of some stress .grd outputs (c8834a2)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/relax

On branch  : master
Link       : https://github.com/geodynamics/relax/compare/db1b5c1faa54651a122f7782d0b1f2b0d621ce8d...3636c032c11009cef9eed4bf22bfbd4cdb1928e0

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

commit c8834a26818ff4a4aae288bab98e5c6770f58f63
Author: Sylvain Barbot <sbarbot at ntu.edu.sg>
Date:   Thu Feb 6 10:29:56 2014 +0800

    fix the polarity of some stress .grd outputs


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

c8834a26818ff4a4aae288bab98e5c6770f58f63
 src/export.f90 | 15 ++++++++++++---
 1 file changed, 12 insertions(+), 3 deletions(-)

diff --git a/src/export.f90 b/src/export.f90
index a982edd..618dbbf 100644
--- a/src/export.f90
+++ b/src/export.f90
@@ -1291,9 +1291,6 @@ END SUBROUTINE exportcreep_vtk
     CALL exportspatial(c2(:,:,int(oz/dx3)+1),sx1,sx2,temp2,doflip=.true.)
     CALL exportspatial(c3(:,:,int(oz/dx3)+1),sx1,sx2,temp3,doflip=.true.)
 
-    ! positive up
-    temp3=-temp3
-    
     pos=INDEX(wdir," ")
     WRITE (digit,'(I3.3)') i
     
@@ -1302,14 +1299,26 @@ END SUBROUTINE exportcreep_vtk
        file1=wdir(1:pos-1) // "/" // digit // "-north.grd"
        file2=wdir(1:pos-1) // "/" // digit // "-east.grd"
        file3=wdir(1:pos-1) // "/" // digit // "-up.grd"
+
+       ! positive up
+       temp3=-temp3
+    
     CASE (2) ! postseismic displacement
        file1=wdir(1:pos-1) // "/" // digit // "-relax-north.grd"
        file2=wdir(1:pos-1) // "/" // digit // "-relax-east.grd"
        file3=wdir(1:pos-1) // "/" // digit // "-relax-up.grd"
+
+       ! positive up
+       temp3=-temp3
+    
     CASE (3) ! equivalent body forces
        file1=wdir(1:pos-1) // "/" // digit // "-eqbf-north.grd"
        file2=wdir(1:pos-1) // "/" // digit // "-eqbf-east.grd"
        file3=wdir(1:pos-1) // "/" // digit // "-eqbf-up.grd"
+
+       ! positive up
+       temp3=-temp3
+    
     CASE (4) ! equivalent body forces
        file1=wdir(1:pos-1) // "/" // digit // "-s11.grd"
        file2=wdir(1:pos-1) // "/" // digit // "-s12.grd"