Index: /issm/trunk-jpl/src/m/plot/applyoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/applyoptions.m	(revision 13835)
+++ /issm/trunk-jpl/src/m/plot/applyoptions.m	(revision 13836)
@@ -265,4 +265,10 @@
 if exist(options,'northarrow'),
 	northarrow(getfieldvalue(options,'northarrow'));
+end
+
+%curved arrow
+if exist(options,'curvedarrow'),
+	curvedoptions=getfieldvalue(options,'curvedarrow');
+	curvedarrow(curvedoptions{:});
 end
 
Index: /issm/trunk-jpl/src/m/plot/colormaps/getcolormap.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/colormaps/getcolormap.m	(revision 13835)
+++ /issm/trunk-jpl/src/m/plot/colormaps/getcolormap.m	(revision 13836)
@@ -21,5 +21,5 @@
 
 if strcmpi(map,'Ala'),
-	map = jet;
+	map = jet(256);
 	map = map(128:end,:);
 elseif strcmpi(map,'redblue'),
Index: /issm/trunk-jpl/src/m/plot/curvedarrow.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/curvedarrow.m	(revision 13836)
+++ /issm/trunk-jpl/src/m/plot/curvedarrow.m	(revision 13836)
@@ -0,0 +1,60 @@
+function curvedarrow(centerx,centery,distance,angle,secondangle,varargin)
+%CURVEDARROW - plot curved arrow, where curvature center is on (centerx,centery), radius is 'distance', 
+%   and angle determines arc length (angle in degrees). secondangle determines by how much the arc 
+%   generated is rotated along the [0,0,1] z axis.30,
+%
+%   Usage:
+%      curevedarrow(x1,y1,r,30,30,options)
+%      where options is a lit of paired arguments of string OR enums
+%      options can be: 
+%            'ratio': default .5 (ratio headarrow/length)
+%            'widthratio': default is 1/10 of length
+
+	%recover options
+	options=pairoptions(varargin{:});
+	ratio=getfieldvalue(options,'ratio',.1);
+	widthratio=getfieldvalue(options,'widthratio',.1);
+	color=getfieldvalue(options,'color','k');
+
+	%transform angle in radians
+	angle=angle/180*pi;
+	nsteps=10;
+
+	%compute some values out of (x1,y1) and (x2,y2)
+	length=distance*angle;
+	width=length*widthratio;
+
+	%buidl the arrow itself: 
+	A=[centerx+distance, centery];
+	B=[centerx+distance+width, centery];
+	BC=[centerx+cos(0:angle/nsteps:angle)*(distance+width); centery+sin(0:angle/nsteps:angle)*(distance+width)]';
+	C=[centerx+cos(angle)*(distance+width), centery+sin(angle)*(distance+width)];
+	D=[centerx+cos(angle)*(distance), centery+sin(angle)*(distance)];
+	DA=[centerx+cos(angle:-angle/nsteps:0)*(distance); centery+sin(angle:-angle/nsteps:0)*(distance)]';
+
+	%Plot arrow
+	hold on
+	p1=patch([A(1) B(1) BC(:,1)' C(1) D(1) DA(:,1)' ],[A(2) B(2) BC(:,2)' C(2) D(2) DA(:,2)'],color);
+	set(p1,'EdgeColor',color); set(p1,'FaceColor',color);
+
+	%Build arrowhead 
+	E=D+2/3*(D-C);
+	F=C+2/3*(C-D);
+	n=(F-E)/norm(F-E,2);
+	m=[-n(2) n(1)];
+	if(angle<0)
+		m=-m;
+	end
+	d=(distance*angle)*ratio;
+	
+	%G is d distance from middle of E and F: 
+	G=(E+F)/2+d*m;
+
+	p2=patch([E(1) F(1) G(1)],[E(2) F(2) G(2)],color);
+	set(p2,'EdgeColor',color); set(p2,'FaceColor',color);
+
+	%now rotate our arrow: 
+	rotate(p1,[0 0  1],secondangle);
+	rotate(p2,[0 0 1],secondangle);
+
+end
Index: /issm/trunk-jpl/src/m/plot/plot_riftrelvel.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_riftrelvel.m	(revision 13835)
+++ /issm/trunk-jpl/src/m/plot/plot_riftrelvel.m	(revision 13836)
@@ -15,4 +15,7 @@
 end
 options=addfielddefault(options,'scaling',2);
+
+%markersize: 
+markersize=getfieldvalue(options,'markersize',1);
 
 %recover vx and vy:
@@ -37,9 +40,4 @@
 subplot(nlines,ncols,index); 
 hold on
-
-%plot mesh boundaries
-for i=1:size(md.mesh.segments,1),
-	plot(x(md.mesh.segments(i,1:2)),y(md.mesh.segments(i,1:2)),'k-');
-end
 
 %plot rifts vel
@@ -81,8 +79,8 @@
 		%if penetration is negative, plot in black, positive, plot in red;: ie: if rift is closing, black, if rift is opening, red.
 		if(penetration>0),
-			p2=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'*'); set(p2,'Color',[140 140 140]/255);
+			p2=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'*','MarkerSize',markersize); set(p2,'Color',[140 140 140]/255);
 			isp2=1;
 		else
-			p1=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'k*');
+			p1=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'k*','MarkerSize',markersize);
 			isp1=1;
 		end
@@ -90,6 +88,6 @@
 
 	%point out the tips
-	h2=plot(x(md.rifts.riftstruct(i).tips(1)),y(md.rifts.riftstruct(i).tips(1)),'g*');
-	plot(x(md.rifts.riftstruct(i).tips(2)),y(md.rifts.riftstruct(i).tips(2)),'g*');
+	h2=plot(x(md.rifts.riftstruct(i).tips(1)),y(md.rifts.riftstruct(i).tips(1)),'g*','MarkerSize',markersize);
+	plot(x(md.rifts.riftstruct(i).tips(2)),y(md.rifts.riftstruct(i).tips(2)),'g*','MarkerSize',markersize);
 	segments=md.rifts.riftstruct(i).segments(:,1:2);
 end