Index: /issm/trunk/src/m/classes/model.m
===================================================================
--- /issm/trunk/src/m/classes/model.m	(revision 9463)
+++ /issm/trunk/src/m/classes/model.m	(revision 9464)
@@ -449,8 +449,5 @@
 
 			 %Field class change
-			 if (size(structmd.pressureload,2)==3 | size(structmd.pressureload,2)==5),
-				 md.pressureload=[structmd.pressureload...
-					 WaterEnum*structmd.elementoniceshelf(structmd.pressureload(:,end))+AirEnum*structmd.elementonicesheet(structmd.pressureload(:,end))];
-			 end
+
 			 if (isfield(structmd,'type') & ischar(structmd.type)), 
 				 if strcmpi(structmd.type,'2d'), md.dim=2; end
@@ -475,4 +472,9 @@
 				 pos=find(structmd.spcvelocity(:,2)); md.spcvy(pos)=structmd.spcvelocity(pos,5); 
 				 pos=find(structmd.spcvelocity(:,3)); md.spcvz(pos)=structmd.spcvelocity(pos,6); 
+			 end
+			 if ismember(structmd.pressureload(end,end),[AirEnum() WaterEnum() IceEnum()]),
+				 pos=find(md.pressureload(:,end)==AirEnum);   md.pressureload(pos,end)=0;
+				 pos=find(md.pressureload(:,end)==WaterEnum); md.pressureload(pos,end)=1;
+				 pos=find(md.pressureload(:,end)==IceEnum);   md.pressureload(pos,end)=2;
 			 end
 
Index: /issm/trunk/src/m/model/display/displaybc.m
===================================================================
--- /issm/trunk/src/m/model/display/displaybc.m	(revision 9463)
+++ /issm/trunk/src/m/model/display/displaybc.m	(revision 9464)
@@ -22,5 +22,5 @@
 fielddisplay(md,'spcvy','y-axis velocity constraint (NaN means no constraint)');
 fielddisplay(md,'spcvz','z-axis velocity constraint (NaN means no constraint)');
-fielddisplay(md,'pressureload','segments on ice front list');
+fielddisplay(md,'pressureload','segments on ice front list (0:Air 1:Water 2:Ice)');
 
 disp(sprintf('\n      prognostic:'));
Index: /issm/trunk/src/m/model/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/model/ismodelselfconsistent.m	(revision 9463)
+++ /issm/trunk/src/m/model/ismodelselfconsistent.m	(revision 9464)
@@ -79,5 +79,5 @@
 	message('dim should be either 2 3');
 end
-checkvalues(md,{'pressureload(:,end)'},[WaterEnum() AirEnum() IceEnum()]);
+checkvalues(md,{'pressureload(:,end)'},[0 1 2]);
 %}}}
 %NO NAN {{{1
Index: /issm/trunk/src/m/model/plot/plot_pressureload.m
===================================================================
--- /issm/trunk/src/m/model/plot/plot_pressureload.m	(revision 9463)
+++ /issm/trunk/src/m/model/plot/plot_pressureload.m	(revision 9464)
@@ -22,14 +22,14 @@
 
 	%highlight elements on neumann
-	pos=find(pressureload(:,end)==WaterEnum());
+	pos=find(pressureload(:,end)==1);
 	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','blue','EdgeColor','black');
-	pos=find(pressureload(:,end)==AirEnum());
+	pos=find(pressureload(:,end)==0);
 	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');
 
-	pos=find(pressureload(:,end)==IceEnum());
+	pos=find(pressureload(:,end)==2);
 	pos=pressureload(pos,end-1);
 	A=elements(pos,1); B=elements(pos,2); C=elements(pos,3); 
@@ -58,5 +58,5 @@
 
 	%highlight elements on neumann
-	pos=find(pressureload(:,end)==WaterEnum());
+	pos=find(pressureload(:,end)==1);
 	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);
@@ -66,5 +66,5 @@
 	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=find(pressureload(:,end)==0);
 	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);
@@ -74,5 +74,5 @@
 	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');
-	pos=find(pressureload(:,end)==IceEnum());
+	pos=find(pressureload(:,end)==2);
 	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);
Index: /issm/trunk/src/m/utils/BC/SetIceShelfBC.m
===================================================================
--- /issm/trunk/src/m/utils/BC/SetIceShelfBC.m	(revision 9463)
+++ /issm/trunk/src/m/utils/BC/SetIceShelfBC.m	(revision 9464)
@@ -59,5 +59,5 @@
 
 %Add water or air enum depending on the element
-pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))];
+pressureload=[pressureload 1*md.elementoniceshelf(pressureload(:,end))];
 
 %plug onto model
Index: /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m
===================================================================
--- /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 9463)
+++ /issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m	(revision 9464)
@@ -69,5 +69,5 @@
 
 %Add water or air enum depending on the element
-pressureload=[pressureload WaterEnum*md.elementoniceshelf(pressureload(:,end))+AirEnum*md.elementonicesheet(pressureload(:,end))];
+pressureload=[pressureload 1*md.elementoniceshelf(pressureload(:,end))+ 0*md.elementonicesheet(pressureload(:,end))];
 
 %plug onto model