[cig-commits] commit: fix calculation of non-periodic loads

[Mercurial]

changeset:   203:d7bc9499126c
tag:         tip
user:        Sylvain Barbot <sbarbot at ntu.edu.sg>
date:        Wed May 15 11:46:11 2013 +0800
files:       src/elastic3d.f90
description:
fix calculation of non-periodic loads


diff -r bb0f18cb5f71 -r d7bc9499126c src/elastic3d.f90
--- a/src/elastic3d.f90	Wed May 15 11:45:48 2013 +0800
+++ b/src/elastic3d.f90	Wed May 15 11:46:11 2013 +0800
@@ -2095,10 +2095,11 @@ CONTAINS
                  omega((x2-y)/W,beta)* &
                  mu*e%l(i)%slip
 
+             period=e%l(i)%period
+
              IF (israte) THEN
                 IF (0 .NE. period) THEN
                    ! surface tractions rate
-                   period=e%l(i)%period
                    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)))