Index: /issm/trunk/src/m/classes/public/partition/AreaAverageOntoPartition.m
===================================================================
--- /issm/trunk/src/m/classes/public/partition/AreaAverageOntoPartition.m	(revision 5216)
+++ /issm/trunk/src/m/classes/public/partition/AreaAverageOntoPartition.m	(revision 5216)
@@ -0,0 +1,10 @@
+function partvector=AreaAverageOntoPartition(md,vector)
+%AREAAVERAGEONTOPARTITION  compute partition values for a certain vector expressed on the vertices of the mesh. Use area weighted average.
+
+partvector=zeros(max(md.part),1);
+weightedvector=vector.*md.vwgt;
+
+for i=1:max(md.part),
+	pos=find(md.part==i);
+	partvector(i)=sum(weightedvector(pos))/sum(md.vwgt(pos));
+end