Index: /issm/trunk/src/m/utils/Mesh/GetAreas.m
===================================================================
--- /issm/trunk/src/m/utils/Mesh/GetAreas.m	(revision 1794)
+++ /issm/trunk/src/m/utils/Mesh/GetAreas.m	(revision 1795)
@@ -16,4 +16,5 @@
 nels=size(index,1);
 nods=length(x);
+if nargin==4, z=varargin{1}; end
 
 %some checks
@@ -46,5 +47,5 @@
 else
 	%V=area(triangle)*1/3(z1+z2+z3)
-	thickness=mean(z(index(:,4:6)))-mean(z(index(:,1:3)));
-	areas=abs(0.5*((x2-x1).*(y3-y1)-(y2-y1).*(x3-x1)))*thickness;
+	thickness=mean(z(index(:,4:6)),2)-mean(z(index(:,1:3)),2);
+	areas=abs(0.5*((x2-x1).*(y3-y1)-(y2-y1).*(x3-x1))).*thickness;
 end