function xb_hydro_shape_analyse(globalx, Sk_mean,As_mean,zb,obs)

Sk_mean = squeeze( mean(Sk_mean,1));

As_mean = squeeze( mean(As_mean,1));
% ---

figure;
subplot('position',[0.1 0.6 0.6 0.3]);
plot(globalx, Sk_mean,'r-','linewidth',2)
hold on
plot(globalx, As_mean,'b-','linewidth',2)
plot(obs.xhrms, obs.Sku','k<')
plot(obs.xhrms, obs.Asu','k>')
grid on
xlabel('Distance [m]');ylabel('Skewness & asymmetry [-]'); title('Wave shape')
legend({'Computed skewness','Computed asymmetry','Observed skewness','Observed asymmetry'},'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
