Index: /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.m
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.m	(revision 15945)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.m	(revision 15946)
@@ -35,5 +35,5 @@
 
 %Ice front position: 
-md.mask.icelevelset(find(nodeonicefront))=0;
+md.mask.ice_levelset(find(nodeonicefront))=0;
 
 %First find segments that are not completely on the front
@@ -43,5 +43,5 @@
 	numbernodesfront=4;
 end
-segmentsfront=md.mask.icelevelset(md.mesh.segments(:,1:numbernodesfront))==0;
+segmentsfront=md.mask.ice_levelset(md.mesh.segments(:,1:numbernodesfront))==0;
 segments=find(sum(segmentsfront,2)~=numbernodesfront);
 %Find all nodes for these segments and spc them
Index: /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 15945)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py	(revision 15946)
@@ -40,5 +40,5 @@
 	#Icefront position
 	pos=numpy.nonzero(nodeonicefront)[0]
-	md.mask.icelevelset[pos]=0
+	md.mask.ice_levelset[pos]=0
 
 	#First find segments that are not completely on the front
@@ -47,6 +47,6 @@
 	else:
 		numbernodesfront=4
-	if any(md.mask.icelevelset<=0):
-		values=md.mask.icelevelset[md.mesh.segments[:,0:-1]-1]
+	if any(md.mask.ice_levelset<=0):
+		values=md.mask.ice_levelset[md.mesh.segments[:,0:-1]-1]
 		segmentsfront=1-values
 		numpy.sum(segmentsfront,axis=1)!=numbernodesfront
Index: /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.m
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.m	(revision 15945)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.m	(revision 15946)
@@ -42,5 +42,5 @@
 
 %Position of ice front
-md.mask.icelevelset(find(vertexonicefront))=0;
+md.mask.ice_levelset(find(vertexonicefront))=0;
 
 %First find segments that are not completely on the front
@@ -50,5 +50,5 @@
 	numbernodesfront=4;
 end
-segmentsfront=md.mask.icelevelset(md.mesh.segments(:,1:numbernodesfront))==0;
+segmentsfront=md.mask.ice_levelset(md.mesh.segments(:,1:numbernodesfront))==0;
 segments=find(sum(segmentsfront,2)~=numbernodesfront);
 %Find all nodes for these segments and spc them
Index: /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py
===================================================================
--- /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 15945)
+++ /issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.py	(revision 15946)
@@ -48,5 +48,5 @@
 	#Position of ice front
 	pos=numpy.nonzero(vertexonicefront)[0]
-	md.mask.icelevelset[pos]=0
+	md.mask.ice_levelset[pos]=0
 
 	#First find segments that are not completely on the front
@@ -55,6 +55,6 @@
 	else:
 		numbernodesfront=4
-	if any(md.mask.icelevelset<=0):
-		values=md.mask.icelevelset[md.mesh.segments[:,0:-1]-1]
+	if any(md.mask.ice_levelset<=0):
+		values=md.mask.ice_levelset[md.mesh.segments[:,0:-1]-1]
 		segmentsfront=1-values
 		numpy.sum(segmentsfront,axis=1)!=numbernodesfront
Index: /issm/trunk-jpl/src/m/plot/plot_icefront.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_icefront.m	(revision 15945)
+++ /issm/trunk-jpl/src/m/plot/plot_icefront.m	(revision 15946)
@@ -12,7 +12,7 @@
 %process mesh and data
 [x y z elements is2d isplanet]=processmesh(md,[],options);
-ice=(md.mask.icelevelset>0);
-noice=(md.mask.icelevelset<=0);
-zeroice=(md.mask.icelevelset==0);
+ice=(md.mask.ice_levelset>0);
+noice=(md.mask.ice_levelset<=0);
+zeroice=(md.mask.ice_levelset==0);
 elementice=sum(ice(md.mesh.elements),2);
 elementnoice=sum(noice(md.mesh.elements),2);
@@ -33,5 +33,5 @@
 	hold on;
 
-	%Plot zero icelevelset line
+	%Plot zero ice_levelset line
 
 else