Index: /issm/trunk-jpl/src/m/mech/mechanicalproperties.py
===================================================================
--- /issm/trunk-jpl/src/m/mech/mechanicalproperties.py	(revision 17817)
+++ /issm/trunk-jpl/src/m/mech/mechanicalproperties.py	(revision 17818)
@@ -76,5 +76,6 @@
 		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])