function PlotPSDCompareAfricaBullard(Struct, OStruct, Bin, D72, dvdlnd72, D120, dvdlnd120)
% PlotPSDCompareAfrica(Struct, OStruct, Bin)
% Input Oceano normalized volume PSD "Struct", normalized volume PSD of
% African in situ measurements "OStruct" (under geometric diameter), the
% geometric diamter bin size "Bin" (unit: um). Plot figures of comparison
% of PSD

% plot PSD averaged over all shear velocity of Oceano Dunes
PlotPSDMean(Struct, Bin, [0,0,0]); % plot the mean PSD of Oceano over all u*

% plot African in situ PSDs
chou_plot = plot(OStruct.cho08.SizeUm,OStruct.cho08.dVdlnD_B160,'v',OStruct.cho08.SizeUm,OStruct.cho08.dVdlnD_B165,'v'); hold on; %Chou et al. (2008)
mcconnell_plot=plot(OStruct.mcc08.SizeUm1,OStruct.mcc08.dVdlnD1,'<',OStruct.mcc08.SizeUm2,OStruct.mcc08.dVdlnD2,'<'); hold on; %McConnell et al. (2008)
osborne_plot = plot(OStruct.osb08.SizeUm_B160, OStruct.osb08.dVdlnD_B160,'^',OStruct.osb08.SizeUm_B165,OStruct.osb08.dVdlnD_B165,'^'); hold on; %Osborne et al. (2008)
weinzierl_plot = plot(OStruct.wei09.SizeUm_L02,OStruct.wei09.dVdlnD_L02,'>',OStruct.wei09.SizeUm_L07,OStruct.wei09.dVdlnD_L07,'>'); hold on; %Weinzierl et al. (2009)
ryder_plot = plot(OStruct.ryd13.SizeUm,OStruct.ryd13.dVnorm,'s'); hold on; %Ryder et al. (2013)

% plot Bullard et al. (2004)'s 2 PSDs
Bu72 = plot(D72, dvdlnd72, 'x'); hold on;
Bu120 = plot(D120, dvdlnd120, 'd'); hold on;
 
set(chou_plot(1),'markeredge',[0.2,0.4,0],'markersize',6);
set(chou_plot(2),'markeredge',[0.2,0.4,0],'markersize',6);
set(mcconnell_plot(1),'markeredge',[1,0.4,0.2],'markersize',6);
set(mcconnell_plot(2),'markeredge',[1,0.4,0.2],'markersize',6);
set(osborne_plot(1),'markeredge',[0,0.6,0.8],'markersize',6);
set(osborne_plot(2),'markeredge',[0,0.6,0.8],'markersize',6);
set(weinzierl_plot(1),'markeredge',[0.8,0.6,0.4],'markersize',6);
set(weinzierl_plot(2),'markeredge',[0.8,0.6,0.4],'markersize',6);
set(ryder_plot,'markeredge',[0.4,0,0],'markersize',6);
set(Bu72,'Markerface','none','markeredge','m','markersize',6);
set(Bu120,'Markerface','none','markeredge','m','markersize',6);
 
%--prepare for the legend
Oceano = plot(-10,-10,'o','MarkerFaceColor','none','MarkerEdgeColor',[0,0,0],'Markersize',7,'linewidth',1); hold on;
%--- add Kok (2011) theory ---
Ds = linspace(0.1,20,1000);
Ds_avge = 3.4;
Sigma = 3.0;
lam = 12;
cv = 12.62;
JK_dvdlnd = Ds/cv.*(1+erf(log(Ds/Ds_avge)/sqrt(2)/log(Sigma))).*exp(-(Ds/lam).^3);
JK2011 = plot(Ds,JK_dvdlnd,'b-.'); hold on;
%--- set the legend ---
thelegend = legend([chou_plot(1) mcconnell_plot(1) osborne_plot(1) weinzierl_plot(1) ryder_plot Bu72 Bu120 Oceano JK2011],...
        'Chou et al. (2008)','McConnell et al. (2008)','Osborne et al. (2008)','Weinzierl et al. (2009)','Ryder et al. (2013)',...
        'Bullard et al. (2004) - 72h','Bullard et al. (2004) - 120h','Oceano Dunes','Brittle fragmentation theory (Kok, 2011)');
    
set(thelegend,'location','southeast','FontSize',8,'FontWeight','bold');
    
%--- plot setting---
set(gca,'xscale','log');
set(gca,'yscale','log');

axis([0.4,19,5*10^(-4),1.9]);
set(gca,'xtick',[0.4,1,10,19]);
set(gca,'xticklabel',{'0.4','1','10',''});
set(gca,'ytick',[10^(-3),10^(-2),10^(-1),10^(0)]);
set(gca,'yticklabel',{'10^{-3}','10^{-2}','10^{-1}','10^{0}'});

% xlabel('Dust geometric diameter \itD_{g,i}, \rm\bf\mum','FontSize',14,'FontWeight','bold');
% ylabel('Norm. volume size distr. \itdV_{i}/dlnD_{g,i}','FontSize',14,'FontWeight','bold');

xlabel('Dust geometric diameter \itD_{g}, \rm\bf\mum','FontSize',14,'FontWeight','bold');
ylabel('Norm. volume size distr. \itdV/dlnD_{g}','FontSize',14,'FontWeight','bold');

end