Index: /issm/trunk/src/m/model/plot/showregion.m
===================================================================
--- /issm/trunk/src/m/model/plot/showregion.m	(revision 6367)
+++ /issm/trunk/src/m/model/plot/showregion.m	(revision 6367)
@@ -0,0 +1,31 @@
+function showregion(md,insetpos)
+global ISSM_DIR
+
+%get inset relative position (x,y,width,height)
+%insetpos=getfieldvalue(options,'insetpos',[0.02 0.70 0.18 0.18]);
+
+%get current plos position
+cplotpos=get(gca,'pos');
+%compute inset position
+PosInset=[cplotpos(1)+insetpos(1)*cplotpos(3),cplotpos(2)+insetpos(2)*cplotpos(4), insetpos(3)*cplotpos(3), insetpos(4)*cplotpos(4)];
+axes('pos',PosInset);
+axis equal off
+%box off
+if strcmpi(md.hemisphere,'n') | strcmpi(md.hemisphere,'north'),
+	A=expread([ ISSM_DIR 'projects/Exp/Greenland.exp']);
+elseif strcmpi(md.hemisphere,'s') | strcmpi(md.hemisphere,'south'),
+	A=expread([ ISSM_DIR '/projects/Exp/Antarctica.exp']);
+else
+	error('applyoptions error message: hemisphere not defined');
+end
+
+Ax=[min(A.x) max(A.x)];
+Ay=[min(A.y) max(A.y)];
+
+mdx=[min(md.x) max(md.x)];
+mdy=[min(md.y) max(md.y)];
+
+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','r','FaceAlpha',0.5)
+colorbar('off');