Index: /issm/trunk/src/m/classes/public/SectionValues.m
===================================================================
--- /issm/trunk/src/m/classes/public/SectionValues.m	(revision 1170)
+++ /issm/trunk/src/m/classes/public/SectionValues.m	(revision 1171)
@@ -72,5 +72,5 @@
 
 	%Interpolation of data on specified points
-	data_interp=griddata_mesh_to_mesh(md.elements,md.x,md.y,data,X,Y);
+	data_interp=InterpFromMesh2d(md.elements,md.x,md.y,data,X,Y);
 
 	%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=griddata_mesh_to_mesh(md.elements2d,md.x2d,md.y2d,project2d(md,md.bed,1),X,Y,'node')+offset;
-	surface=griddata_mesh_to_mesh(md.elements2d,md.x2d,md.y2d,project2d(md,md.surface,1),X,Y,'node')+offset;
+	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;
 
 	%Some useful parameters
@@ -109,5 +109,5 @@
 
 	%Interpolation of data on specified points
-	data_interp=griddata_mesh_to_mesh_3d(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);
 
 	%build outputs
Index: /issm/trunk/src/m/classes/public/plot/plot_section.m
===================================================================
--- /issm/trunk/src/m/classes/public/plot/plot_section.m	(revision 1170)
+++ /issm/trunk/src/m/classes/public/plot/plot_section.m	(revision 1171)
@@ -19,5 +19,4 @@
 	index1=i;
 end
-
 
 %process data and model
Index: /issm/trunk/src/m/classes/public/plot/processdata.m
===================================================================
--- /issm/trunk/src/m/classes/public/plot/processdata.m	(revision 1170)
+++ /issm/trunk/src/m/classes/public/plot/processdata.m	(revision 1171)
@@ -11,4 +11,7 @@
 	error('plotmodel error message: data provided is empty');
 end
+
+%make a column out of data
+data=data(:);
 
 %check length