Index: /issm/trunk-jpl/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 17392)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 17393)
@@ -71,4 +71,5 @@
 	case 272, name='SquareShelfCMZSSA2dDamage';
 	case 274, name='SquareShelfStressSSA2dDamageRift';
+	case 275, name='SquareShelfDamageEvolutionSSA2d';
 	case 280, name='SquareShelfStressSSA2dHigherOrder';
 	case 285, name='SquareShelfStressHOHigherOrder';
Index: /issm/trunk-jpl/test/NightlyRun/IdToName.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 17392)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 17393)
@@ -72,4 +72,5 @@
 	272  : 'SquareShelfCMZSSA2dDamage',
 	274  : 'SquareShelfStressSSA2dDamageRift',
+	275  : 'SquareShelfDamageEvolutionSSA2d',
 	280  : 'SquareShelfStressSSA2dHigherOrder',
 	285  : 'SquareShelfStressHOHigherOrder',
Index: /issm/trunk-jpl/test/NightlyRun/runme.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/runme.py	(revision 17392)
+++ /issm/trunk-jpl/test/NightlyRun/runme.py	(revision 17393)
@@ -141,7 +141,8 @@
 			archive_name='Archive'+str(id)
 			if strcmpi(procedure,'update'):
-
-				os.remove(os.path.join('..','Archives',archive_name+'.nc'))
-				f = netCDF4.Dataset(os.path.join('..','Archives',archive_name+'.nc'),'w',format='NETCDF3_CLASSIC')
+				archive=os.path.join('..','Archives',archive_name+'.nc')
+				if os.path.isfile(archive):
+					os.remove(archive)
+				f = netCDF4.Dataset(archive,'w',format='NETCDF3_CLASSIC')
 				for k,fieldname in enumerate(field_names):
 					field=numpy.array(field_values[k],dtype=float)
Index: /issm/trunk-jpl/test/NightlyRun/test275.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test275.m	(revision 17393)
+++ /issm/trunk-jpl/test/NightlyRun/test275.m	(revision 17393)
@@ -0,0 +1,26 @@
+md=triangle(model(),'../Exp/Square.exp',150000.);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md.damage.D=0.1*ones(md.mesh.numberofvertices,1);
+md.damage.spcdamage=NaN*ones(md.mesh.numberofvertices,1);
+md.damage.law='pralong';
+
+md.damage.c1=1.e-11;
+md.damage.c2=0.4;
+md.damage.c3=1.e-3;
+md.damage.healing=0.4;
+md.damage.stress_threshold=1e5;
+md.damage.stabilization=1;
+
+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={1e-13,1e-13}
+field_values={...
+		(md.results.DamageEvolutionSolution.DamageD),...
+		(md.results.DamageEvolutionSolution.DamageF),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test275.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test275.py	(revision 17393)
+++ /issm/trunk-jpl/test/NightlyRun/test275.py	(revision 17393)
@@ -0,0 +1,38 @@
+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=0.1*npy.ones((md.mesh.numberofvertices,1))
+md.damage.spcdamage=npy.float('nan')*npy.ones((md.mesh.numberofvertices,1))
+md.damage.law='pralong'
+
+md.damage.c1=1.e-11
+md.damage.c2=0.4
+md.damage.c3=1.e-3
+md.damage.healing=0.4
+md.damage.stress_threshold=1.e5
+md.damage.stabilization=1
+
+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,\
+		]