Index: /issm/trunk-jpl/test/NightlyRun/test4002.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test4002.m	(revision 22725)
+++ /issm/trunk-jpl/test/NightlyRun/test4002.m	(revision 22726)
@@ -846,5 +846,6 @@
 	loaddata(org,'Bathymetry');
 	loaddata(org,'IceSheetGeometry');
-        endtime = round(MITgcmDeltaT * floor(final_time*y2s/MITgcmDeltaT));
+		endtime = round(MITgcmDeltaT * floor(final_time*y2s/MITgcmDeltaT));
+		outputtime = round(MITgcmDeltaT * floor(time_step*y2s/MITgcmDeltaT));
 
         % {{{ prepare MITgcm 
@@ -967,5 +968,5 @@
             end
             if strcmpi(tline,' taveFreq=2592000.,'),
-                fprintf(fido,'%s%i%s\n',' taveFreq= ',endtime,',');
+                fprintf(fido,'%s%i%s\n',' taveFreq= ',outputtime,',');
                 continue;
             end
@@ -1002,17 +1003,12 @@
 
 %Fields and tolerances to track changes
-fnm=['run/SHICE_fwFluxtave.0000000146.data'];
+fnm=['run/SHICE_fwFluxtave.0000000006.data'];
 melting_rate_1=readbin(fnm,[Nx Ny]);
-melting_rate_2=melting_rate_1;
-melting_rate_3=melting_rate_1;
-melting_rate_4=melting_rate_1;
-%fnm=['run/SHICE_fwFluxtave_0000004380.data'];
-%melting_rate_1=readbin(fnm,[Nx Ny]);
-%fnm=['run/SHICE_fwFluxtave_0000008760.data'];
-%melting_rate_2=readbin(fnm,[Nx Ny]);
-%fnm=['run/SHICE_fwFluxtave_0000013140.data'];
-%melting_rate_3=readbin(fnm,[Nx Ny]);
-%fnm=['run/SHICE_fwFluxtave_0000017520.data'];
-%melting_rate_4=readbin(fnm,[Nx Ny]);
+fnm=['run/SHICE_fwFluxtave.0000000012.data'];
+melting_rate_2=readbin(fnm,[Nx Ny]);
+fnm=['run/SHICE_fwFluxtave.0000000018.data'];
+melting_rate_3=readbin(fnm,[Nx Ny]);
+fnm=['run/SHICE_fwFluxtave.0000000024.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',...
@@ -1046,27 +1042,2 @@
 	(melting_rate_4(:)),...
 	};
-%field_values={...
-%	(md.results(1).TransientSolution(end).Base),...
-%	(melting_rate_1(:)),...
-%	(md.results(2).TransientSolution(end).Vx),...
-%	(md.results(2).TransientSolution(end).Vy),...
-%	(md.results(2).TransientSolution(end).Thickness),...
-%	(md.results(2).TransientSolution(end).Base),...
-%	(md.results(2).TransientSolution(end).MaskGroundediceLevelset),...
-%	(md.results(2).TransientSolution(end).FloatingiceMeltingRate),...
-%	(melting_rate_2(:)),...
-%	(md.results(3).TransientSolution(end).Vx),...
-%	(md.results(3).TransientSolution(end).Vy),...
-%	(md.results(3).TransientSolution(end).Thickness),...
-%	(md.results(3).TransientSolution(end).Base),...
-%	(md.results(3).TransientSolution(end).MaskGroundediceLevelset),...
-%	(md.results(3).TransientSolution(end).FloatingiceMeltingRate),...
-%	(melting_rate_3(:)),...
-%	(md.results(4).TransientSolution(end).Vx),...
-%	(md.results(4).TransientSolution(end).Vy),...
-%	(md.results(4).TransientSolution(end).Thickness),...
-%	(md.results(4).TransientSolution(end).Base),...
-%	(md.results(4).TransientSolution(end).MaskGroundediceLevelset),...
-%	(md.results(4).TransientSolution(end).FloatingiceMeltingRate),...
-%	(melting_rate_4(:)),...
-%	};