fix a bug for nonlinear viscous zones.

geodynamics · Aug 25, 2014 · 6b5c2f2 · 6b5c2f2
1 parent abc4c6a
commit 6b5c2f2
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git 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, &

diff --git a/src/viscoelastic3d.f90 b/src/viscoelastic3d.f90
@@ -272,7 +272,7 @@ SUBROUTINE builddgammadot0(sx1,sx2,sx3,dx1,dx2,dx3,beta, &
 
     INTEGER :: i1,i2,i3
     REAL*8 :: x1,x2,x3,dum
-    
+
 !$omp parallel do private(i1,i2,x1,x2,x3,dum)
     DO i3=1,sx3
        x3=DBLE(i3-1)*dx3