Index: /issm/trunk/src/m/model/plot/applyoptions.m
===================================================================
--- /issm/trunk/src/m/model/plot/applyoptions.m	(revision 7315)
+++ /issm/trunk/src/m/model/plot/applyoptions.m	(revision 7316)
@@ -6,4 +6,11 @@
 %
 %   See also: PLOTMODEL, PARSE_OPTIONS
+		
+
+%some defaults
+if strcmpi(md.hemisphere,'n'), options=addfielddefault(options,'hemisphere','n'); 
+elseif strcmpi(md.hemisphere,'s'), options=addfielddefault(options,'hemisphere','s'); 
+else error('hemisphere field not present');
+end
 
 %fontsize
@@ -89,12 +96,10 @@
 %Basinzoom
 if exist(options,'basin');
-	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));
+	basinzoom(options);
 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);
+	showbasins(options);
 end
 
@@ -364,6 +369,5 @@
 	%if we are zooming on a basin, don't take the mesh for the boundaries!
 	if exist(options,'basin'),
-		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));
+		[mdx mdy]=basinzoom(options);
 	else
 		mdx=[min(md.x) max(md.x)];