[cig-commits] r12826 - short/3D/PyLith/trunk/libsrc/materials

[willic3 at geodynamics.org]

Author: willic3
Date: 2008-09-06 13:58:26 -0700 (Sat, 06 Sep 2008)
New Revision: 12826

Modified:
   short/3D/PyLith/trunk/libsrc/materials/GenMaxwellIsotropic3D.cc
Log:
Maxwell time for an element should be set to a large value (rather than
zero) if the element is not being used (shear_ratio = 0 for the element).



Modified: short/3D/PyLith/trunk/libsrc/materials/GenMaxwellIsotropic3D.cc
===================================================================
--- short/3D/PyLith/trunk/libsrc/materials/GenMaxwellIsotropic3D.cc	2008-09-06 14:18:44 UTC (rev 12825)
+++ short/3D/PyLith/trunk/libsrc/materials/GenMaxwellIsotropic3D.cc	2008-09-06 20:58:26 UTC (rev 12826)
@@ -161,7 +161,7 @@
     double muRatio = dbValues[_GenMaxwellIsotropic3D::didShearRatio1 + model];
     double viscosity = dbValues[_GenMaxwellIsotropic3D::didViscosity1 + model];
     double muFac = muRatio*mu;
-    double maxwellTime = 0.0;
+    double maxwellTime = 1.0e30;
     if (muFac > 0.0) maxwellTime = viscosity / muFac;
     propValues[_GenMaxwellIsotropic3D::pidShearRatio + model] = muRatio;
     propValues[_GenMaxwellIsotropic3D::pidMaxwellTime + model] = maxwellTime;
@@ -514,7 +514,7 @@
   double shearRatio = 0.0;
   for (int model = 0; model < numMaxwellModels; ++model) {
     shearRatio = properties[_GenMaxwellIsotropic3D::pidShearRatio + model];
-    double maxwellTime = 0.0;
+    double maxwellTime = 1.0e30;
     visFrac += shearRatio;
     if (shearRatio != 0.0) {
       maxwellTime = properties[_GenMaxwellIsotropic3D::pidMaxwellTime + model];