[cig-commits] commit: fix the export of Coulomb stress on receiver faults in .xy GMT format

Mercurial hg at geodynamics.org
Wed Oct 17 19:39:50 PDT 2012


changeset:   150:b04ce120d6d3
tag:         tip
user:        Sylvain Barbot <sbarbot at caltech.edu>
date:        Wed Oct 17 19:39:46 2012 -0700
files:       src/export.f90
description:
fix the export of Coulomb stress on receiver faults in .xy GMT format


diff -r 8e013252c3a4 -r b04ce120d6d3 src/export.f90
--- a/src/export.f90	Wed Oct 17 01:01:55 2012 -0700
+++ b/src/export.f90	Wed Oct 17 19:39:46 2012 -0700
@@ -1806,10 +1806,10 @@ END SUBROUTINE exportcreep_vtk
        taucoulomb=taus+friction*taun
 
        WRITE (15,'("> -Z",5ES11.2)') taus, taun, taucoulomb, taustrike, taudip
-       WRITE (15,'(3ES11.2)') x1-d(1)*W/2-s(1)*L/2, x2-d(2)*W/2-s(2)*L/2
-       WRITE (15,'(3ES11.2)') x1-d(1)*W/2+s(1)*L/2, x2-d(2)*W/2+s(2)*L/2
-       WRITE (15,'(3ES11.2)') x1+d(1)*W/2+s(1)*L/2, x2+d(2)*W/2+s(2)*L/2
-       WRITE (15,'(3ES11.2)') x1+d(1)*W/2-s(1)*L/2, x2+d(2)*W/2-s(2)*L/2
+       WRITE (15,'(3ES11.2)') x2-d(2)*W/2-s(2)*L/2, x1-d(1)*W/2-s(1)*L/2
+       WRITE (15,'(3ES11.2)') x2-d(2)*W/2+s(2)*L/2, x1-d(1)*W/2+s(1)*L/2
+       WRITE (15,'(3ES11.2)') x2+d(2)*W/2+s(2)*L/2, x1+d(1)*W/2+s(1)*L/2
+       WRITE (15,'(3ES11.2)') x2+d(2)*W/2-s(2)*L/2, x1+d(1)*W/2-s(1)*L/2
 
     END DO
 



More information about the CIG-COMMITS mailing list