[cig-commits] commit: fixes the man page to specify convention about zero period for surface loads.

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


changeset:   211:eccc38abafa0
tag:         tip
user:        Sylvain Barbot <sbarbot at ntu.edu.sg>
date:        Sat Aug 10 21:59:03 2013 +0800
files:       man/man1/relax.1 src/relax.f90
description:
fixes the man page to specify convention about zero period for surface loads.


diff -r 18e421c92579 -r eccc38abafa0 man/man1/relax.1
--- a/man/man1/relax.1	Sat Aug 10 21:41:20 2013 +0800
+++ b/man/man1/relax.1	Sat Aug 10 21:59:03 2013 +0800
@@ -270,7 +270,7 @@
 
 # nb x1 x2 length width t3 T phi
 
-where nb is an index running from 1 to nsl, x1, x2, length and width define the position and dimension of the load, t3 is in units of stress (force/surface), positive down, and T can be a period (T > 0 implies stress=t3*sin(2 pi/T + phi) or not (T <= 0 implies stress = t3 H(t), with H(t) the Heaviside function).
+where nb is an index running from 1 to nsl, x1, x2, length and width define the position and dimension of the load, t3 is in units of stress (force/surface), positive down, and T can be a period (T != 0 implies stress=t3*sin(2 pi/T + phi) or not (T = 0 implies stress = t3 H(t), with H(t) the Heaviside function).
 
 .TP
 .B time of next event (te)
diff -r 18e421c92579 -r eccc38abafa0 src/relax.f90
--- a/src/relax.f90	Sat Aug 10 21:41:20 2013 +0800
+++ b/src/relax.f90	Sat Aug 10 21:59:03 2013 +0800
@@ -648,7 +648,7 @@
      CALL equivalentbodyforce(moment,in%dx1,in%dx2,in%dx3,in%sx1,in%sx2,in%sx3/2,v1,v2,v3,t1,t2,t3)
 
      ! add time-dependent surface loads
-     CALL traction(in%mu,in%events(e),in%sx1,in%sx2,in%dx1,in%dx2,t,Dt,t3,rate=.true.)
+     CALL traction(in%mu,in%events(e),in%sx1,in%sx2,in%dx1,in%dx2,t,Dt,t3,rate=.TRUE.)
 
      ! export equivalent body forces
      IF (isoutput(in%skip,t,i,in%odt,oi,in%events(e)%time)) THEN



More information about the CIG-COMMITS mailing list