Index: /issm/trunk/src/m/classes/public/plot/applyoptions.m
===================================================================
--- /issm/trunk/src/m/classes/public/plot/applyoptions.m	(revision 3173)
+++ /issm/trunk/src/m/classes/public/plot/applyoptions.m	(revision 3174)
@@ -264,9 +264,11 @@
 	axes('pos',PosInset);
 	axis equal off
-	box off
+	%box off
 	if strcmpi(md.hemisphere,'n') | strcmpi(md.hemisphere,'north'),
 		A=expread([ ISSM_DIR '/../models/Greenland/Exp_Par/DomainOutline.exp']);
 	elseif strcmpi(md.hemisphere,'s') | strcmpi(md.hemisphere,'south'),
 		A=expread([ ISSM_DIR '/../models/Antarctica/Exp_Par/DomainOutline.exp']);
+	else
+		error('applyoptions error message: hemisphere not defined');
 	end
 	Ax=[min(A.x) max(A.x)];
@@ -274,6 +276,6 @@
 	mdx=[min(md.x) max(md.x)];
 	mdy=[min(md.y) max(md.y)];
-	patch([Ax(1)  Ax(2)  Ax(2)  Ax(1)],[Ay(1)  Ay(1)  Ay(2)  Ay(2)],[1 1 1],'edgecolor','none')
 	line(A.x,A.y,'color','b');
+	patch([Ax(1)  Ax(2)  Ax(2)  Ax(1) Ax(1)],[Ay(1)  Ay(1)  Ay(2)  Ay(2) Ay(1)],[1 1 1],'EdgeColor',[0 0 0],'LineWidth',1,'FaceLighting','none')
 	patch( [mdx(1) mdx(2) mdx(2) mdx(1)],[mdy(1) mdy(1) mdy(2) mdy(2)],ones(4,1),'EdgeColor',[0 0 0],'FaceColor',[0.5 0.5 0.5],'FaceAlpha',0.5)
 	colorbar('off');
Index: /issm/trunk/src/m/classes/public/plot/plot_overlay.m
===================================================================
--- /issm/trunk/src/m/classes/public/plot/plot_overlay.m	(revision 3173)
+++ /issm/trunk/src/m/classes/public/plot/plot_overlay.m	(revision 3174)
@@ -20,6 +20,6 @@
 
 %get xlim and ylim
-xlim=getfieldvalue(options,'xlim',[min(md.x) max(md.x)]);
-ylim=getfieldvalue(options,'ylim',[min(md.y) max(md.y)]);
+xlim=getfieldvalue(options,'xlim',[min(x) max(x)])/getfieldvalue(options,'unit',1);
+ylim=getfieldvalue(options,'ylim',[min(y) max(y)])/getfieldvalue(options,'unit',1);
 
 %radar power
@@ -31,4 +31,8 @@
 	t2=clock;fprintf('%s\n',[' done (' num2str(etime(t2,t1)) ' seconds)']);
 end
+
+
+xlim=getfieldvalue(options,'xlim',[min(x) max(x)]);
+ylim=getfieldvalue(options,'ylim',[min(y) max(y)]);
 
 %mesh2grid
@@ -101,5 +105,5 @@
 
 %Plot: 
-imagesc(md.sarxm,md.sarym,image_rgb);set(gca,'YDir','normal');
+imagesc(md.sarxm*getfieldvalue(options,'unit',1),md.sarym*getfieldvalue(options,'unit',1),image_rgb);set(gca,'YDir','normal');
 
 %last step: mesh overlay?