Index: /issm/trunk/src/m/model/plot/plot_overlay.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_overlay.m	(revision 9225)
+++ /issm/trunk/src/m/model/plot/plot_overlay.m	(revision 9226)
@@ -115,11 +115,4 @@
 %Select plot area 
 subplot(plotlines,plotcols,i);
-%uncomment following lines to have more space
-%P=get(gca,'pos');
-%P(3)=P(3)+0.07;
-%P(4)=P(4)+0.07;
-%%P(2)=P(2)+0.08;
-%P(1)=P(1)-0.03;
-%set(gca,'pos',P);
 
 %Plot: 
Index: /issm/trunk/src/m/model/radarpower.m
===================================================================
--- /issm/trunk/src/m/model/radarpower.m	(revision 9225)
+++ /issm/trunk/src/m/model/radarpower.m	(revision 9226)
@@ -11,5 +11,5 @@
 
 %If gdal does not work, uncomment the following line
-%setenv('LD_LIBRARY_PATH','');
+setenv('LD_LIBRARY_PATH','');
 
 %Parse inputs
@@ -21,4 +21,5 @@
 xlim=getfieldvalue(options,'xlim',[min(md.x) max(md.x)]);
 ylim=getfieldvalue(options,'ylim',[min(md.y) max(md.y)]);
+posting=getfieldvalue(options,'posting',0); % 0 -> image posting default
 
 %find gdal coordinates
@@ -45,4 +46,5 @@
 		numcols=floor((xmax-xmin)/geom(1)); % x posting
 		numrows=floor((ymax-ymin)/geom(4)); % y posting
+		pixelskip=max(1,ceil(posting/geom(1)));
 
 		%Read and crop file
@@ -50,5 +52,5 @@
 		im=imread(jpgim);
 		im=im(firstrow:firstrow+numrows-1,firstcol:firstcol+numcols-1);
-		md.sarpwr=double(flipud(im));
+		md.sarpwr=double(flipud(im(1:pixelskip:end,1:pixelskip:end)));
 		md.sarxm=(xmin:(xmax-xmin)/(size(md.sarpwr,2)-1):xmax);
 		md.sarym=(ymin:(ymax-ymin)/(size(md.sarpwr,1)-1):ymax);
@@ -72,5 +74,7 @@
 
 		%Read in temp.tif:
-		md.sarpwr=double(flipud(imread('temp.tif','TIFF')));
+		im=imread('temp.tif','TIFF');
+		pixelskip=max(1,ceil(posting/((x1-x0)/(size(im,2)))));
+		md.sarpwr=double(flipud(im(1:pixelskip:end,1:pixelskip:end)));
 		md.sarxm=(x0:(x1-x0)/(size(md.sarpwr,2)-1):x1);
 		md.sarym=(y0:(y1-y0)/(size(md.sarpwr,1)-1):y1);