Index: /issm/trunk-jpl/src/m/plot/plot_overlay.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_overlay.m	(revision 17223)
+++ /issm/trunk-jpl/src/m/plot/plot_overlay.m	(revision 17224)
@@ -43,4 +43,8 @@
 radar = (md.radaroverlay.pwr).^(contrast);
 radar = radar./max(radar(:));
+if size(radar,3)>1,
+	disp('WARNING: color image converted to greyscale intensity image');
+	radar=sum(radar,3)/3;
+end
 %radar(find(radar==0))=1; %Change background from black to white
 
Index: /issm/trunk-jpl/src/m/plot/radarpower.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/radarpower.m	(revision 17223)
+++ /issm/trunk-jpl/src/m/plot/radarpower.m	(revision 17224)
@@ -90,7 +90,14 @@
 		im=imread('temp.tif','TIFF');
 		%adjust contrast and brightness
-		im=imadjust(im,[a b],[c d]);
+		%im=imadjust(im,[a b],[c d]);
 		pixelskip=max(1,ceil(posting/((x1-x0)/(size(im,2)))));
-		md.radaroverlay.pwr=double(flipud(im(1:pixelskip:end,1:pixelskip:end)));
+		if size(im,3)==1,
+			md.radaroverlay.pwr=double(flipud(im(1:pixelskip:end,1:pixelskip:end)));
+		else
+			md.radaroverlay.pwr=double(im(1:pixelskip:end,1:pixelskip:end,1:pixelskip:end));
+			md.radaroverlay.pwr(:,:,1)=flipud(md.radaroverlay.pwr(:,:,1));
+			md.radaroverlay.pwr(:,:,2)=flipud(md.radaroverlay.pwr(:,:,2));
+			md.radaroverlay.pwr(:,:,3)=flipud(md.radaroverlay.pwr(:,:,3));
+		end
 		md.radaroverlay.x=(x0:(x1-x0)/(size(md.radaroverlay.pwr,2)-1):x1);
 		md.radaroverlay.y=(y0:(y1-y0)/(size(md.radaroverlay.pwr,1)-1):y1);
@@ -153,5 +160,5 @@
 	im=imread(overlay_image);
 	%adjust contrast and brightness
-	im=imadjust(im,[a b],[c d]);
+	%im=imadjust(im,[a b],[c d]);
 	im=im(toplefty:toplefty+sizey,topleftx:topleftx+sizex);
 	md.radaroverlay.pwr=double(flipud(im));