0001 function vector_average=DepthAverage(md,vector);
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011 if ~strcmpi(md.type,'3d');
0012 error('DepthAverage error message: the model given in input must be 3d')
0013 end
0014
0015 vector_average=zeros(md.numberofgrids2d,1);
0016 for i=1:md.numlayers-1,
0017 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));
0018 end
0019 vector_average=vector_average./project2d(md,md.thickness,1);