Changeset 10432

Timestamp:

11/02/11 12:05:27 (13 years ago)

Author:

Mathieu Morlighem

Message:

Support for 3d models

File:

• issm/trunk/src/m/model/plot/plot_referential.m (modified) (1 diff)

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

@@ -48 +48 @@
         yend=ystart+len*Yhat(:,2);
         hy=quiver(xstart,ystart,xend-xstart,yend-ystart,'Color','red','ShowArrowHead','on','AutoScale','off');
+
+        legend([hx,hy],'local X direction','local Y direction')
 else
-        error('3d not supported yet')
+        %plot mesh
+        A=elements(:,1); B=elements(:,2); C=elements(:,3); D=elements(:,4); E=elements(:,5); F=elements(:,6);
+        h1=patch( 'Faces', [A B C],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor','black');
+        patch( 'Faces', [D E F],  'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor','black');
+        patch( 'Faces', [A B E D],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor','black');
+        patch( 'Faces', [B E F C],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor','black');
+        patch( 'Faces', [C A D F],'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','none','EdgeColor','black');
+        hold on;
+
+        xstart=x;
+        ystart=y;
+        zstart=z;
+        edgex=max(md.mesh.x(elements),[],2)-min(md.mesh.x(elements),[],2);
+        edgez=max(md.mesh.z(elements),[],2)-min(md.mesh.z(elements),[],2);
+        len=min(edgex)/1.5;
+        lenz=min(edgez)/1.5;
+        %plot X axis
+        xend=xstart+len*Xhat(:,1);
+        yend=ystart+len*Xhat(:,2);
+        zend=zstart+len*Xhat(:,3);
+        hx=quiver3(xstart,ystart,zstart,xend-xstart,yend-ystart,zend-zstart,'Color','blue','ShowArrowHead','on','AutoScale','off');
+        %plot Y axis
+        xend=xstart+len*Yhat(:,1);
+        yend=ystart+len*Yhat(:,2);
+        zend=zstart+len*Yhat(:,3);
+        hy=quiver3(xstart,ystart,zstart,xend-xstart,yend-ystart,zend-zstart,'Color','red','ShowArrowHead','on','AutoScale','off');
+        %plot Z axis
+        xend=xstart+lenz*Zhat(:,1);
+        yend=ystart+lenz*Zhat(:,2);
+        zend=zstart+lenz*Zhat(:,3);
+        hz=quiver3(xstart,ystart,zstart,xend-xstart,yend-ystart,zend-zstart,'Color','green','ShowArrowHead','on','AutoScale','off');
+
+        legend([hx,hy,hz],'local X direction','local Y direction','local Z direction')
 end
-
-%legend (disable warnings)
-legend([hx,hy],'local X direction','local Y direction')
 
 %apply options