function PlotPSDComparison(Struct, OStruct, Bin) % PlotPSDComparison(Struct, OStruct, Bin) % Input Oceano normalized volume PSD "Struct", normalized volume PSD of % other non-sandy soils "OStruct" (under geometric diameter), the % geometric diamter bin size "Bin" (unit: um). Plot figures of comparison % of PSD PlotPSDMean(Struct, Bin, [0,0,0]); % plot the mean PSD of Oceano over all u* % plot other non-sandy PSDs Gil72 = plot(OStruct.Gil72etal.d,OStruct.Gil72etal.dvdlnd,'p'); hold on; %Gillette et al. (1972) Gil74a = plot(OStruct.Gil74.d,OStruct.Gil74.dvdlnd,'h'); hold on; %Gillette (1974) Gil74b = plot(OStruct.Gil74etal.d,OStruct.Gil74etal.dvdlnd,'s'); hold on; % Gillette et al. (1974) Frat = plot(OStruct.Fra07etal.d,OStruct.Fra07etal.dvdlnd,'>'); hold on;% Fratini et al. (2007) Sow = plot(OStruct.sow09etal.d,OStruct.sow09etal.dvdlnd,'^'); hold on;%Sow et al. (2009) Shao = plot(OStruct.sha11etal.d,OStruct.sha11etal.dvdlnd,'v'); hold on;% Shao et al. (2011) Ros = plot(OStruct.ros14etal.d,OStruct.ros14etal.dvdlnd,'<'); hold on;% Roseberg et al. (2014) set(Gil72,'Markerface','none','markeredge',[1,0.4,0.2],'markersize',6); set(Gil74a,'Markerface','none','markeredge',[1,0,1],'markersize',6); set(Gil74b,'Markerface','none','markeredge',[0.2,0.6,0.2],'markersize',6); set(Frat,'Markerface','none','markeredge',[0,0.6,0.8],'markersize',6); set(Sow,'Markerface','none','markeredge',[0,0.4,0],'markersize',6); set(Shao,'Markerface','none','markeredge',[0.4,0.4,0.8],'markersize',6); set(Ros,'Markerface','none','markeredge',[0.8,0.6,0.4],'markersize',6); %--- plot error in the PSDs for i = 1:4 tempx_Gil72 = [OStruct.Gil72etal.d(i),OStruct.Gil72etal.d(i)]; tempy_Gil72 = [OStruct.Gil72etal.dvdlnd(i)-OStruct.Gil72etal.dvdlnd_SE(i),OStruct.Gil72etal.dvdlnd(i)+OStruct.Gil72etal.dvdlnd_SE(i)]; plot(tempx_Gil72,tempy_Gil72,'-','color',[1,0.4,0.2]); hold on; tempx_Gil74a = [OStruct.Gil74.d(i),OStruct.Gil74.d(i)]; tempy_Gil74a = [OStruct.Gil74.dvdlnd(i)-OStruct.Gil74.dvdlnd_SE(i),OStruct.Gil74.dvdlnd(i)+OStruct.Gil74.dvdlnd_SE(i)]; plot(tempx_Gil74a,tempy_Gil74a,'-','color',[1,0,1]); hold on; tempx_Gil74b = [OStruct.Gil74etal.d(i),OStruct.Gil74etal.d(i)]; tempy_Gil74b = [OStruct.Gil74etal.dvdlnd(i)-OStruct.Gil74etal.dvdlnd_SE(i),OStruct.Gil74etal.dvdlnd(i)+OStruct.Gil74etal.dvdlnd_SE(i)]; plot(tempx_Gil74b,tempy_Gil74b,'-','color',[0.2,0.6,0.2]); hold on; end for i = 1:length(OStruct.Fra07etal.d) tempx_Frat = [OStruct.Fra07etal.d(i),OStruct.Fra07etal.d(i)]; tempy_Frat = [OStruct.Fra07etal.dvdlnd(i)-OStruct.Fra07etal.dvdlnd_SE(i),OStruct.Fra07etal.dvdlnd(i)+OStruct.Fra07etal.dvdlnd_SE(i)]; plot(tempx_Frat,tempy_Frat,'-','color',[0,0.6,0.8]); hold on; end for i = 1:length(OStruct.sow09etal.d) tempx_Sow = [OStruct.sow09etal.d(i),OStruct.sow09etal.d(i)]; tempy_Sow = [OStruct.sow09etal.dvdlnd(i)-OStruct.sow09etal.dvdlnd_SE(i),OStruct.sow09etal.dvdlnd(i)+OStruct.sow09etal.dvdlnd_SE(i)]; plot(tempx_Sow,tempy_Sow,'-','color',[0,0.4,0]); hold on; end for i = 1:length(OStruct.sha11etal.d) tempx_Shao = [OStruct.sha11etal.d(i),OStruct.sha11etal.d(i)]; tempy_Shao = [OStruct.sha11etal.dvdlnd(i)-OStruct.sha11etal.dvdlnd_SE(i),OStruct.sha11etal.dvdlnd(i)+OStruct.sha11etal.dvdlnd_SE(i)]; plot(tempx_Shao,tempy_Shao,'-','color',[0.4,0.4,0.8]); hold on; end for i = 1:length(OStruct.ros14etal.d) tempx_Ros = [OStruct.ros14etal.d(i),OStruct.ros14etal.d(i)]; tempy_Ros = [OStruct.ros14etal.dvdlnd(i)-OStruct.ros14etal.dvdlnd_SE(i),OStruct.ros14etal.dvdlnd(i)+OStruct.ros14etal.dvdlnd_SE(i)]; plot(tempx_Ros,tempy_Ros,'-','color',[0.8,0.6,0.4]); hold on; end %--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 --- thethe = legend([Gil72 Gil74a Gil74b Frat Sow Shao Ros Oceano JK2011],{'Gillette et al. (1972)',... 'Gillette (1974)', 'Gillette et al. (1974)','Fratini et al. (2007)','Sow et al. (2009)','Shao et al. (2011)',... 'Rosenberg et al. (2014)','Oceano Dunes','Brittle fragmentation theory (Kok, 2011)'}); set(thethe,'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