Changeset 3202 for issm/trunk/src/m/classes/public/plot/plot_riftrelvel.m

Timestamp:

03/05/10 14:49:35 (15 years ago)

Author:

Eric.Larour

Message:

corrections for plots

File:

• issm/trunk/src/m/classes/public/plot/plot_riftrelvel.m (modified) (7 diffs)

issm/trunk/src/m/classes/public/plot/plot_riftrelvel.m

@@ -1 @@
-function plot_riftrelvel(md,options,width,i);
+function plot_riftrelvel(md,options,nlines,ncols,index);
 %PLOT_RIFTRELVEL - plot rift relative velocities
 %
 %   Usage:
-%      plot_riftrelvel(md,options,width,i);
+%      plot_riftrelvel(md,options,nlines,ncols,i);
 %
 %   See also: PLOTMODEL
@@ -17 +17 @@
 
 %set as NaN all velocities not on rifts
-for i=1:size(md.rifts,1),
+u=NaN*ones(md.numberofgrids,1);
+v=NaN*ones(md.numberofgrids,1);
+for i=1:md.numrifts,
         penaltypairs=md.rifts(i).penaltypairs(:,[1 2]);
-        u=NaN*ones(md.numberofgrids,1);
-        v=NaN*ones(md.numberofgrids,1);
         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));
@@ -31 +31 @@
 
 %prepare plot
-subplot(width,width,i); 
+subplot(nlines,ncols,index); 
 hold on
 
@@ -37 +37 @@
 for i=1:size(md.segments,1),
         plot(x(md.segments(i,1:2)),y(md.segments(i,1:2)),'k.-');
+end
+
+%plot rifts vel
+h3=[];
+for i=1:quivers.numcolors
+        pos=find(quivers.colorind==i);
+        hprime=quiver(quivers.x(pos),quivers.y(pos),quivers.u(pos),quivers.v(pos),...
+                'Color',palette(i,:),'ShowArrowHead','on','AutoScale','off');
+        hprime=quiver(quivers.x(pos),quivers.y(pos),-quivers.u(pos),-quivers.v(pos),...
+                'Color',palette(i,:),'ShowArrowHead','on','AutoScale','off');
+        h3=[h3;hprime];
 end
 
@@ -47 +58 @@
         h1=plot(x(md.segments(i,1:2)),y(md.segments(i,1:2)),'b-');
 end
-
-for i=1:size(md.rifts,1),
+for i=1:md.numrifts,
+        
         %get grids on rift
         penaltypairs=md.rifts(i).penaltypairs;
@@ -66 +77 @@
                 %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)),'r*');
+                        p2=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'*'); set(p2,'Color',[140 140 140]/255);
                         isp2=1;
                 else
@@ -79 +90 @@
         segments=md.rifts(i).segments(:,1:2);
 end
-
-%plot rifts vel
-h3=[];
-for i=1:quivers.numcolors
-        pos=find(quivers.colorind==i);
-        hprime=quiver(quivers.x(pos),quivers.y(pos),quivers.u(pos),quivers.v(pos),...
-                'Color',palette(i,:),'ShowArrowHead','on','AutoScale','off');
-        hprime=quiver(quivers.x(pos),quivers.y(pos),-quivers.u(pos),-quivers.v(pos),...
-                'Color',palette(i,:),'ShowArrowHead','on','AutoScale','off');
-        h3=[h3;hprime];
-end
-
 %legend
 if isp1 & isp2