Index: /issm/trunk/src/m/model/plot/plot_none.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_none.m	(revision 6659)
+++ /issm/trunk/src/m/model/plot/plot_none.m	(revision 6660)
@@ -10,4 +10,9 @@
 options=addfielddefault(options,'axis','equal');
 
+if exist(options,'overlay'),
+	plot_overlay(md,'none',options,nlines,ncols,i);
+	return;
+end
+
 %apply options
 applyoptions(md,[],options);
Index: /issm/trunk/src/m/model/plot/plot_overlay.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_overlay.m	(revision 6659)
+++ /issm/trunk/src/m/model/plot/plot_overlay.m	(revision 6660)
@@ -9,5 +9,12 @@
 %process mesh and data
 [x y z elements is2d]=processmesh(md,[],options);
-[data datatype]=processdata(md,data,options);
+if strcmpi(data,'none'),
+	radaronly=1;
+	data=NaN*ones(md.numberofgrids,1);
+	datatype=1;
+else
+	radaronly=0;
+	[data datatype]=processdata(md,data,options);
+end
 
 %check is2d
@@ -30,5 +37,7 @@
 else
 	disp('Extracting radar image...');
-	md=radarpower(md,xlim2,ylim2,getfieldvalue(options,'highres',0));
+	options=addfielddefault(options,'xlim',xlim2);
+	options=addfielddefault(options,'ylim',ylim2);
+	md=radarpower(md,options);
 end
 
@@ -64,4 +73,7 @@
 	%hue (H)
 	h_data=(data_grid-data_min)/(data_max-data_min);
+	if radaronly,
+		h_data(:)=0;
+	end
 	%saturation (S)
 	s_data=max(min((0.1+h_data).^(1/transparency),1),0);
@@ -110,5 +122,7 @@
 %Apply options, without colorbar and without grid
 options=changefieldvalue(options,'colormap',colorm);           %We used an HSV colorbar
-options=changefieldvalue(options,'caxis',[data_min data_max]); %force caxis so that the colorbar is ready
+if ~isnan(data_min),
+	options=changefieldvalue(options,'caxis',[data_min data_max]); %force caxis so that the colorbar is ready
+end
 options=addfielddefault(options,'axis','equal off');           %default axis
 applyoptions(md,data,options);
Index: /issm/trunk/src/m/model/plot/plotdoc.m
===================================================================
--- /issm/trunk/src/m/model/plot/plotdoc.m	(revision 6659)
+++ /issm/trunk/src/m/model/plot/plotdoc.m	(revision 6660)
@@ -96,4 +96,5 @@
 disp('       ''yticklabel'': specifiy yticklabel');
 disp('       ''overlay'': yes or no. This will overlay a radar amplitude image behind');
+disp('       ''overlay_image'': path to overlay image. provide overlay_xlim, overlay_ylim, overlay_xposting and overlay_yposting options also');
 disp('       ''contrast'': (default 1) coefficient to add contrast to the radar amplitude image used in overlays');
 disp('       ''highres'': resolution of overlayed radar amplitude image (default is 0, high resolution is 1).');