0001 function velocity_average_g=HorizontalVelocityDepthAverage(md,velocity_g);
0002
0003
0004
0005
0006
0007 global gridset
0008
0009 vx_average2d=DepthAverage(md,velocity_g(1:6:gridset.gsize));
0010 vy_average2d=DepthAverage(md,velocity_g(2:6:gridset.gsize));
0011
0012
0013 vx_average=project3d(md,vx_average2d,'node',1);
0014 vy_average=project3d(md,vy_average2d,'node',1);
0015
0016
0017 velocity_average_g=zeros(gridset.gsize,1);velocity_average_g(1:6:gridset.gsize)=vx_average;velocity_average_g(2:6:gridset.gsize)=vy_average;
0018
0019