Index: /issm/trunk/src/m/classes/public/DepthAverage.m
===================================================================
--- /issm/trunk/src/m/classes/public/DepthAverage.m	(revision 2556)
+++ /issm/trunk/src/m/classes/public/DepthAverage.m	(revision 2557)
@@ -13,7 +13,21 @@
 end
 
-vector_average=zeros(md.numberofgrids2d,1);
-for i=1:md.numlayers-1,
-	vector_average=vector_average+(project2d(md,vector,i)+project2d(md,vector,i+1))/2.*(project2d(md,md.z,i+1)-project2d(md,md.z,i));
+%nods data
+if (length(vector)==md.numberofgrids),
+	vector_average=zeros(md.numberofgrids2d,1);
+	for i=1:md.numlayers-1,
+		vector_average=vector_average+(project2d(md,vector,i)+project2d(md,vector,i+1))/2.*(project2d(md,md.z,i+1)-project2d(md,md.z,i));
+	end
+	vector_average=vector_average./project2d(md,md.thickness,1);
+
+%element data
+elseif (length(vector)==md.numberofelements),
+	vector_average=zeros(md.numberofelements2d,1);
+	for i=1:md.numlayers-1,
+		vector_average=vector_average+project2d(md,vector,i).*(project2d(md,md.z,i+1)-project2d(md,md.z,i));
+	end
+	vector_average=vector_average./project2d(md,md.thickness,1);
+
+else
+	error('vector size not supported yet');
 end
-vector_average=vector_average./project2d(md,md.thickness,1);