Changeset 5241 for issm/trunk

Timestamp:

08/13/10 10:58:49 (15 years ago)

Author:

Eric.Larour

Message:

More comments, more robust

File:

• issm/trunk/src/m/classes/public/partition/AreaAverageOntoPartition.m (modified) (1 diff)

issm/trunk/src/m/classes/public/partition/AreaAverageOntoPartition.m

@@ -1 +1 @@
 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