Changeset 3103

Timestamp:

02/23/10 16:27:41 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added nice colors

File:

• issm/trunk/src/m/classes/public/plot/plot_pressureload.m (modified) (4 diffs)

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

@@ -21 +21 @@
         hold on;
 
-        %highlight elements on neumann in Green
-        pos=pressureload(:,end-1);
+        %highlight elements on neumann
+        pos=find(pressureload(:,end)==WaterEnum());
+        pos=pressureload(pos,end-1);
         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');
-        hold on;
+        h2=patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','blue','EdgeColor','black');
+        pos=find(pressureload(:,end)==AirEnum());
+        pos=pressureload(pos,end-1);
+        A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); 
+        h3=patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','green','EdgeColor','black');
 
         %display arrows pointing outward
@@ -36 +40 @@
         yend=ystart+length.*normal(:,2);
         q=quiver(xstart,ystart,xend-xstart,yend-ystart); hold on;
-        h3=plot(xstart,ystart,'r*');
+        h4=plot(xstart,ystart,'r*');
 else
 
@@ -48 +52 @@
         hold on;
 
-        %highlight elements on neumann in Green
-        pos=pressureload(:,end-1);
+        %highlight elements on neumann
+        pos=find(pressureload(:,end)==WaterEnum());
+        pos=pressureload(pos,end-1);
         A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); D=elements(pos,4); E=elements(pos,5); F=elements(pos,6);
-        h2=patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','green','EdgeColor','black');
+        h2=patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','blue','EdgeColor','black');
+        patch( 'Faces', [D E F], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','blue','EdgeColor','black');
+        patch( 'Faces', [A B E D], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','blue','EdgeColor','black');
+        patch( 'Faces', [B E F C ], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','blue','EdgeColor','black');
+        patch( 'Faces', [C A D F ], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','blue','EdgeColor','black');
+        pos=find(pressureload(:,end)==AirEnum());
+        pos=pressureload(pos,end-1);
+        A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); D=elements(pos,4); E=elements(pos,5); F=elements(pos,6);
+        h3=patch( 'Faces', [A B C], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','green','EdgeColor','black');
         patch( 'Faces', [D E F], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','green','EdgeColor','black');
         patch( 'Faces', [A B E D], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','green','EdgeColor','black');
         patch( 'Faces', [B E F C ], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','green','EdgeColor','black');
         patch( 'Faces', [C A D F ], 'Vertices', [x y z],'FaceVertexCData', [1 1 1],'FaceColor','green','EdgeColor','black');
-        hold on;
 
         %display arrows pointing outward
@@ -73 +85 @@
 %legend (disable warnings)
 warning off
-legend([h2,q],'element on ice front','normal vectors')
+legend([h2,h3,q],'element on ice front (Water)','element on ice front (Air)','normal vectors')
 warning on