Index: /issm/trunk/src/m/model/plot/plot_riftpenetration.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_riftpenetration.m	(revision 7081)
+++ /issm/trunk/src/m/model/plot/plot_riftpenetration.m	(revision 7082)
@@ -15,5 +15,5 @@
 %plot mesh boundaries
 for i=1:size(md.segments,1),
-	plot(x(md.segments(i,1:2)),y(md.segments(i,1:2)),'k.-');
+	plot(x(md.segments(i,1:2)),y(md.segments(i,1:2)),'k-');
 end
 
Index: /issm/trunk/src/m/model/plot/plot_riftrelvel.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_riftrelvel.m	(revision 7081)
+++ /issm/trunk/src/m/model/plot/plot_riftrelvel.m	(revision 7082)
@@ -16,4 +16,8 @@
 options=addfielddefault(options,'scaling',2);
 
+%recover vx and vy:
+vx=getfieldvalue(options,'riftrelvel_vx',md.vx);
+vy=getfieldvalue(options,'riftrelvel_vy',md.vy);
+
 %set as NaN all velocities not on rifts
 u=NaN*ones(md.numberofgrids,1);
@@ -21,6 +25,6 @@
 for i=1:md.numrifts,
 	penaltypairs=md.rifts(i).penaltypairs(:,[1 2]);
-	u(md.rifts(i).penaltypairs(:,1))=md.vx(penaltypairs(:,1))-md.vx(penaltypairs(:,2));
-	v(md.rifts(i).penaltypairs(:,1))=md.vy(penaltypairs(:,1))-md.vy(penaltypairs(:,2));
+	u(md.rifts(i).penaltypairs(:,1))=vx(penaltypairs(:,1))-vx(penaltypairs(:,2));
+	v(md.rifts(i).penaltypairs(:,1))=vy(penaltypairs(:,1))-vy(penaltypairs(:,2));
 end
 
@@ -36,5 +40,5 @@
 %plot mesh boundaries
 for i=1:size(md.segments,1),
-	plot(x(md.segments(i,1:2)),y(md.segments(i,1:2)),'k.-');
+	plot(x(md.segments(i,1:2)),y(md.segments(i,1:2)),'k-');
 end
 
@@ -73,20 +77,20 @@
 		normal(2)=penaltypairs(j,6);
 
-		vx1=md.vx(penaltypairs(j,1)); vx2=md.vx(penaltypairs(j,2)); vy1=md.vy(penaltypairs(j,1)); vy2=md.vy(penaltypairs(j,2));
+		vx1=vx(penaltypairs(j,1)); vx2=vx(penaltypairs(j,2)); vy1=vy(penaltypairs(j,1)); vy2=vy(penaltypairs(j,2));
 		penetration=(vx2-vx1)*normal(1)+(vy2-vy1)*normal(2);
 		%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);
-			isp2=1;
-		else
-			p1=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'k*');
-			isp1=1;
-		end
+		%if(penetration>0),
+		%	p2=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'*'); set(p2,'Color',[140 140 140]/255);
+		%	isp2=1;
+		%else
+		%	p1=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'k*');
+		%	isp1=1;
+		%end
 	end
 
 	%point out the tips
-	h2=plot(x(md.rifts(i).tips(1)),y(md.rifts(i).tips(1)),'g*');
-	plot(x(md.rifts(i).tips(2)),y(md.rifts(i).tips(2)),'g*');
-	segments=md.rifts(i).segments(:,1:2);
+	%h2=plot(x(md.rifts(i).tips(1)),y(md.rifts(i).tips(1)),'g*');
+	%plot(x(md.rifts(i).tips(2)),y(md.rifts(i).tips(2)),'g*');
+	%segments=md.rifts(i).segments(:,1:2);
 end
 %legend