Changeset 6751

Timestamp:

12/17/10 17:42:36 (14 years ago)

Author:

Eric.Larour

Message:

Added plot_riftnumbering: to identify segments in each rift.
Fixed plot_riftrelvel bugs.

Location:

issm/trunk/src/m/model/plot

Files:

• plot_manager.m (modified) (1 diff)
• plot_riftnumbering.m (added)
• plot_riftrelvel.m (modified) (1 diff)

issm/trunk/src/m/model/plot/plot_manager.m

@@ -84 +84 @@
                         plot_riftvel(md,options,nlines,ncols,i);
                         return;
+                case 'riftnumbering',
+                        plot_riftnumbering(md,options,nlines,ncols,i);
+                        return;
+
                 case 'rifts',
                         plot_rifts(md,options,nlines,ncols,i);

issm/trunk/src/m/model/plot/plot_riftrelvel.m

@@ -76 +76 @@
                 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