Index: /issm/trunk-jpl/test/NightlyRun/test4005.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test4005.m	(revision 27600)
+++ /issm/trunk-jpl/test/NightlyRun/test4005.m	(revision 27601)
@@ -666,12 +666,9 @@
 		 system(['cp run/shelfice_dmdt.bin run/shelfice_dmdt_' int2str(t) '.bin'])
 
-
 		 % {{{ system call to run MITgcm
 		 cd run
-
-                 newline = [' niter0 = ' num2str(t*y2s/MITgcmDeltaT)];
+       newline = [' niter0 = ' num2str(t*y2s/MITgcmDeltaT)];
 		 command=['!sed "s/.*niter0.*/' newline '/" data > data.temp; mv data.temp data'];
 		 eval(command)
-
 
        ds=round(endtime/MITgcmDeltaT);
@@ -696,5 +693,4 @@
 
        melt=readbin('./run/melt.data',[Nx,Ny]);
-
        md.basalforcings.floatingice_melting_rate=-melt(:)*y2s/rho_ice;
 		 md=solve(md,'Transient');
@@ -720,43 +716,43 @@
 % }}}
 
-%Fields and tolerances to track changes
-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]);
-field_names     ={'Base1','Melting1','Vx2','Vy2','Thickness2','Base2','MaskOceanLevelset2','FloatingiceMeltingRate2',...
-    'Melting2','Vx3','Vy3','Thickness3','Base3','MaskOceanLevelset3','FloatingiceMeltingRate3',...
-    'Melting3','Vx4','Vy4','Thickness4','Base4','MaskOceanLevelset4','FloatingiceMeltingRate4','Melting4'};
-field_tolerances={2e-13,1e-13,...
-    8e-06,7e-06,2e-07,2e-08,3e-08,1e-13,1e-13,...
-    8e-06,7e-06,4e-07,3e-08,5e-08,1e-13,6e-08,...
-    8e-06,7e-06,5e-07,4e-08,8e-08,6e-08,3e-07};
-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).MaskOceanLevelset),...
-    (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).MaskOceanLevelset),...
-    (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).MaskOceanLevelset),...
-    (md.results(4).TransientSolution(end).FloatingiceMeltingRate),...
-    (melting_rate_4(:)),...
-    };
+%%Fields and tolerances to track changes
+%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]);
+%field_names     ={'Base1','Melting1','Vx2','Vy2','Thickness2','Base2','MaskOceanLevelset2','FloatingiceMeltingRate2',...
+%    'Melting2','Vx3','Vy3','Thickness3','Base3','MaskOceanLevelset3','FloatingiceMeltingRate3',...
+%    'Melting3','Vx4','Vy4','Thickness4','Base4','MaskOceanLevelset4','FloatingiceMeltingRate4','Melting4'};
+%field_tolerances={2e-13,1e-13,...
+%    8e-06,7e-06,2e-07,2e-08,3e-08,1e-13,1e-13,...
+%    8e-06,7e-06,4e-07,3e-08,5e-08,1e-13,6e-08,...
+%    8e-06,7e-06,5e-07,4e-08,8e-08,6e-08,3e-07};
+%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).MaskOceanLevelset),...
+%    (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).MaskOceanLevelset),...
+%    (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).MaskOceanLevelset),...
+%    (md.results(4).TransientSolution(end).FloatingiceMeltingRate),...
+%    (melting_rate_4(:)),...
+%    };