[cig-commits] commit: fix the amplitude unit of the periodic loads.

[Mercurial]

changeset:   204:8afab31fef81
tag:         tip
user:        Sylvain Barbot <sbarbot at ntu.edu.sg>
date:        Wed May 15 12:05:16 2013 +0800
files:       src/elastic3d.f90
description:
fix the amplitude unit of the periodic loads.


diff -r d7bc9499126c -r 8afab31fef81 src/elastic3d.f90
--- a/src/elastic3d.f90	Wed May 15 11:46:11 2013 +0800
+++ b/src/elastic3d.f90	Wed May 15 12:05:16 2013 +0800
@@ -2102,7 +2102,7 @@ CONTAINS
                    ! surface tractions rate
                    phi=e%l(i)%phase
 
-                   t3(i1,i2)=REAL(t3(i1,i2)-amp*(sin(2*pi*(t+Dt)/period+phi)-sin(2*pi*t/period+phi)))
+                   t3(i1,i2)=REAL(t3(i1,i2)-amp*(sin(2*pi*(t+Dt)/period+phi)-sin(2*pi*t/period+phi))/Dt)
                 END IF
              ELSE
                 IF (e%l(i)%period .LE. 0) THEN