Changeset 7219

Timestamp:

01/28/11 12:33:46 (14 years ago)

Author:

Mathieu Morlighem

Message:

Added extrusion of elementconnectivity

File:

• issm/trunk/src/m/model/extrude.m (modified) (3 diffs)

issm/trunk/src/m/model/extrude.m ¶

@@ -111 +111 @@
 md.upperelements=upperelements;
 
-
-
 %Save old mesh 
 md.x2d=md.x;
@@ -196 +194 @@
 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 +200 @@
         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