


compute depth averaged horizontal velocity


0001 function velocity_average_g=CieloHorizontalVelocityDepthAverage(md,velocity_g); 0002 %compute depth averaged horizontal velocity 0003 0004 vx_average2d=DepthAverage(md,velocity_g(1:6:uset.gsize)); 0005 vy_average2d=DepthAverage(md,velocity_g(2:6:uset.gsize)); 0006 0007 %project onto layer 1 of 3d mesh 0008 vx_average=project3d(md,vx_average2d,'node',1); 0009 vy_average=project3d(md,vy_average2d,'node',1); 0010 0011 %output on g-set 0012 velocity_average_g=zeros(uset.gsize,1);velocity_average_g(1:6:uset.gsize)=vx_average;velocity_average_g(2:6:uset.gsize)=vy_average;