Index: /issm/branches/trunk-jpl-damage/src/m/model/damage.m
===================================================================
--- /issm/branches/trunk-jpl-damage/src/m/model/damage.m	(revision 12901)
+++ /issm/branches/trunk-jpl-damage/src/m/model/damage.m	(revision 12901)
@@ -0,0 +1,12 @@
+function md=damage(md)
+%DAMAGE - calculated damage D from inverse of damage effect tensor Z^-1
+%
+%   Z^-1 = (1 - D) for isotropic scalar damage 
+%	 Z^-1 is MaterialsRheologyZbar
+%
+%   Usage:
+%      md=damage(md)
+
+damage=struct('D',[]);
+damage.D=1-md.results.DiagnosticSolution.MaterialsRheologyZbar;
+md.results.damage=damage;