Index: /issm/trunk-jpl/test/NightlyRun/test271.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test271.m	(revision 17399)
+++ /issm/trunk-jpl/test/NightlyRun/test271.m	(revision 17400)
@@ -2,5 +2,5 @@
 md=setmask(md,'all','');
 md=parameterize(md,'../Par/SquareShelf.par');
-md.damage.D=0.5*ones(md.mesh.numberofvertices,1);
+md.damage.D=zeros(md.mesh.numberofvertices,1);
 md.damage.spcdamage=NaN*ones(md.mesh.numberofvertices,1);
 md.damage.law='pralong';
@@ -10,7 +10,8 @@
 boundary(md.mesh.segments(:,1:2))=1;
 pos=find(~md.mask.ice_levelset==0 & boundary==1);
-md.damage.spcdamage(pos)=0;
+md.damage.spcdamage(pos)=0.5;
+md.damage.stabilization=1;
 
-md.verbose=verbose('solution',true);
+md.damage.requested_outputs={'default','DamageF'}
 
 md=setflowequation(md,'SSA','all');
@@ -19,5 +20,8 @@
 
 %Fields and tolerances to track changes
-field_names     ={'D'};
-field_tolerances={1e-13};
-field_values={md.results.DamageEvolutionSolution.DamageD};
+field_names={'D','F'}
+field_tolerances={1e-13,1e-13};
+field_values={...
+		(md.results.DamageEvolutionSolution.DamageD),...
+		(md.results.DamageEvolutionSolution.DamageF),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test271.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test271.py	(revision 17400)
+++ /issm/trunk-jpl/test/NightlyRun/test271.py	(revision 17400)
@@ -0,0 +1,35 @@
+import numpy as npy
+from triangle import triangle
+from model import model
+from setmask import setmask
+from parameterize import parameterize
+from verbose import verbose
+from setflowequation import setflowequation
+from generic import generic
+from socket import gethostname as oshostname
+from solve import solve
+from EnumDefinitions import DamageEvolutionSolutionEnum
+
+md=triangle(model(),'../Exp/Square.exp',150000.)
+md=setmask(md,'all','')
+md=parameterize(md,'../Par/SquareShelf.py')
+md.damage.D=npy.zeros((md.mesh.numberofvertices,1))
+md.damage.spcdamage=npy.float('nan')*npy.ones((md.mesh.numberofvertices,1))
+md.damage.law='pralong'
+md.damage.stabilization=1;
+
+pos=npy.nonzero(npy.logical_and(md.mask.ice_levelset[:,0]!=0,md.mesh.vertexonboundary==1))
+md.damage.spcdamage[pos]=0.5;
+
+md.damage.requested_outputs=['default','DamageF']
+
+md=setflowequation(md,'SSA','all')
+md.cluster=generic('name',oshostname(),'np',3)
+md=solve(md,DamageEvolutionSolutionEnum())
+
+field_names=['D','F']
+field_tolerances=[1.e-13,1.e-13]
+field_values=[\
+		md.results.DamageEvolutionSolution.DamageD,\
+		md.results.DamageEvolutionSolution.DamageF,\
+		]