[cig-commits] commit: change the surface traction units from force/area/rigidity to force/area.

Mercurial hg at geodynamics.org
Sat Aug 10 06:41:41 PDT 2013


changeset:   210:18e421c92579
tag:         tip
user:        Sylvain Barbot <sbarbot at ntu.edu.sg>
date:        Sat Aug 10 21:41:20 2013 +0800
files:       src/elastic3d.f90 src/input.f90
description:
change the surface traction units from force/area/rigidity to force/area.


diff -r 9dcb3ece04f7 -r 18e421c92579 src/elastic3d.f90
--- a/src/elastic3d.f90	Sat Aug 10 21:19:38 2013 +0800
+++ b/src/elastic3d.f90	Sat Aug 10 21:41:20 2013 +0800
@@ -2093,7 +2093,7 @@
 
              amp=omega((x1-x)/L,beta)* &
                  omega((x2-y)/W,beta)* &
-                 mu*e%l(i)%slip
+                 e%l(i)%slip
 
              period=e%l(i)%period
 
diff -r 9dcb3ece04f7 -r 18e421c92579 src/input.f90
--- a/src/input.f90	Sat Aug 10 21:19:38 2013 +0800
+++ b/src/input.f90	Sat Aug 10 21:41:20 2013 +0800
@@ -1239,7 +1239,7 @@
                STAT=iostatus)
           IF (iostatus>0) STOP "could not allocate the load list"
           PRINT 2000
-          PRINT '(a)',"t3 in units of force/surface/rigidity, positive down"
+          PRINT '(a)',"t3 in units of force/surface, positive down"
           PRINT '(a)',"T>0 for t3 sin(2pi/T+phi), T<=0 for t3 H(t)"
           PRINT '(a)',"no.       xs       ys   length    width       t3        T      phi"
           PRINT 2000



More information about the CIG-COMMITS mailing list