Index: /issm/trunk/src/m/classes/public/plot/plot_pressureload.m
===================================================================
--- /issm/trunk/src/m/classes/public/plot/plot_pressureload.m	(revision 3102)
+++ /issm/trunk/src/m/classes/public/plot/plot_pressureload.m	(revision 3103)
@@ -21,9 +21,13 @@
 	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,5 +40,5 @@
 	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,13 +52,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); 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,5 +85,5 @@
 %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