[cig-commits] [commit] master: small update the documentation of grdcoulomb.sh (85f235b)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/relax

On branch  : master
Link       : https://github.com/geodynamics/relax/compare/d80f922cd4775ef98aa4eccadd02e196d5db8410...85f235bc32a892310a95e9a9f885243c5c696f83

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

commit 85f235bc32a892310a95e9a9f885243c5c696f83
Author: Sylvain Barbot <sbarbot at ntu.edu.sg>
Date:   Fri Jul 4 00:23:18 2014 -0700

    small update the documentation of grdcoulomb.sh


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

85f235bc32a892310a95e9a9f885243c5c696f83
 util/extrude.sh    |  2 +-
 util/grdcoulomb.sh | 13 ++++++++-----
 2 files changed, 9 insertions(+), 6 deletions(-)

diff --git a/util/extrude.sh b/util/extrude.sh
index cf4b40b..21002a8 100755
--- a/util/extrude.sh
+++ b/util/extrude.sh
@@ -55,7 +55,7 @@ if (test "$hset" == "1"); then
 	exit
 fi
 
-echo "# nb   x1   x2   x3 length width thickness strike dip"
+echo "#  n   x1   x2   x3 length width thickness strike dip"
 while read line
 do
 	if [ "$line" == "" ]; then
diff --git a/util/grdcoulomb.sh b/util/grdcoulomb.sh
index 53d9147..cf38a18 100755
--- a/util/grdcoulomb.sh
+++ b/util/grdcoulomb.sh
@@ -62,7 +62,7 @@ if [ $# -lt "1" -o "$hset" == "1" ]; then
   index-tc.grd (Coulomb stress)
 
 description:
-  using the coordinate system
+  Using the coordinate system
 
              (x1, north)
             /
@@ -84,9 +84,11 @@ if [ $# -lt "1" -o "$hset" == "1" ]; then
      s = [cos(str), sin(str), 0]
      d = [+sin(str) cos(dip), -cos(str) cos(dip), -sin(dip)]
 
-  and form a right-handed coordinate system.
+  and form a right-handed coordinate system. The rake-direction vector is
 
-  the local traction components are:
+     b = [ cos(r) s(1) + sin(r) d(1), cos(r) s(2) + sin(r) d(2), sin(r) d(3) ]
+
+  The local traction components are:
 
      t1 = S11 n(1) + S12 n(2) + S13 n(3)
      t2 = S12 n(1) + S22 n(2) + S23 n(3)
@@ -98,13 +100,14 @@ if [ $# -lt "1" -o "$hset" == "1" ]; then
      ts = t(1) s(1) + t(2) s(2)
      td = t(1) d(1) + t(2) d(2) + t(3) d(3)
 
-  the Coulomb stress is
+  The Coulomb stress is
 
      tc = t(1) b(1) + t(2) b(2) + t(3) b(3) + mu * tn
 
   by convention, ts is positive for left lateral strike shear. dip direction
   is up dip. the 1, 2 and 3 directions are north, east and down, respectively.
-  the vector s, d, and n form a right-handed coordinate system.
+  The vector s, d, and n form a right-handed coordinate system.
+
 EOF
 
 	exit