Index: /issm/trunk/src/m/classes/mesh.m
===================================================================
--- /issm/trunk/src/m/classes/mesh.m	(revision 9734)
+++ /issm/trunk/src/m/classes/mesh.m	(revision 9735)
@@ -6,38 +6,43 @@
 classdef mesh
 	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);
-		vertexonboundary            = modelfield('default',NaN,'marshall',false);
+		x                 = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1);
+		y                 = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1);
+		z                 = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1);
+		elements          = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',2);
+		dimension         = modelfield('default',0,'marshall',true,'format','Integer');
+		numberoflayers    = modelfield('default',0,'marshall',true,'format','Integer');
+		numberofelements  = modelfield('default',0,'marshall',true,'format','Integer');
+		numberofvertices  = modelfield('default',0,'marshall',true,'format','Integer');
+		numberofedges     = modelfield('default',0,'marshall',true,'format','Integer');
+		
 		lat        = modelfield('default',NaN,'marshall',false);
 		long       = modelfield('default',NaN,'marshall',false);
 		hemisphere = modelfield('default',NaN,'marshall',false);
+
+		elementonbed      = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',2);
+		elementonsurface  = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',2);
+		vertexonbed       = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',1);
+		vertexonsurface   = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',1);
+		lowerelements     = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',2);
+		lowervertex       = modelfield('default',NaN,'marshall',false);
+		upperelements     = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',2);
+		uppervertex       = modelfield('default',NaN,'marshall',false);
+		vertexonboundary  = modelfield('default',NaN,'marshall',false);
+
+		edges               = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',3);
 		segments            = modelfield('default',NaN,'marshall',false);
 		segmentmarkers      = modelfield('default',NaN,'marshall',false);
-		dimension           = modelfield('default',0,'marshall',true,'format','Integer');
-		numberofelements    = modelfield('default',0,'marshall',true,'format','Integer');
-		numberofvertices    = modelfield('default',0,'marshall',true,'format','Integer');
-		numberofedges       = modelfield('default',0,'marshall',true,'format','Integer');
-		numberoflayers      = modelfield('default',0,'marshall',true,'format','Integer');
-		numberofvertices2d    = modelfield('default',0,'marshall',true,'format','Integer');
-		numberofelements2d = modelfield('default',0,'marshall',true,'format','Integer');
-		vertexconnectivity    = modelfield('default',NaN,'marshall',false);
+		vertexconnectivity  = modelfield('default',NaN,'marshall',false);
 		elementconnectivity = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',3);
-		uppervertex           = modelfield('default',NaN,'marshall',false);
-		upperelements        = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',2);
-		lowerelements        = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',2);
-		lowervertex           = modelfield('default',NaN,'marshall',false);
-		elementonbed      = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',2);
-		elementonsurface  = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',2);
-		vertexonbed         = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',1);
-		vertexonsurface     = modelfield('default',NaN,'marshall',true,'format','BooleanMat','mattype',1);
-		elements2d         = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',3);
+		average_vertex_connectivity = modelfield('default',0,'marshall',true,'format','Integer');
+
 		x2d                = modelfield('default',NaN,'marshall',false);
 		y2d                = modelfield('default',NaN,'marshall',false);
-		elements            = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',2);
-		edges               = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',3);
-		x                   = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1);
-		y                   = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1);
-		z                   = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',1);
+		elements2d         = modelfield('default',NaN,'marshall',true,'format','DoubleMat','mattype',3);
+		numberofvertices2d = modelfield('default',0,'marshall',true,'format','Integer');
+		numberofelements2d = modelfield('default',0,'marshall',true,'format','Integer');
+
+		extractedvertices  = modelfield('default',NaN,'marshall',false);
+		extractedelements  = modelfield('default',NaN,'marshall',false);
 	end
 	methods