Index: /issm/trunk-jpl/src/m/plot/plot_overlay.py
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_overlay.py	(revision 17655)
+++ /issm/trunk-jpl/src/m/plot/plot_overlay.py	(revision 17656)
@@ -52,6 +52,15 @@
 	gtif=gdal.Open(inputname)
 	arr=gtif.ReadAsArray()
-	#os.system('rm -rf ./temp.tif')
+	os.system('rm -rf ./temp.tif')
 	
+	if gtif.RasterCount>=3:  # RGB array
+		r=gtif.GetRasterBand(1).ReadAsArray()
+		g=gtif.GetRasterBand(2).ReadAsArray()
+		b=gtif.GetRasterBand(3).ReadAsArray()
+		arr=0.299*r+0.587*g+0.114*b
+
+	# normalize array
+	arr=arr/npy.max(arr.ravel())
+
 	# get parameters from cropped geotiff
 	trans=gtif.GetGeoTransform()
@@ -67,5 +76,5 @@
 	xg,yg=npy.meshgrid(xarr,yarr)
 	if options.exist('basemap'):
-		# get handle to basemap instance
+		# TODO get handle to basemap instance 
 		# handle= functiontogethandle()
 		# create coordinate grid in map projection units (for plotting)