function xb_plot_concentration(globalx, ccg_mean,obs)

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

% ---

figure;
subplot('position',[0.1 0.1 0.6 0.8]);
plot(globalx, ccg_mean*1000,'r-','linewidth',2)
hold on
plot(obs.xum, obs.Cm,'ko')
grid on
xlabel('Distance [m]');ylabel('Concentration [g/l]'); title('Sediment concentration')
legend({'Computed','Observed '},'Location',[.75 .1 .2 .3])




end
