Index: /issm/trunk-jpl/test/NightlyRun/test4003.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test4003.m	(revision 24127)
+++ /issm/trunk-jpl/test/NightlyRun/test4003.m	(revision 24128)
@@ -371,4 +371,5 @@
 		  md.cluster.np=2;
 		  md.cluster.executionpath=[pwd '/run'];
+		  md.transient.requested_outputs={'default','MaskGroundediceLevelset'};
 
 		  md=solveiceocean(md,'Transient','runtimename',false);
@@ -667,2 +668,43 @@
 end
 % }}}
+
+%Fields and tolerances to track changes
+fnm=['run/SHICE_fwFlux.0000011400.data'];
+melting_rate_1=readbin(fnm,[Nx Ny]);
+fnm=['run/SHICE_fwFlux.0000011450.data'];
+melting_rate_2=readbin(fnm,[Nx Ny]);
+fnm=['run/SHICE_fwFlux.0000011500.data'];
+melting_rate_3=readbin(fnm,[Nx Ny]);
+fnm=['run/SHICE_fwFlux.0000011550.data'];
+melting_rate_4=readbin(fnm,[Nx Ny]);
+field_names     ={'Base1','Melting1','Vx2','Vy2','Thickness2','Base2','MaskGroundediceLevelset2','FloatingiceMeltingRate2',...
+	'Melting2','Vx3','Vy3','Thickness3','Base3','MaskGroundediceLevelset3','FloatingiceMeltingRate3',...
+	'Melting3','Vx4','Vy4','Thickness4','Base4','MaskGroundediceLevelset4','FloatingiceMeltingRate4','Melting4'};
+field_tolerances={2e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-12,...
+	1e-13, 1e-12, 1e-12, 1e-13, 1e-13, 1e-13, 1e-12,...
+	1e-13, 1e-12, 1e-12, 1e-13, 1e-13, 1e-13, 1e-12, 1e-13 };
+field_values={...
+	(md.results.TransientSolution(1).Base),...
+	(melting_rate_1(:)),...
+	(md.results.TransientSolution(2).Vx),...
+	(md.results.TransientSolution(2).Vy),...
+	(md.results.TransientSolution(2).Thickness),...
+	(md.results.TransientSolution(2).Base),...
+	(md.results.TransientSolution(2).MaskGroundediceLevelset),...
+	(md.results.TransientSolution(2).BasalforcingsFloatingiceMeltingRate),...
+	(melting_rate_2(:)),...
+	(md.results.TransientSolution(3).Vx),...
+	(md.results.TransientSolution(3).Vy),...
+	(md.results.TransientSolution(3).Thickness),...
+	(md.results.TransientSolution(3).Base),...
+	(md.results.TransientSolution(3).MaskGroundediceLevelset),...
+	(md.results.TransientSolution(3).BasalforcingsFloatingiceMeltingRate),...
+	(melting_rate_3(:)),...
+	(md.results.TransientSolution(4).Vx),...
+	(md.results.TransientSolution(4).Vy),...
+	(md.results.TransientSolution(4).Thickness),...
+	(md.results.TransientSolution(4).Base),...
+	(md.results.TransientSolution(4).MaskGroundediceLevelset),...
+	(md.results.TransientSolution(4).BasalforcingsFloatingiceMeltingRate),...
+	(melting_rate_4(:)),...
+	};