function PlotDustMassFlux(Struct, N)
% PlotDustMassFlux(Struct, i)
% plot mass flux (kg/m^2/s) of the N^{th} bin of the input Struct, 
% as a function of shear velocity. Plot mass flux in ug/m^2/s

x = Struct.ustar;
y = Struct.Fd(:, N)*10^9; % unit: ug/m^2/s
y_err = Struct.Fd_err(:, N)*10^9; % unit: ug/m^2/s

plot(x, y,'o','MarkerFaceColor','none','MarkerEdgeColor','k','Markersize',7,'linewidth',1); hold on;
%--- plot error bar
for i = 1:length(x)
    tempx = [x(i), x(i)];
    tempy = [max(y(i)-y_err(i),0), y(i)+y_err(i)];
    plot(tempx,tempy,'k-'); hold on;
end

%--- plot setting
set(gca,'xscale','log');
set(gca,'yscale','log');
title(['Bin - ',num2str(N)],'FontSize',12,'FontWeight','bold');  
xlabel('Shear velocity \itu_{*}\rm\bf, m/s','FontSize',12,'FontWeight','bold');
ylabel('Vertical dust mass flux \itF_{d}\rm\bf, \mug/m^{2}/s','FontSize',12,'FontWeight','bold');

end