Index: /issm/trunk-jpl/src/m/model/plot/plot_manager.m
===================================================================
--- /issm/trunk-jpl/src/m/model/plot/plot_manager.m	(revision 11178)
+++ /issm/trunk-jpl/src/m/model/plot/plot_manager.m	(revision 11179)
@@ -78,6 +78,8 @@
 			return;
 		case 'none',
-			plot_none(md,options,nlines,ncols,i);
-			return;
+			if ~exist(options,'overlay'),
+				plot_none(md,options,nlines,ncols,i);
+				return;
+			end
 		case 'penalties',
 			plot_penalties(md,options,subplotwidth,i);
Index: /issm/trunk-jpl/src/m/model/plot/plot_overlay.m
===================================================================
--- /issm/trunk-jpl/src/m/model/plot/plot_overlay.m	(revision 11178)
+++ /issm/trunk-jpl/src/m/model/plot/plot_overlay.m	(revision 11179)
@@ -50,7 +50,13 @@
 ncols =length(md.radaroverlay.x);
 disp('Interpolating data on grid...');
-[x_m y_m data_grid]=InterpFromMeshToGrid(elements,x/getfieldvalue(options,'unit',1),y/getfieldvalue(options,'unit',1),data,xmin,ymax,xspacing,yspacing,nlines,ncols,NaN);
+[x_m y_m data_grid]=InterpFromMeshToGrid(elements,x/getfieldvalue(options,'unit',1),y/getfieldvalue(options,'unit',1),...
+	data,xmin,ymax,xspacing,yspacing,nlines,ncols,NaN);
 
 %Process data_grid
+pos=find(isinf(data_grid));
+if ~isempty(pos),
+	disp('Warning: removing Infs from vector (probably log(0)?)');
+	data_grid(pos)=NaN;
+end
 if exist(options,'caxis'),
 	caxis_opt=getfieldvalue(options,'caxis');
@@ -73,5 +79,5 @@
 if strcmpi(colorm,'Rignot'),
 	%hue (H)
-	h_data=(data_grid-data_min)/(data_max-data_min);
+	h_data=(data_grid-data_min)/(data_max-data_min+eps);
 	if radaronly,
 		h_data(:)=0;
@@ -81,5 +87,5 @@
 elseif strcmpi(colorm,'Seroussi'),
 	%hue (H)
-	h_data=1-(data_grid-data_min)/(data_max-data_min)*0.7;
+	h_data=1-(data_grid-data_min)/(data_max-data_min+eps)*0.7;
 	%h_data=(data_grid-data_min)/(data_max-data_min)*2/3;
 	if radaronly,