Index: /issm/trunk/src/m/classes/public/modelextract.m
===================================================================
--- /issm/trunk/src/m/classes/public/modelextract.m	(revision 4347)
+++ /issm/trunk/src/m/classes/public/modelextract.m	(revision 4348)
@@ -68,5 +68,5 @@
 elements_2(:,2)=Pgrid(elements_2(:,2));
 elements_2(:,3)=Pgrid(elements_2(:,3));
-if strcmpi(md1.type,'3d'),
+if md1.dim==3,
 	elements_2(:,4)=Pgrid(elements_2(:,4));
 	elements_2(:,5)=Pgrid(elements_2(:,5));
@@ -107,5 +107,5 @@
 
 	%uppernodes lowernodes
-	if strcmpi(md1.type,'3d')
+	if md1.dim==3
 		md2.uppergrids=md1.uppergrids(pos_grid);
 		pos=find(~isnan(md2.uppergrids));
@@ -118,5 +118,5 @@
 
 	%Initial 2d mesh 
-	if strcmpi(md1.type,'3d')
+	if md1.dim==3
 		flag_elem_2d=flag_elem(1:md1.numberofelements2d);
 		pos_elem_2d=find(flag_elem_2d);
@@ -162,5 +162,5 @@
 
 	%recreate segments
-	if strcmpi(md1.type,'2d')
+	if md1.dim==2
 		md2.nodeconnectivity=NodeConnectivity(md2.elements,md2.numberofgrids);
 		md2.elementconnectivity=ElementConnectivity(md2.elements,md2.nodeconnectivity);
@@ -215,5 +215,5 @@
 		md2.pressureload(:,2)=Pgrid(md1.pressureload(:,2)); 
 		md2.pressureload(:,end-1)=Pelem(md1.pressureload(:,end-1));
-		if strcmpi(md1.type,'3d')
+		if md1.dim==3
 			md2.pressureload(:,3)=Pgrid(md1.pressureload(:,3)); 
 			md2.pressureload(:,4)=Pgrid(md1.pressureload(:,4));