Changeset 17818

Timestamp:

04/22/14 19:38:22 (11 years ago)

Author:

cborstad

Message:

CHG: adding a print statement to remind user that damage-dependent properties are being calculated

File:

• issm/trunk-jpl/src/m/mech/mechanicalproperties.py (modified) (1 diff)

issm/trunk-jpl/src/m/mech/mechanicalproperties.py

@@ -76 +76 @@
                 nu[location]=10^18
         elif 'matdamageice' in md.materials.__module__:
-                Zinv=npy.dot(md.materials.rheology_Z[index-1],summation/3.).reshape(-1,)
+                print 'computing damage-dependent properties!'
+                Zinv=npy.dot(1-md.damage.D[index-1],summation/3.).reshape(-1,)
                 location=npy.nonzero(second_inv)
                 nu[location]=Zinv[location]*B_bar[location]/npy.power(second_inv[location],power[location])