function analyze_this_nonh(info,dirs)
idtime      = [100    1600];

if strcmp(info.run(1:2),'1A')
    Tp = 2.1;
elseif strcmp(info.run(1:2),'1B')
    Tp = 2.1;
elseif strcmp(info.run(1:2),'1C')
    Tp = 3.4;
end

%% read xb
fname            = fullfile(dirs.output, 'xboutput.nc');
point_zs         = nc_varget(fname, 'point_zs'); 
globalx          = nc_varget(fname, 'globalx'); 
t                = nc_varget(fname, 'globaltime');
zb               = nc_varget(fname, 'zb'); zb = squeeze(zb(1,:,:));
zs               = nc_varget(fname, 'zs'); 
pointtime        = nc_varget(fname, 'pointtime');
pointx           = nc_varget(fname, 'pointx');

idwanted = pointtime>idtime(1) & pointtime<idtime(2);

result = xb_analyse_nonh(pointtime(idwanted), point_zs(idwanted,:),pointx,(1/Tp)/2, true);

%%
data = load([dirs.data filesep info.run(1:2) '.txt']);

datatime    = data(2:end,1);
data_zs     = data(2:end,2:end);

zs0_data = mean(data_zs(1:50,:));

data_zs = data_zs-zs0_data;

idwanted = datatime>idtime(1) & datatime<idtime(2);

result_data = xb_analyse_nonh(datatime(idwanted), data_zs(idwanted,:),pointx,(1/Tp)/2, true);

%%
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')

subplot(4,2,2)
plot(pointx, [result_data(:).Hm0_tot],'ko','linewidth',2);hold on
plot(pointx, [result(:).Hm0_tot],'r','linewidth',2)
grid on; title('$H_{m0}$ [m]','interpreter','latex')
ylim([0 0.3])
legend('Observed','xb')
subplot(4,2,3)
plot(pointx, [result_data(:).Hm0_LF],'ko','linewidth',2);hold on
plot(pointx, [result(:).Hm0_LF],'r','linewidth',2)
grid on; title('$H_{m0,LF}$ [m]','interpreter','latex')
ylim([0 0.06])
subplot(4,2,4)
plot(pointx, [result_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])
subplot(4,2,5)
plot(pointx, [result_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])
subplot(4,2,6)
plot(pointx, [result_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])
subplot(4,2,7)
plot(pointx, [result_data(:).skewness],'ko','linewidth',2);hold on
plot(pointx, [result(:).skewness],'r','linewidth',2)
grid on; title('skewness [-]','interpreter','latex')
xlabel('x [m]')
ylim([0 1.5])
subplot(4,2,8)
plot(pointx, [result_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 1])
xb_write_plot(gcf, dirs.output, 'hydro');


result_data = rmfield(result_data,{'fx','vy','A','H','groupiness'}) 



save(['data' info.run(1:2)],'result_data','pointx')
save(['xb' info.run(1:2)],'result','pointx')
end