Changeset 512

Timestamp:

05/20/09 08:05:00 (16 years ago)

Author:

Mathieu Morlighem

Message:

fixed extrusion (second layer was skipped)

File:

• issm/trunk/src/m/classes/public/extrude.m (modified) (1 diff)

issm/trunk/src/m/classes/public/extrude.m

@@ -48 +48 @@
 thickness3d=md.thickness; %thickness and bed for these grids
 bed3d=md.bed;
+
 %Create the new layers
-for i=2:numlayers,
-        x3d=[x3d; md.x]; %build the grids of the other layers
+for i=1:numlayers-1,
+        x3d=[x3d; md.x]; 
         y3d=[y3d; md.y];
-        z3d=[z3d; bed3d+thickness3d*(i/numlayers)^extrusionexponent]; %grids are distributed between bed and surface accordingly to the given exponent
+        %grids are distributed between bed and surface accordingly to the given exponent
+        z3d=[z3d; bed3d+thickness3d*(i/(numlayers-1))^extrusionexponent]; 
 end
 number_grids3d=size(x3d,1); %number of 3d grids for the non extruded part of the mesh