Index: /issm/trunk-jpl/test/NightlyRun/test273.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test273.m	(revision 19382)
+++ /issm/trunk-jpl/test/NightlyRun/test273.m	(revision 19382)
@@ -0,0 +1,27 @@
+%Test Name: SquareShelfStressSSA2dDamageUpdate
+md=triangle(model(),'../Exp/Square.exp',150000.);
+md=setmask(md,'all','');
+md.materials=matdamageice();
+md=parameterize(md,'../Par/SquareShelf.par');
+md.damage.isdamage=1;
+md.damage.D=0.*ones(md.mesh.numberofvertices,1);
+md.damage.spcdamage=NaN*ones(md.mesh.numberofvertices,1);
+md=setflowequation(md,'SSA','all');
+md.cluster=generic('name',oshostname(),'np',3);
+
+md.stressbalance.requested_outputs={'default','NewDamage'}
+md.damage.stress_threshold=1.3e5;
+md.damage.kappa=2.8;
+
+md=solve(md,StressbalanceSolutionEnum());
+
+%Fields and tolerances to track changes
+field_names     ={'Vx','Vy','Vel','Pressure','NewDamage'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	(md.results.StressbalanceSolution.Vx),...
+	(md.results.StressbalanceSolution.Vy),...
+	(md.results.StressbalanceSolution.Vel),...
+	(md.results.StressbalanceSolution.Pressure),...
+	(md.results.StressbalanceSolution.NewDamage),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test273.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test273.py	(revision 19382)
+++ /issm/trunk-jpl/test/NightlyRun/test273.py	(revision 19382)
@@ -0,0 +1,39 @@
+#Test Name: SquareShelfStressSSA2dDamageUpdate
+import numpy
+from model import model
+from triangle import triangle
+from setmask import setmask
+from parameterize import parameterize
+from setflowequation import setflowequation
+from EnumDefinitions import StressbalanceSolutionEnum
+from solve import solve
+from matdamageice import matdamageice
+from generic import generic
+from socket import gethostname as oshostname
+
+md=triangle(model(),'../Exp/Square.exp',150000.)
+md=setmask(md,'all','')
+md.materials=matdamageice()
+md=parameterize(md,'../Par/SquareShelf.py')
+md.damage.isdamage=1
+md.damage.D=0.*numpy.ones(md.mesh.numberofvertices)
+md.damage.spcdamage=numpy.nan*numpy.ones(md.mesh.numberofvertices)
+md=setflowequation(md,'SSA','all')
+md.cluster=generic('name',oshostname(),'np',3)
+
+md.stressbalance.requested_outputs=['default','NewDamage']
+md.damage.stress_threshold=1.3e5
+md.damage.kappa=2.8
+
+md=solve(md,StressbalanceSolutionEnum())
+
+#Fields and tolerances to track changes
+field_names     =['Vx','Vy','Vel','Pressure','NewDamage']
+field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13]
+field_values=[\
+	md.results.StressbalanceSolution.Vx,\
+	md.results.StressbalanceSolution.Vy,\
+	md.results.StressbalanceSolution.Vel,\
+	md.results.StressbalanceSolution.Pressure,\
+	md.results.StressbalanceSolution.NewDamage,\
+	]