[cig-commits] r5341 - short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d

[willic3 at geodynamics.org]

Author: willic3
Date: 2006-11-21 18:06:37 -0800 (Tue, 21 Nov 2006)
New Revision: 5341

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/mat_7.f
Log:
Altered relaxation time to be consistent with uniaxial test problem.



Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/mat_7.f
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/mat_7.f	2006-11-22 00:48:18 UTC (rev 5340)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/mat_7.f	2006-11-22 02:06:37 UTC (rev 5341)
@@ -276,14 +276,13 @@
       do i=1,3
         rmue=prop(2*(i-1)+4)*rmu
         vise=prop(2*(i-1)+5)
-ctest        if(rmue.ne.zero) tmax=min(tmax,two*vise/rmue)
-        if(rmue.ne.zero) tmax=min(tmax,two*vise/rmue)
+        if(rmue.ne.zero) tmax=min(tmax,vise/rmue)
       end do
       return
       end
 c
 c
-      function compdq(deltp,vis,tau)
+      function compdq(deltp,tau)
 c
 c...  function to compute the viscous state variable increment delta-q.
 c
@@ -299,7 +298,7 @@
 c
 c...  subroutine arguments
 c
-      double precision compdq,deltp,vis,tau
+      double precision compdq,deltp,tau
 c
 c...  intrinsic functions
 c
@@ -411,9 +410,9 @@
         vise=prop(2*(i-1)+5)
         tau=zero
         if(rmue.ne.zero) then
-          tau=two*vise/rmut
+          tau=vise/rmut
           tmax=min(tmax,tau)
-          shfac=shfac+rmue*compdq(deltp,vise,tau)
+          shfac=shfac+rmue*compdq(deltp,tau)
         end if
       end do
 c
@@ -520,9 +519,9 @@
         vise=prop(2*(i-1)+5)
         rtime(i)=zero
         if(rmue.ne.zero) then
-          rtime(i)=two*vise/rmut
+          rtime(i)=vise/rmut
           tmax=min(tmax,rtime(i))
-          dq(i)=compdq(deltp,vise,rtime(i))
+          dq(i)=compdq(deltp,rtime(i))
         end if
       end do
 c