Index: /issm/trunk/src/m/classes/public/display/displayparameters.m
===================================================================
--- /issm/trunk/src/m/classes/public/display/displayparameters.m	(revision 2529)
+++ /issm/trunk/src/m/classes/public/display/displayparameters.m	(revision 2530)
@@ -11,5 +11,5 @@
 disp(sprintf('   Parameters:'));
 
-disp(sprintf('\n      geometry:'));
+disp(sprintf('\n      geometrical parameters:'));
 fielddisplay(md,'surface','surface height [m]');
 fielddisplay(md,'thickness','thickness [m]');
Index: /issm/trunk/src/m/classes/public/modelextract.m
===================================================================
--- /issm/trunk/src/m/classes/public/modelextract.m	(revision 2529)
+++ /issm/trunk/src/m/classes/public/modelextract.m	(revision 2530)
@@ -179,5 +179,5 @@
 		md2.nodeconnectivity=NodeConnectivity(md2.elements,md2.numberofgrids);
 		md2.elementconnectivity=ElementConnectivity(md2.elements,md2.nodeconnectivity);
-		md2.segments=findsegments(md2);
+		md2.segments=contourenvelope(md2);
 		md2.gridonboundary=zeros(numberofgrids2,1); md2.gridonboundary(md2.segments(:,1:2))=1;
 	end
@@ -248,4 +248,5 @@
 	%Results fields
 	if isstruct(md1.results),
+		md2.results=struct();
 		solutionfields=fields(md1.results);
 		for i=1:length(solutionfields),
@@ -258,4 +259,6 @@
 				elseif length(field)==numberofelements1,
 					md2.results.(solutionfields(i)).(solutionsubfields(j))=field(pos_elem);
+				else
+					md2.results.(solutionfields(i)).(solutionsubfields(j))=field;
 				end
 			end