Changeset 5900

Timestamp:

09/20/10 07:46:00 (14 years ago)

Author:

seroussi

Message:

minor

File:

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

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

@@ -92 +92 @@
 
         %Get bed and surface for each 2d point, offset to make sure that it is inside the glacier system
-        offset=10^-10;
+        offset=10^-3;
         bed=InterpFromMeshToMesh2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.bed,1),X,Y)+offset;
         surface=InterpFromMeshToMesh2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.surface,1),X,Y)-offset;
 
         %Some useful parameters
-        layers=ceil(res_v/md.numlayers);
+        layers=ceil(mean(md.thickness)/res_v);
         gridsperlayer=numberofgrids;
         gridstot=gridsperlayer*layers;