Changeset 7317 for issm/trunk
- Timestamp:
- 02/03/11 16:23:09 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/m/model/plot/plot_riftrelvel.m
r7082 r7317 80 80 penetration=(vx2-vx1)*normal(1)+(vy2-vy1)*normal(2); 81 81 %if penetration is negative, plot in black, positive, plot in red;: ie: if rift is closing, black, if rift is opening, red. 82 %if(penetration>0),83 %p2=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'*'); set(p2,'Color',[140 140 140]/255);84 %isp2=1;85 %else86 %p1=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'k*');87 %isp1=1;88 %end82 if(penetration>0), 83 p2=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'*'); set(p2,'Color',[140 140 140]/255); 84 isp2=1; 85 else 86 p1=plot(x(penaltypairs(j,1)) ,y(penaltypairs(j,1)),'k*'); 87 isp1=1; 88 end 89 89 end 90 90 91 91 %point out the tips 92 %h2=plot(x(md.rifts(i).tips(1)),y(md.rifts(i).tips(1)),'g*');93 %plot(x(md.rifts(i).tips(2)),y(md.rifts(i).tips(2)),'g*');94 %segments=md.rifts(i).segments(:,1:2);92 h2=plot(x(md.rifts(i).tips(1)),y(md.rifts(i).tips(1)),'g*'); 93 plot(x(md.rifts(i).tips(2)),y(md.rifts(i).tips(2)),'g*'); 94 segments=md.rifts(i).segments(:,1:2); 95 95 end 96 96 %legend
Note:
See TracChangeset
for help on using the changeset viewer.