[cig-commits] [commit] master: no smoothing for the ductile zones (beta=0). (17d2ea6)

[cig_noreply at geodynamics.org]

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

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

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

commit 17d2ea614ff358822f81bb94d3774925a60710fb
Author: Sylvain Barbot <sbarbot at ntu.edu.sg>
Date:   Sun Jan 25 18:16:26 2015 +0800

    no smoothing for the ductile zones (beta=0).


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

17d2ea614ff358822f81bb94d3774925a60710fb
 src/relax.f90 | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/relax.f90 b/src/relax.f90
index 4b80fc7..05dd762 100644
--- a/src/relax.f90
+++ b/src/relax.f90
@@ -468,7 +468,7 @@ PROGRAM relax
      IF (0 .LT. in%nlwz) THEN
         ALLOCATE(lineardgammadot0(in%sx1,in%sx2,in%sx3/2),STAT=iostatus)
         IF (iostatus.GT.0) STOP "could not allocate lineardgammadot0"
-        CALL builddgammadot0(in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3,in%beta, &
+        CALL builddgammadot0(in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3,0._8, &
                              in%nlwz,in%linearweakzone,lineardgammadot0)
      END IF
   END IF
@@ -483,7 +483,7 @@ PROGRAM relax
      IF (0 .LT. in%nnlwz) THEN
         ALLOCATE(nonlineardgammadot0(in%sx1,in%sx2,in%sx3/2),STAT=iostatus)
         IF (iostatus.GT.0) STOP "could not allocate nonlineardgammadot0"
-        CALL builddgammadot0(in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3,in%beta, &
+        CALL builddgammadot0(in%sx1,in%sx2,in%sx3/2,in%dx1,in%dx2,in%dx3,0._8, &
                              in%nnlwz,in%nonlinearweakzone,nonlineardgammadot0)
      END IF
   END IF