[cig-commits] [commit] master: fix a bug for nonlinear viscous zones. (6b5c2f2)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/relax/compare/392758611720e5bf2bede2fa63db44922d54edda...6b5c2f26e0e2cb532d632dff1e43c559bc030db4

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

commit 6b5c2f26e0e2cb532d632dff1e43c559bc030db4
Author: Sylvain Barbot <sbarbot at ntu.edu.sg>
Date:   Sun Aug 24 22:25:48 2014 -0700

    fix a bug for nonlinear viscous zones.


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

6b5c2f26e0e2cb532d632dff1e43c559bc030db4
 src/relax.f90          | 4 ++--
 src/viscoelastic3d.f90 | 0
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/relax.f90 b/src/relax.f90
index 227f74c..f9d5bd4 100644
--- a/src/relax.f90
+++ b/src/relax.f90
@@ -640,10 +640,10 @@ PROGRAM relax
      IF (ALLOCATED(in%nonlinearstruc)) THEN
         ! powerlaw viscosity
         v1=0
-        IF (0 .LT. in%nlwz) THEN
+        IF (0 .LT. in%nnlwz) THEN
            CALL viscouseigenstress(in%mu,in%nonlinearstruc, &
                 sig,in%sx1,in%sx2,in%sx3/2, &
-                in%dx1,in%dx2,in%dx3,moment,DGAMMADOT0=lineardgammadot0,GAMMA=v1)
+                in%dx1,in%dx2,in%dx3,moment,DGAMMADOT0=nonlineardgammadot0,GAMMA=v1)
         ELSE
            CALL viscouseigenstress(in%mu,in%nonlinearstruc, &
                 sig,in%sx1,in%sx2,in%sx3/2, &