function analyse_this(info, dirs)

data_folder = dirs.data;
output_folder = dirs.output;

sim_str = info.run;
match = regexp(sim_str, '\d+', 'match');
sim = str2double(match{1});
name2 = "_nh+2DV";

Analyse_Buckly(sim,name2,output_folder, data_folder);

ppath = dirs.output;
data = load([ppath '/data']);
result = load([ppath '/result']);
model = load([ppath '/model']);

data = data.data;
result = result.result;
model = model.model;

globalx = model.globalx;
zb = model.zb;
zs = model.zs;
pointx = model.pointx;

%%
figure('pos',[100 100 500 600]);
subplot(4,2,1)
plot(globalx, zb,'k','linewidth',2)
hold on
plot(globalx, squeeze( zs(100,1,:)),'b','linewidth',2)
grid on; title('elevation [m]','interpreter','latex')
%ylim([-0.2,0.2])
xlim([7.37,40])

subplot(4,2,2)
plot(pointx, [data(:).Hs],'ko','linewidth',2);hold on
plot(pointx, [result(:).Hs],'r','linewidth',2)
grid on; title('$H_{m0}$ [m]','interpreter','latex')
xlim([7.37,40])
legend('Observed','XBeach')

subplot(4,2,3)
plot(pointx, [data(:).Hs_low],'ko','linewidth',2);hold on
plot(pointx, [result(:).Hs_low],'r','linewidth',2)
grid on; title('$H_{m0,LF}$ [m]','interpreter','latex')
%ylim([0 0.06])
xlim([7.37,40])

subplot(4,2,4)
plot(pointx, [data(:).eta],'ko','linewidth',2);hold on
plot(pointx, [result(:).eta]-mean([result(:).eta]),'r','linewidth',2)
grid on; title('$\eta$ [m]','interpreter','latex')
%ylim([-0.01 0.02])
xlim([7.37,40])

subplot(4,2,5)
plot(pointx, [data(:).Tm_10],'ko','linewidth',2);hold on
plot(pointx, [result(:).Tm_10],'r','linewidth',2)
grid on; title('$T_{m-1,0}$ [s]','interpreter','latex')
%ylim([0 10])
xlim([7.37,40])

subplot(4,2,6)
plot(pointx, [data(:).Tm02],'ko','linewidth',2);hold on
plot(pointx, [result(:).Tm02],'r','linewidth',2)
grid on; title('$T_{m02}$ [s]','interpreter','latex')
%ylim([1 3])
xlim([7.37,40])

subplot(4,2,7)
plot(pointx, [data(:).skew],'ko','linewidth',2);hold on
plot(pointx, [result(:).skew],'r','linewidth',2)
grid on; title('skewness [-]','interpreter','latex')
xlabel('x [m]')
%ylim([0 1.5])
xlim([7.37,40])

subplot(4,2,8)
plot(pointx, [data(:).asymmetry],'ko','linewidth',2);hold on
plot(pointx, [result(:).asymmetry],'r','linewidth',2)
grid on; title('asymmetry [-]','interpreter','latex')
xlabel('x [m]')
%ylim([-1.5 1])
xlim([7.37,40])

xb_write_plot(gcf, dirs.output, 'hydro');


end