Index: /issm/trunk/src/m/classes/public/SectionValues.m
===================================================================
--- /issm/trunk/src/m/classes/public/SectionValues.m	(revision 5899)
+++ /issm/trunk/src/m/classes/public/SectionValues.m	(revision 5900)
@@ -92,10 +92,10 @@
 
 	%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;