Index: /issm/trunk-jpl/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 13162)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 13163)
@@ -186,4 +186,5 @@
 	case 1601, name='SquareShelfM2dRotation';
 	case 1602, name='SquareSheetShelfP3dRotation';
+	case 3001, name='SquareShelfConstrainedDiagM2dAdolc';
 	otherwise, name='N/A';
 end
Index: /issm/trunk-jpl/test/NightlyRun/IdToName.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 13162)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 13163)
@@ -186,5 +186,6 @@
 		1502 : 'SquareShelfTranSawTooth3d',
 		1601 : 'SquareShelfM2dRotation',
-		1602 : 'SquareSheetShelfP3dRotation',
+		1602 : 'SquareSheetShelfP3dRotation', 
+		3001 : 'SquareShelfConstrainedDiagM2dAdolc',
 	}
 
Index: /issm/trunk-jpl/test/NightlyRun/test3001.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3001.m	(revision 13163)
+++ /issm/trunk-jpl/test/NightlyRun/test3001.m	(revision 13163)
@@ -0,0 +1,23 @@
+md=triangle(model,'../Exp/Square.exp',50000);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setflowequation(md,'macayeal','all');
+md.cluster=generic('name',oshostname(),'np',3);
+md.diagnostic.requested_outputs=StressTensorEnum;
+md.autodiff.isautodiff=false;
+md=solve(md,DiagnosticSolutionEnum);
+
+%Fields and tolerances to track changes
+field_names     ={'Vx','Vy','Vel','Pressure',...
+	'StressTensorxx','StressTensoryy','StressTensorxy'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,...
+	1e-13,1e-13,1e-13};
+field_values={...
+	(md.results.DiagnosticSolution.Vx),...
+	(md.results.DiagnosticSolution.Vy),...
+	(md.results.DiagnosticSolution.Vel),...
+	(md.results.DiagnosticSolution.Pressure),...
+	(md.results.DiagnosticSolution.StressTensorxx),...
+	(md.results.DiagnosticSolution.StressTensoryy),...
+	(md.results.DiagnosticSolution.StressTensorxy),...
+	};