Index: /issm/trunk/src/m/classes/mesh.m
===================================================================
--- /issm/trunk/src/m/classes/mesh.m	(revision 9705)
+++ /issm/trunk/src/m/classes/mesh.m	(revision 9706)
@@ -7,4 +7,6 @@
 	properties (SetAccess=public) 
 		average_vertex_connectivity = modelfield('default',0,'marshall',true,'format','Integer');
+		extractedvertices           = modelfield('default',NaN,'marshall',false);
+		extractedelements           = modelfield('default',NaN,'marshall',false);
 	end
 	methods
Index: /issm/trunk/src/m/classes/model/model.m
===================================================================
--- /issm/trunk/src/m/classes/model/model.m	(revision 9705)
+++ /issm/trunk/src/m/classes/model/model.m	(revision 9706)
@@ -84,8 +84,4 @@
 		 nodeonsurface     = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',1);
 
-		 %Extraction
-		 extractednodes    = modelfield('default',NaN,'marshall',false);
-		 extractedelements = modelfield('default',NaN,'marshall',false);
-
 		 %Boundary conditions
 		 nodeonboundary = modelfield('default',NaN,'marshall',false);
@@ -234,5 +230,5 @@
 			 if isfield(structmd,'gridonbed'), md.nodeonbed=structmd.gridonbed; end
 			 if isfield(structmd,'gridonsurface'), md.nodeonsurface=structmd.gridonsurface; end
-			 if isfield(structmd,'extractedgrids'), md.extractednodes=structmd.extractedgrids; end
+			 if isfield(structmd,'extractedgrids'), md.mesh.extractedvertices=structmd.extractedgrids; end
 			 if isfield(structmd,'gridoniceshelf'), md.mask.vertexonfloatingice=structmd.gridoniceshelf; end
 			 if isfield(structmd,'gridonicesheet'), md.mask.vertexongroundedice=structmd.gridonicesheet; end
@@ -330,4 +326,6 @@
 			 if isfield(structmd,'thickness_coeff'), md.geometry.hydrostatic_ratio=structmd.thickness_coeff; end
 			 if isfield(structmd,'connectivity'), md.mesh.average_vertex_connectivity=structmd.connectivity; end
+			 if isfield(structmd,'extractednodes'), md.mesh.extractedvertices=structmd.extractednodes; end
+			 if isfield(structmd,'extractedelements'), md.mesh.extractedelements=structmd.extractedelements; end
 
 			 %Field changes
Index: /issm/trunk/src/m/model/modelextract.m
===================================================================
--- /issm/trunk/src/m/model/modelextract.m	(revision 9705)
+++ /issm/trunk/src/m/model/modelextract.m	(revision 9706)
@@ -267,4 +267,4 @@
 
 %Keep track of pos_node and pos_elem
-md2.extractednodes=pos_node;
-md2.extractedelements=pos_elem;
+md2.mesh.extractedvertices=pos_node;
+md2.mesh.extractedelements=pos_elem;