[1] | 1 | function vector_average=DepthAverage(md,vector);
|
---|
| 2 | %DEPTHAVERAGE - computes depth average of 3d vector, and return value on 2d mesh.
|
---|
| 3 | %
|
---|
| 4 | % Usage:
|
---|
| 5 | % vector_average=DepthAverage(md,vector);
|
---|
| 6 | %
|
---|
| 7 | % Example:
|
---|
[9684] | 8 | % vel_bar=DepthAverage(md,md.initialization.vel);
|
---|
[1] | 9 |
|
---|
| 10 | %check that the model given in input is 3d
|
---|
[9719] | 11 | if ~md.mesh.dimension==3;
|
---|
[1] | 12 | error('DepthAverage error message: the model given in input must be 3d')
|
---|
| 13 | end
|
---|
| 14 |
|
---|
[2557] | 15 | %nods data
|
---|
[9725] | 16 | if (length(vector)==md.mesh.numberofvertices),
|
---|
| 17 | vector_average=zeros(md.mesh.numberofvertices2d,1);
|
---|
| 18 | for i=1:md.mesh.numberoflayers-1,
|
---|
[9734] | 19 | vector_average=vector_average+(project2d(md,vector,i)+project2d(md,vector,i+1))/2.*(project2d(md,md.mesh.z,i+1)-project2d(md,md.mesh.z,i));
|
---|
[2557] | 20 | end
|
---|
[9691] | 21 | vector_average=vector_average./project2d(md,md.geometry.thickness,1);
|
---|
[2557] | 22 |
|
---|
| 23 | %element data
|
---|
[9725] | 24 | elseif (length(vector)==md.mesh.numberofelements),
|
---|
| 25 | vector_average=zeros(md.mesh.numberofelements2d,1);
|
---|
| 26 | for i=1:md.mesh.numberoflayers-1,
|
---|
[9734] | 27 | vector_average=vector_average+project2d(md,vector,i).*(project2d(md,md.mesh.z,i+1)-project2d(md,md.mesh.z,i));
|
---|
[2557] | 28 | end
|
---|
[9691] | 29 | vector_average=vector_average./project2d(md,md.geometry.thickness,1);
|
---|
[2557] | 30 |
|
---|
| 31 | else
|
---|
| 32 | error('vector size not supported yet');
|
---|
[1] | 33 | end
|
---|