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

[willic3 at geodynamics.org]

Author: willic3
Date: 2007-01-12 13:20:20 -0800 (Fri, 12 Jan 2007)
New Revision: 5780

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/matmod_def.f
Log:
Replaced original power-law model (model 6) with new implementation and
moved original implementation to material model 9.



Modified: short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/matmod_def.f
===================================================================
--- short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/matmod_def.f	2007-01-12 21:19:06 UTC (rev 5779)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/matmod_def.f	2007-01-12 21:20:20 UTC (rev 5780)
@@ -143,9 +143,18 @@
       infmatmod(5,8) = izero
       infmatmod(6,8) = 6
 c
+c...  Definition for power-law Maxwell viscoelastic material (ESF)
+c
+      infmatmod(1,9) = ione
+      infmatmod(2,9) = 18
+      infmatmod(3,9) = ifive
+      infmatmod(4,9) = ione
+      infmatmod(5,9) = ione
+      infmatmod(6,9) = 6
+c
 c...  Dummy definitions for remaining materials
 c
-      do i=9,nmatmodmax
+      do i=10,nmatmodmax
         do j=2,6
           infmatmod(j,i)=infmatmod(j,5)
         end do