Index: /issm/trunk/test/NightlyRun/test223.m
===================================================================
--- /issm/trunk/test/NightlyRun/test223.m	(revision 5021)
+++ /issm/trunk/test/NightlyRun/test223.m	(revision 5021)
@@ -0,0 +1,20 @@
+md=mesh(model,'../Exp/Square.exp',200000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=setelementstype(md,'macayeal','all');
+
+%control parameters
+md.control_analysis=1;
+md.control_type='rheology_B';
+md.cm_min=10^6;
+md.cm_max=2*10^9;
+md.nsteps=2;
+md.fit=0*ones(md.nsteps,1);
+md.weights=ones(md.numberofgrids,1);
+md.optscal=10^7*ones(md.nsteps,1);
+md.maxiter=2*ones(md.nsteps,1);
+md.cm_jump=0.3*ones(md.nsteps,1);
+md.vx_obs=md.vx; md.vy_obs=md.vy;
+
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test223_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test223_nightly.m	(revision 5021)
+++ /issm/trunk/test/NightlyRun/test223_nightly.m	(revision 5021)
@@ -0,0 +1,11 @@
+field_names     ={'Gradient' 'Misfits' 'RheologyB2d' 'Pressure' 'Vel' 'Vx' 'Vy'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.DiagnosticSolution.Gradient),...
+	md.results.DiagnosticSolution.J,...
+	PatchToVec(md.results.DiagnosticSolution.RheologyB2d),...
+	PatchToVec(md.results.DiagnosticSolution.Pressure),...
+	PatchToVec(md.results.DiagnosticSolution.Vel),...
+	PatchToVec(md.results.DiagnosticSolution.Vx),...
+	PatchToVec(md.results.DiagnosticSolution.Vy)
+};