Index: /issm/trunk/src/m/classes/public/partition/AreaAverageOntoPartition.m
===================================================================
--- /issm/trunk/src/m/classes/public/partition/AreaAverageOntoPartition.m	(revision 5240)
+++ /issm/trunk/src/m/classes/public/partition/AreaAverageOntoPartition.m	(revision 5241)
@@ -1,10 +1,18 @@
 function partvector=AreaAverageOntoPartition(md,vector)
 %AREAAVERAGEONTOPARTITION  compute partition values for a certain vector expressed on the vertices of the mesh. Use area weighted average.
+%
+%   Usage: average=AreaAverageOntoPartition(md,vector)
+%
 
-partvector=zeros(max(md.part),1);
+%ok, first check that part is Matlab matlab indexed
+part=md.part+1;
+
+%initialize output
+partvector=zeros(max(part),1);
+
+%start weight average
 weightedvector=vector.*md.vwgt;
-
-for i=1:max(md.part),
-	pos=find(md.part==i);
+for i=1:max(part),
+	pos=find(part==i);
 	partvector(i)=sum(weightedvector(pos))/sum(md.vwgt(pos));
 end