Changeset 7082

Timestamp:

01/14/11 08:10:45 (14 years ago)

Author:

Eric.Larour

Message:

better options for these plots

Location:

issm/trunk/src/m/model/plot

Files:

• plot_riftpenetration.m (modified) (1 diff)
• plot_riftrelvel.m (modified) (4 diffs)

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

@@ -15 +15 @@
 %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
 

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

@@ -16 +16 @@
 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 +25 @@
 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 +40 @@
 %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 +77 @@
                 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