Index: /issm/trunk-jpl/test/NightlyRun/test1501.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test1501.m	(revision 21504)
+++ /issm/trunk-jpl/test/NightlyRun/test1501.m	(revision 21505)
@@ -108,5 +108,5 @@
 	for t=starttime:endtime
 		thickness = [thickness (md.results.TransientSolution(t).Thickness)];
-		volume = [volume mean(md.results.TransientSolution(t).Thickness.value,2).*areas];
+		volume = [volume mean(md.results.TransientSolution(t).Thickness(md.mesh.elements),2).*areas];
 		massbal = [massbal (md.results.TransientSolution(t).SmbMassBalance)];
 		velocity = [velocity (md.results.TransientSolution(t).Vel)];
@@ -132,5 +132,5 @@
 		[data datatype]=processdata(md,results(i).(field),options);
 
-		titlestring=[field ' at time ' num2str(results(i).time/md.constants.yts) ' year'];
+		titlestring=[field ' at time ' num2str(results(i).time) ' year'];
 		plot_unit(x,y,z,elements,data,is2d,isplanet,datatype,options)
 		options=changefieldvalue(options,'title',titlestring);
@@ -147,5 +147,5 @@
 		[data datatype]=processdata(md,results(i).(field),options);
 
-		titlestring=[field ' at time ' num2str(results(i).time/md.constants.yts) ' year'];
+		titlestring=[field ' at time ' num2str(results(i).time) ' year'];
 		plot_unit(x,y,z,elements,data,is2d,isplanet,datatype,options)
 		options=changefieldvalue(options,'title',titlestring);
@@ -196,5 +196,5 @@
 		else
 			frame=getframe(gcf);
-			images(:,:,1,count) = rgb2ind(frame.cdata,map,'nodither');
+			images(:,:,1,count) = rgb2ind(frame.cdata,[map; map],'nodither');
 		end