Index: /issm/trunk/src/m/model/plot/applyoptions.m
===================================================================
--- /issm/trunk/src/m/model/plot/applyoptions.m	(revision 7141)
+++ /issm/trunk/src/m/model/plot/applyoptions.m	(revision 7142)
@@ -89,6 +89,14 @@
 %Basinzoom
 if exist(options,'basin');
-	basinzoom(getfieldvalue(options,'basin'),getfieldvalue(options,'unit',1));
-end
+	if strcmpi(md.hemisphere,'n'), hemisphere=+1; else hemisphere=-1; end
+	basinzoom('basin',getfieldvalue(options,'basin'),'hemisphere',getfieldvalue(options,'hemisphere',hemisphere),'basindeltax',getfieldvalue(options,'basindeltax',300000),'basindeltay',getfieldvalue(options,'basindeltay',300000),'unit',getfieldvalue(options,'unit',1));
+end
+
+%ShowBasins
+if strcmpi(getfieldvalue(options,'showbasins','off'),'on')
+	if strcmpi(md.hemisphere,'n'), hemisphere=+1; else hemisphere=-1; end
+	showbasins('hemisphere',getfieldvalue(options,'hemisphere',hemisphere),'fontsize',fontsize);
+end
+
 
 %Caxis
@@ -202,5 +210,5 @@
 	end
 	if exist(options,'colorbartitle'),
-		set(get(c,'title'),'String',getfieldvalue(options,'colorbartitle'),'Color',getfieldvalue(options,'FontColor','k'));
+		set(get(c,'title'),'FontSize',getfieldvalue(options,'colorbarfontsize',fontsize),'String',getfieldvalue(options,'colorbartitle'),'Color',getfieldvalue(options,'FontColor','k'));
 	end
 	if exist(options,'colorbarwidth'),
@@ -356,26 +364,27 @@
 	%if we are zooming on a basin, don't take the mesh for the boundaries!
 	if exist(options,'basin'),
-		[mdx mdy]=basinzoom(getfieldvalue(options,'basin'));
-	else
-		mdx=[min(md.x) max(md.x)];
-		mdy=[min(md.y) max(md.y)];
-	end
-	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');
+		if strcmpi(md.hemisphere,'n'), hemisphere=+1; else hemisphere=-1; end
+		[mdx mdy]=basinzoom('basin',getfieldvalue(options,'basin'),'hemisphere',getfieldvalue(options,'hemisphere',hemisphere),'basindeltax',getfieldvalue(options,'basindeltax',300000),'basindeltay',getfieldvalue(options,'basindeltay',300000),'unit',getfieldvalue(options,'unit',1));
+else
+	mdx=[min(md.x) max(md.x)];
+	mdy=[min(md.y) max(md.y)];
+end
+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');
 end
 
 %flag edges of a partition
 if exist(options,'partitionedges')
-	[xsegments ysegments]=flagedges(md.elements,md.x,md.y,md.part);
-	xsegments=xsegments*getfieldvalue(options,'unit',1);
-	ysegments=ysegments*getfieldvalue(options,'unit',1);
-	color=getfieldvalue(options,'partitionedgescolor','r-');
-	linewidth=getfieldvalue(options,'linewidth',1);
-	hold on;
-	for i=1:length(xsegments),
-		plot(xsegments(i,:),ysegments(i,:),color,'LineWidth',linewidth);
-	end
+[xsegments ysegments]=flagedges(md.elements,md.x,md.y,md.part);
+xsegments=xsegments*getfieldvalue(options,'unit',1);
+ysegments=ysegments*getfieldvalue(options,'unit',1);
+color=getfieldvalue(options,'partitionedgescolor','r-');
+linewidth=getfieldvalue(options,'linewidth',1);
+hold on;
+for i=1:length(xsegments),
+	plot(xsegments(i,:),ysegments(i,:),color,'LineWidth',linewidth);
+end
 end