DepthAverage

PURPOSE ^

DEPTHAVERAGE - computes depth average of 3d vector, and return value on 2d mesh.

SYNOPSIS ^

function vector_average=DepthAverage(md,vector);

DESCRIPTION ^

DEPTHAVERAGE - computes depth average of 3d vector, and return value on 2d mesh. 

   Usage:
      vector_average=DepthAverage(md,vector);

   Example:
      vel_bar=DepthAverage(md,md.vel);

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function  vector_average=DepthAverage(md,vector);
0002 %DEPTHAVERAGE - computes depth average of 3d vector, and return value on 2d mesh.
0003 %
0004 %   Usage:
0005 %      vector_average=DepthAverage(md,vector);
0006 %
0007 %   Example:
0008 %      vel_bar=DepthAverage(md,md.vel);
0009 
0010 %check that the model given in input is 3d
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);

Generated on Sun 29-Mar-2009 20:22:55 by m2html © 2003