Changeset 2557

Timestamp:

10/28/09 11:47:27 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added Element size support of DepthAverage.m

File:

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

issm/trunk/src/m/classes/public/DepthAverage.m

@@ -13 +13 @@
 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);