Index: /issm/trunk-jpl/src/m/classes/model.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.m	(revision 19956)
+++ /issm/trunk-jpl/src/m/classes/model.m	(revision 19957)
@@ -511,5 +511,5 @@
 				md2.mesh.vertexconnectivity=NodeConnectivity(md2.mesh.elements,md2.mesh.numberofvertices);
 				md2.mesh.elementconnectivity=ElementConnectivity(md2.mesh.elements,md2.mesh.vertexconnectivity);
-				md2.mesh.segments=contourenvelope(md2);
+				md2.mesh.segments=contourenvelope(md2.mesh);
 				md2.mesh.vertexonboundary=zeros(numberofvertices2,1); md2.mesh.vertexonboundary(md2.mesh.segments(:,1:2))=1;
 			else
@@ -517,5 +517,5 @@
 				md2.mesh.vertexconnectivity=NodeConnectivity(md2.mesh.elements2d,md2.mesh.numberofvertices2d);
 				md2.mesh.elementconnectivity=ElementConnectivity(md2.mesh.elements2d,md2.mesh.vertexconnectivity);
-				segments=contourenvelope(md2);
+				segments=contourenvelope(md2.mesh);
 				md2.mesh.vertexonboundary=zeros(numberofvertices2/md2.mesh.numberoflayers,1); md2.mesh.vertexonboundary(segments(:,1:2))=1;
 				md2.mesh.vertexonboundary=repmat(md2.mesh.vertexonboundary,md2.mesh.numberoflayers,1);
Index: /issm/trunk-jpl/src/m/contrib/massbalance/contourmassbalance.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/massbalance/contourmassbalance.m	(revision 19956)
+++ /issm/trunk-jpl/src/m/contrib/massbalance/contourmassbalance.m	(revision 19957)
@@ -18,5 +18,5 @@
 
 %Get segments enveloping contour
-segments=contourenvelope(md,file);
+segments=contourenvelope(md.mesh,file);
 %md.stressbalance.icefront=segments; plotmodel(md,'data','pressureload','expdisp',file);
 
Index: /issm/trunk-jpl/src/m/mesh/rifts/meshaddrifts.m
===================================================================
--- /issm/trunk-jpl/src/m/mesh/rifts/meshaddrifts.m	(revision 19956)
+++ /issm/trunk-jpl/src/m/mesh/rifts/meshaddrifts.m	(revision 19957)
@@ -34,5 +34,5 @@
 
 	%create domain of md2 model: 
-	md2.mesh.segments=contourenvelope(md2,'Meshaddrifts.Contour.exp');
+	md2.mesh.segments=contourenvelope(md2.mesh,'Meshaddrifts.Contour.exp');
 	domain_index=md2.mesh.segments(1,1:2);
 	while (domain_index(end)~=domain_index(1)),