function xb_hydro_velocity_analyse(globalx, ue_mean,urms_mean,u_var, zb, obs)

Urms_hf = squeeze( mean(urms_mean,1));

Urms_lf = squeeze( sqrt(mean(u_var,1)));

U_mean = squeeze(mean(ue_mean,1));
% ---

figure;
subplot('position',[0.1 0.6 0.6 0.3]);
plot(globalx, Urms_hf,'r-','linewidth',2)
hold on
plot(globalx, Urms_lf,'b-','linewidth',2)
plot(globalx, U_mean,'g-','linewidth',2)
plot(obs.xurms, obs.Urms_hf,'ko')
plot(obs.xurms, obs.Urms_lf,'k<')
plot(obs.xum, obs.Um,'ks')
grid on
xlabel('Distance [m]');ylabel('Velocity [m]'); title('Flow velocities')
legend({'Computed U_{rms,HF}','Computed U_{rms,LF}','Computed U_{mean}','Observed U_{rms,HF}','Observed U_{rms,LF}','Observed U_{mean}', },'Location',[.75 .6 .2 .3])

subplot('position',[0.1 0.1 0.6 0.3]);
plot(globalx, zb(:,1),'k--','linewidth',2)
hold on
plot(globalx, zb(:,end),'k-','linewidth',2)

grid on
xlabel('z [m]');ylabel('Distance [m]');
legend({'Initial profile','Final profile'},'Location',[.75 .1 .2 .3])



end
