Index: sm/trunk/src/m/utils/Plot/basinzoom.m
===================================================================
--- /issm/trunk/src/m/utils/Plot/basinzoom.m	(revision 7143)
+++ 	(revision )
@@ -1,159 +1,0 @@
-function varargout=basinzoom(region,varargin)
-%ANTZOOM - zoom on a basin in Antarctica or Greenland.
-%
-%   This function zooms on an existing figure describing Antarctica or Greenland 
-%   The zooming depends on the region name provided as input. 
-%
-%   Usage:
-%      varargout=basinzoom(region,unitmultiplier)
-
-
-%Initialize regions
-available_regions=InitializeRegionInfo();
-
-%recover unitmultiplier if available
-if nargin==2,
-	unitmultiplier=varargin{1};
-else
-	unitmultiplier=NaN;
-end
-
-
-%Check arguments
-if nargin==0,
-	help basinzoom;
-	disp('   Available regions include:');
-	for i=1:length(available_regions),
-		disp(['     ' available_regions(i).name]);
-	end
-	return;
-end
-
-
-found=0;
-for i=1:length(available_regions),
-	if strcmp(available_regions(i).name,region),
-		%Make sure the region is a square, enlarge it otherwise, except for Antarctica. 
-		x0=available_regions(i).x0;
-		x1=available_regions(i).x1;
-		y0=available_regions(i).y0;
-		y1=available_regions(i).y1;
-		if ~strcmp(region,'antarctica'),
-			dx=x1-x0;
-			dy=y1-y0;
-			if dy>dx,
-				delta=dy-dx;
-				x0=x0-delta/2;
-				x1=x1+delta/2;
-			elseif dy<dx,
-				delta=dx-dy;
-				y0=y0-delta/2;
-				y1=y1+delta/2;
-			else
-				%do nothing;
-			end
-		end
-
-		if ~isnan(unitmultiplier)
-			x0=x0*unitmultiplier;
-			x1=x1*unitmultiplier;
-			y0=y0*unitmultiplier;
-			y1=y1*unitmultiplier;
-		end
-
-		%if output arguments are present, return the limits, 
-		%otherwise, set them on the current graphic. 
-		if nargout==2,
-			found=1;
-			varargout{1}=[x0 x1];
-			varargout{2}=[y0 y1];
-		else
-			xlim([x0 x1]);
-			ylim([y0 y1]);
-			found=1;
-			daspect([1;1;1]);
-		end
-	else
-		continue;
-	end
-end
-
-if found==0,
-	disp(['Unknow region: ',region]);
-	disp('Available regions include:');
-	for i=1:length(available_regions),
-		disp(['   ' available_regions(i).name]);
-	end
-end
-
-end
-
-function incrementedregions=AddAvailableRegion(regions,name,x0,x1,y0,y1);
-
-	if isempty(regions),
-		incrementedregions(1).name=name;
-	else
-		incrementedregions=regions;
-		incrementedregions(end+1).name=name;
-	end
-	incrementedregions(end).x0=x0;
-	incrementedregions(end).x1=x1;
-	incrementedregions(end).y0=y0;
-	incrementedregions(end).y1=y1;
-end
-
-function regions=InitializeRegionInfo()
-	regions=AddAvailableRegion([],'abbot',-2.1473*10^6,-1.6815*10^6,-4.519*10^5,.3785*10^5);
-	regions=AddAvailableRegion(regions,'adelieland',2.699*10^4,2.4636*10^6,-2.1848*10^6,-5.2865*10^5);
-	regions=AddAvailableRegion(regions,'ameryiceshelf',1.5945*10^6,2.3807*10^6,5.1372*10^5,9.0186*10^5);
-	regions=AddAvailableRegion(regions,'amundsencoast',-2.8*10^5,-1*10^5,-5.8*10^5,-4.6*10^5);
-	regions=AddAvailableRegion(regions,'amundsensea',-2.1*10^6,-1.45*10^6,-7.5*10^5,-2.5*10^5);
-	regions=AddAvailableRegion(regions,'antarctica',-2507187,2743299,-2143434,2241152);
-	regions=AddAvailableRegion(regions,'banzarecoast',1.6*10^6,2.3*10^6,-1.9*10^6,-1.2*10^6);
-	regions=AddAvailableRegion(regions,'berknerisland',-1.0709*10^6,-.6405*10^6,5.1823*10^5,9.8177*10^5);
-	regions=AddAvailableRegion(regions,'berkner',-1.0709*10^6,-.6405*10^6,5.1823*10^5,9.8177*10^5);
-	regions=AddAvailableRegion(regions,'borchgrevinkcoast',3*10^5,7*10^5,-2.05*10^6,-1.6*10^6);
-	regions=AddAvailableRegion(regions,'brunt',-8.36*10^5,-4.4*10^5,1.33*10^6,1.64*10^6);
-	regions=AddAvailableRegion(regions,'bryan',-1.9727*10^6,-1.7421*10^6,-0.4427*10^5,1.9010*10^5);
-	regions=AddAvailableRegion(regions,'buddcoast',2.1*10^6,2.7*10^6,-1.3*10^6,-.7*10^6);
-	regions=AddAvailableRegion(regions,'byrdglacier',2.5*10^5,4.5*10^5,-10.8*10^5,-9.2*10^5);
-	regions=AddAvailableRegion(regions,'drygalski',3.7916*10^5,7.2246*10^5,-1.7266*10^6,-1.3724*10^6);
-	regions=AddAvailableRegion(regions,'dufekcoast',-1.6*10^5,0*10^5,-6.2*10^5,-5.5*10^5);
-	regions=AddAvailableRegion(regions,'edouardpeninsula',-7.5*10^5,-4*10^5,-1.4*10^6,-1.05*10^6);
-	regions=AddAvailableRegion(regions,'filchnericeshelf',-9.3781*10^5,-3.6671*10^5,.7389*10^6,1.0360*10^6);
-	regions=AddAvailableRegion(regions,'fimbulisen',-2.8693*10^5,3.2792*10^5,1.9818*10^6,2.2891*10^6);
-	regions=AddAvailableRegion(regions,'fordranges',-8.5*10^5,-5.7*10^5,-1.33*10^6,-1.12*10^6);
-	regions=AddAvailableRegion(regions,'georgevcoast',.9*10^6,1.7*10^6,-2.4*10^6,-1.7*10^6);
-	regions=AddAvailableRegion(regions,'getziceshelf',-1.7*10^6,-1.15*10^6,-1.23*10^6,-.5*10^6);
-	regions=AddAvailableRegion(regions,'hillarycoast',2.6*10^5,4.4*10^5,-1.24*10^6,-1.06*10^6);
-	regions=AddAvailableRegion(regions,'icestreams',-1.2274*10^6,-.137*10^6,-1.1172*10^6,-.138*10^6);
-	regions=AddAvailableRegion(regions,'icestreamE',2.2*10^6,2.7*10^6,0,5*10^5);
-	regions=AddAvailableRegion(regions,'jks',-2.4*10^5,-.8*10^5,-2.3*10^6,-2.2*10^6);
-	regions=AddAvailableRegion(regions,'jakobshavn',-2.4*10^5,-.8*10^5,-2.3*10^6,-2.2*10^6);
-	regions=AddAvailableRegion(regions,'larseniceshelf',-2.3855*10^6,-1.9649*10^6,0.9498*10^6,1.2996*10^6);
-	regions=AddAvailableRegion(regions,'larc',-2.3855*10^6,-1.9649*10^6,0.9498*10^6,1.2996*10^6);
-	regions=AddAvailableRegion(regions,'larc2',-2.332760*10^6,-1.995350*10^6,0.9895*10^6,1.28559*10^6);
-	regions=AddAvailableRegion(regions,'mcmurdo',1.3750*10^5,4.8255*10^5,-1.4569*10^6,-1.1119*10^6);
-	regions=AddAvailableRegion(regions,'nimrodglacier',1*10^5,3.8*10^5,-9*10^5,-7*10^5);
-	regions=AddAvailableRegion(regions,'oatescoast',5.3*10^5,8.5*10^5,-2.15*10^6,-1.9*10^6);
-	regions=AddAvailableRegion(regions,'peninsula',-2.7617*10^6,-1.2831*10^6,.4492*10^6,1.6940*10^6);
-	regions=AddAvailableRegion(regions,'pennelcoast',3*10^5,6*10^5,-2.1*10^6,-1.85*10^6); 
-	regions=AddAvailableRegion(regions,'petermann',-3.8*10^5,-.9*10^5,-1.1*10^6,-.88*10^6); 
-	regions=AddAvailableRegion(regions,'pet',-3.8*10^5,-.9*10^5,-1.1*10^6,-.88*10^6); 
-	regions=AddAvailableRegion(regions,'pig',-1.72*10^6,-1.45*10^6,-3.5*10^5,.5*10^5);
-	regions=AddAvailableRegion(regions,'ronneiceshelf',-1.5*10^6,-.7*10^6,1*10^5,9.3*10^5);
-	regions=AddAvailableRegion(regions,'ris',-1.5*10^6,-.7*10^6,1*10^5,9.3*10^5);
-	regions=AddAvailableRegion(regions,'roosevelt',-6.2797*10^5,-1.5414*10^5,-1.2734*10^6,-.9609*10^6);
-	regions=AddAvailableRegion(regions,'rossiceshelf',-8*10^5,4.5*10^5,-13.8*10^5,-5*10^5);
-	regions=AddAvailableRegion(regions,'rois',-8*10^5,4.5*10^5,-13.8*10^5,-5*10^5);
-	regions=AddAvailableRegion(regions,'ruppertcoast',-11*10^5,-8.5*10^5,-1.34*10^6,-1.14*10^6);
-	regions=AddAvailableRegion(regions,'scottcoast',2.5*10^5,6*10^5,-1.6*10^6,-1.2*10^6);
-	regions=AddAvailableRegion(regions,'shackletoncoast',2*10^5,3.8*10^5,-9.8*10^5,-8.4*10^5);
-	regions=AddAvailableRegion(regions,'shackletoniceshelf',2.4*10^6,2.8*10^6,-7*10^5,-1*10^5);
-	regions=AddAvailableRegion(regions,'shirasecoast',-6*10^5,-4.4*10^5,-10.9*10^5,-8.2*10^5);
-	regions=AddAvailableRegion(regions,'siplecoast',-4.3*10^5,-1.5*10^5,-9*10^5,-6.5*10^5);
-	regions=AddAvailableRegion(regions,'umq',-5*10^5,-3*10^5,-2.1*10^6,-1.8*10^6);
-	regions=AddAvailableRegion(regions,'wais',-1.8804*10^6,-1.0965*10^6,-6.5885*10^5,.8594*10^5);
-	regions=AddAvailableRegion(regions,'westiceshelf',2.2*10^6,2.7*10^6,0,5*10^5);
-
-end