Index: /issm/trunk-jpl/src/m/plot/plot_scatter.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_scatter.m	(revision 16052)
+++ /issm/trunk-jpl/src/m/plot/plot_scatter.m	(revision 16053)
@@ -1,5 +1,19 @@
 function plot_scatter(x,y,level,varargin),
+%PLOT_SCATTER - scatter plot
+%
+%   Usage:
+%      plot_scatter(x,y,level,options);
+%
+%   Available options:
+%      'MarkerSize' : default is 3
+%      'Line'       : use line instead of circles
+%      'Cutoff'     : cut the line if the distance between 2 points is
+%                     greater than Cutoff (default is 1000)
 
-options=pairoptions(varargin{:}); 
+if nargin == 4,
+	options = varargin{1};
+else
+	options=pairoptions(varargin{:}); 
+end
 
 %check input
@@ -47,10 +61,26 @@
 hold on
 hp=[];
-for i=1:numcolors
-	pos=find(colorind==i);
-	hprime=plot3(x(pos),y(pos),ones(size(x(pos))),...
-		'o','MarkerSize',getfieldvalue(options,'MarkerSize',3),'MarkerEdgeColor',palette(i,:),...
-		'MarkerFaceColor',palette(i,:));
-	hp=[hp;hprime];
+if ~exist(options,'line'),
+	for i=1:numcolors
+		pos=find(colorind==i);
+		hprime=plot3(x(pos),y(pos),ones(size(x(pos))),...
+			'o','MarkerSize',getfieldvalue(options,'MarkerSize',3),'MarkerEdgeColor',palette(i,:),...
+			'MarkerFaceColor',palette(i,:));
+		hp=[hp;hprime];
+	end
+else
+	distances = sqrt( (x(1:end-1)-x(2:end)).^2 + (y(1:end-1)-y(2:end)).^2);
+	pos=find(distances>getfieldvalue(options,'Cutoff',1000));
+	x(pos,:)=NaN;
+	y(pos,:)=NaN;
+	for j=1:numcolors;
+		pos=find(colorind==j);
+		if(~isempty(pos) & pos(1)==1), pos(1)=[]; end
+		if ~isempty(pos),
+			tempx = [x(pos-1) x(pos) NaN(size(pos))]';
+			tempy = [y(pos-1) y(pos) NaN(size(pos))]';
+			line(tempx(1:end-1),tempy(1:end-1),'color',palette(j,:),'linewidth',getfieldvalue(options,'LineWidth',2));
+		end
+	end
 end