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

[willic3 at geodynamics.org]

Author: willic3
Date: 2006-10-25 07:12:38 -0700 (Wed, 25 Oct 2006)
New Revision: 5090

Modified:
   short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/matmod_def.f
Log:
Added power-law Maxwell definitions.



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	2006-10-25 14:11:57 UTC (rev 5089)
+++ short/3D/PyLith/branches/pylith-0.8/pylith3d/libpylith3d/matmod_def.f	2006-10-25 14:12:38 UTC (rev 5090)
@@ -107,7 +107,7 @@
       infmatmod(6,3) =6
       infmatmod(6,4) =6
 c
-c...  Definition for Maxwell viscoelastic material
+c...  Definition for linear Maxwell viscoelastic material
 c
       infmatmod(1,5) = ione
       infmatmod(2,5) = 18
@@ -116,9 +116,18 @@
       infmatmod(5,5) = izero
       infmatmod(6,5) = 6
 c
+c...  Definition for power-law Maxwell viscoelastic material
+c
+      infmatmod(1,6) = ione
+      infmatmod(2,6) = 18
+      infmatmod(3,6) = ifive
+      infmatmod(4,6) = ione
+      infmatmod(5,6) = ione
+      infmatmod(6,6) = 6
+c
 c...  Dummy definitions for remaining materials
 c
-      do i=6,nmatmodmax
+      do i=7,nmatmodmax
         do j=2,6
           infmatmod(j,i)=infmatmod(j,5)
         end do