[cig-commits] commit: fix treatment of out-of-phase periodic loads

[Mercurial]

changeset:   212:45d392635687
tag:         tip
user:        Sylvain Barbot <sbarbot at ntu.edu.sg>
date:        Sat Aug 10 23:17:04 2013 +0800
files:       src/elastic3d.f90
description:
fix treatment of out-of-phase periodic loads


diff -r eccc38abafa0 -r 45d392635687 src/elastic3d.f90
--- a/src/elastic3d.f90	Sat Aug 10 21:59:03 2013 +0800
+++ b/src/elastic3d.f90	Sat Aug 10 23:17:04 2013 +0800
@@ -2074,6 +2074,9 @@
        x=e%l(i)%x
        y=e%l(i)%y
 
+       period=e%l(i)%period
+       phi=e%l(i)%phase
+
        L=e%l(i)%length
        W=e%l(i)%width
 
@@ -2095,19 +2098,17 @@
                  omega((x2-y)/W,beta)* &
                  e%l(i)%slip
 
-             period=e%l(i)%period
-
              IF (israte) THEN
+                ! surface tractions rate
                 IF (0 .NE. period) THEN
-                   ! 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))/Dt)
                 END IF
              ELSE
-                IF (e%l(i)%period .LE. 0) THEN
-                   ! surface tractions
+                ! surface tractions
+                IF (period .EQ. 0) THEN
                    t3(i1,i2)=REAL(t3(i1,i2)-amp)
+                ELSE
+                   t3(i1,i2)=REAL(t3(i1,i2)-amp*sin(2*pi*t/period+phi))
                 END IF
              END IF
           END DO