%% Open data file
datafs1 = importdata('FigS1.csv');
timefs1 = duration(datafs1.textdata(2:end),"InputFormat","hh:mm");
diam_smps = datafs1.data(1,:);
dndlogdp_norm = datafs1.data(2:end,:);
%% Call new figure
figure('Position',[1063,340,473,347])
%% Plot the figure
imagesc(datenum(timefs1),diam_smps,dndlogdp_norm');
ax=gca;
set(ax,'YDir','normal','XTick',0:3:24,'YScale','log','box','off','YTick',[10 20 50 100 200 500]),ylim([10 635])
ylabel('Diameter (nm)')
h1=colorbar();
h1.Label.String = 'Normalized Concentration';
datetick(ax)
colormap(viridis);
xticklabels({'12:00','18:00','00:00','06:00','12:00'})
ax.XTickLabelRotation = 0;
xlabel('Local Time','FontSize',14)
title('Non Event Days - ZF2 Site');