[cig-commits] [commit] master: add a test for zero length of width in subroutine traction. cycle through for point tractions. (80e49e8)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/relax

On branch  : master
Link       : https://github.com/geodynamics/relax/compare/293b281689facd77a00ce8ffa744eb7e8e21339e...17d2ea614ff358822f81bb94d3774925a60710fb

>---------------------------------------------------------------

commit 80e49e8836bf09a29f68cfeb3952352b8e7e729b
Author: Sylvain Barbot <sbarbot at ntu.edu.sg>
Date:   Sat Jan 24 22:00:41 2015 +0800

    add a test for zero length of width in subroutine traction. cycle through for point tractions.


>---------------------------------------------------------------

80e49e8836bf09a29f68cfeb3952352b8e7e729b
 src/elastic3d.f90 | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/src/elastic3d.f90 b/src/elastic3d.f90
index e573adc..454d0c7 100644
--- a/src/elastic3d.f90
+++ b/src/elastic3d.f90
@@ -2093,6 +2093,9 @@ CONTAINS
        Lp=L*(1._8+2._8*beta)/2._8
        Wp=W*(1._8+2._8*beta)/2._8
 
+       IF (0 .EQ. Lp) CYCLE
+       IF (0 .EQ. Wp) CYCLE
+
        i3=1
        DO i2=1,sx2
           DO i1=1,sx1
@@ -3088,7 +3091,7 @@ CONTAINS
   END SUBROUTINE fieldrep
 
   !-----------------------------------------------------------------
-  ! subroutine SliveAdd
+  ! subroutine SliceAdd
   ! computes in place the sum of two scalar fields
   !
   !   u = c1 * u + c2 * v