Index: /issm/trunk-jpl/src/m/plot/radarpower.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/radarpower.m	(revision 16436)
+++ /issm/trunk-jpl/src/m/plot/radarpower.m	(revision 16437)
@@ -25,4 +25,8 @@
 ylim    = getfieldvalue(options,'ylim',[min(md.mesh.y) max(md.mesh.y)]);
 posting = getfieldvalue(options,'posting',0); % 0 -> image posting default
+a = getfieldvalue(options,'overlay_adjust_a',0);
+b = getfieldvalue(options,'overlay_adjust_b',1);
+c = getfieldvalue(options,'overlay_adjust_c',0);
+d = getfieldvalue(options,'overlay_adjust_d',1);
 
 %find gdal coordinates
@@ -85,4 +89,6 @@
 		%Read in temp.tif:
 		im=imread('temp.tif','TIFF');
+		%adjust contrast and brightness
+		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)));
@@ -114,4 +120,6 @@
 		%Read in temp.tif:
 		im=imread('temp.tif','TIFF');
+		%adjust contrast and brightness
+		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)));
@@ -144,4 +152,6 @@
 	disp('Warning: expecting coordinates in polar stereographic (Std Latitude: 70ºN Meridian: 45º)');
 	im=imread(overlay_image);
+	%adjust contrast and brightness
+	im=imadjust(im,[a b],[c d]);
 	im=im(toplefty:toplefty+sizey,topleftx:topleftx+sizex);
 	md.radaroverlay.pwr=double(flipud(im));