Index: /issm/trunk/src/m/classes/public/SectionValues.m
===================================================================
--- /issm/trunk/src/m/classes/public/SectionValues.m	(revision 1180)
+++ /issm/trunk/src/m/classes/public/SectionValues.m	(revision 1181)
@@ -72,5 +72,5 @@
 
 	%Interpolation of data on specified points
-	data_interp=InterpFromMesh2d(md.elements,md.x,md.y,data,X,Y);
+	data_interp=InterpFromMesh2d(md.elements,md.x,md.y,data,X,Y,NaN);
 
 	%Compute index
@@ -83,6 +83,6 @@
 	%Get bed and surface for each 2d point, offset to make sure that it is inside the glacier system
 	offset=10^-10;
-	bed=InterpFromMesh2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.bed,1),X,Y)+offset;
-	surface=InterpFromMesh2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.surface,1),X,Y)-offset;
+	bed=InterpFromMesh2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.bed,1),X,Y,NaN)+offset;
+	surface=InterpFromMesh2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.surface,1),X,Y,NaN)-offset;
 
 	%Some useful parameters
@@ -109,5 +109,5 @@
 
 	%Interpolation of data on specified points
-	data_interp=InterpFromMesh3d(md.elements,md.x,md.y,md.z,data,X3,Y3,Z3);
+	data_interp=InterpFromMesh3d(md.elements,md.x,md.y,md.z,data,X3,Y3,Z3,NaN);
 
 	%build outputs