Changeset 3114

Timestamp:

02/24/10 10:37:19 (15 years ago)

Author:

Mathieu Morlighem

Message:

Changed segments plot (much faster, pressureload-like)

File:

• issm/trunk/src/m/classes/public/plot/plot_segments.m (modified) (1 diff)

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

@@ -10 +10 @@
 subplot(width,width,i); 
 
-%units
-if exist(options,'unit'),
-        unit=getfieldvalue(options,'unit');
-        md.x=md.x*unit;
-        md.y=md.y*unit;
-        md.z=md.z*unit;
-end
-
+%process mesh and data
+[x y z elements is2d]=processmesh(md,options);
 segments=md.segments;
 
-length=sqrt( (md.x(segments(:,1))-md.x(segments(:,2))).^2 + (md.y(segments(:,1))-md.y(segments(:,2))).^2 );
-%normal(:,1)=cos(atan2( (md.x(segments(:,1))-md.x(segments(:,2))) , (md.y(segments(:,2))-md.y(segments(:,1))) ) );
-%normal(:,2)=sin(atan2( (md.x(segments(:,1))-md.x(segments(:,2))) , (md.y(segments(:,2))-md.y(segments(:,1))) ) );
+if strcmpi(md.type,'2d'),
+        %plot mesh
+        A=elements(:,1); B=elements(:,2); C=elements(:,3); 
+        h1=patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor','black');
+        hold on;
 
-if strcmpi(md.type,'2d'),
-        for n=1:size(segments,1),
+        %highlight elements on neumann
+        pos=segments(:,end);
+        A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); 
+        h2=patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','green','EdgeColor','black');
 
-                plot(md.x(segments(n,1:2)),md.y(segments(n,1:2)),'b-','linewidth',2);
-                hold on
-                plot(md.x(segments(n,1:2)),md.y(segments(n,1:2)),'r.','markersize',12);
-
-                %xstart=mean(md.x(segments(n,1:2)));
-                %ystart=mean(md.y(segments(n,1:2)));
-
-                %xend=xstart+length(n)*normal(n,1);
-                %yend=ystart+length(n)*normal(n,2);
-
-                %plot([xstart xend],[ystart yend],'r-')
-                legend('element edges on ice front','grids on ice front','normal vectors')
-        end
+        %display arrows pointing outward
+        xstart=mean(x(segments(:,1:end-1)),2);
+        ystart=mean(y(segments(:,1:end-1)),2);
+        length=sqrt((x(segments(:,1))-x(segments(:,2))).^2 + (y(segments(:,1))-y(segments(:,2))).^2 );
+        normal(:,1)=cos(atan2((x(segments(:,1))-x(segments(:,2))) , (y(segments(:,2))-y(segments(:,1)))));
+        normal(:,2)=sin(atan2((x(segments(:,1))-x(segments(:,2))) , (y(segments(:,2))-y(segments(:,1)))));
+        xend=xstart+length.*normal(:,1);
+        yend=ystart+length.*normal(:,2);
+        q=quiver(xstart,ystart,xend-xstart,yend-ystart); hold on;
+        h3=plot(xstart,ystart,'r*');
 else
         error('plot_segments: 3d plot of segments not supported yet!');
 end
+
+%legend (disable warnings)
+warning off
+legend([h2,q],'element on segment','normal vectors')
+warning on
 
 %apply options