Index: /issm/trunk/src/m/model/extrude.m
===================================================================
--- /issm/trunk/src/m/model/extrude.m	(revision 7218)
+++ /issm/trunk/src/m/model/extrude.m	(revision 7219)
@@ -111,6 +111,4 @@
 md.upperelements=upperelements;
 
-
-
 %Save old mesh 
 md.x2d=md.x;
@@ -196,6 +194,5 @@
 md.spcthickness=project3d(md,md.spcthickness,'node');
 
-%Extrusion of Neumann BC
-%in 3d, segmentonnumann is: [grid1 grid2 grid3 grid4 element]
+%in 3d, pressureload: [grid1 grid2 grid3 grid4 element]
 pressureload_layer1=[md.pressureload(:,1:2)  md.pressureload(:,2)+md.numberofgrids2d  md.pressureload(:,1)+md.numberofgrids2d  md.pressureload(:,3:4)]; %Add two columns on the first layer 
 pressureload=[];
@@ -203,7 +200,14 @@
 	pressureload=[pressureload ;pressureload_layer1(:,1:4)+(i-1)*md.numberofgrids2d pressureload_layer1(:,5)+(i-1)*md.numberofelements2d pressureload_layer1(:,6)];
 end
-
-%plug into md
 md.pressureload=pressureload;
+
+%connectivity
+md.elementconnectivity=repmat(md.elementconnectivity,numlayers-1,1);
+md.elementconnectivity(find(md.elementconnectivity==0))=NaN;
+for i=2:numlayers-1,
+	md.elementconnectivity((i-1)*md.numberofelements2d+1:(i)*md.numberofelements2d,:)...
+	=md.elementconnectivity((i-1)*md.numberofelements2d+1:(i)*md.numberofelements2d,:)+md.numberofelements2d;
+end
+md.elementconnectivity(find(isnan(md.elementconnectivity)))=0;
 
 %materials