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

if strcmp(info.run(1:3),'jsn')
    Tp = 1;
elseif strcmp(info.run(1:3),'jln')
    Tp = 1/0.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);

%%
startRow = 3;
formatSpec = '%8s%9s%7s%7s%7s%7s%7s%7s%s%[^\n\r]';
fileID      = fopen([dirs.data filesep info.run(1:3) '.dat'],'r');
data = textscan(fileID,'%f%f%f%f%f%f%f%f%f','HeaderLines',1, 'EndOfLine', '\r\n');
fclose(fileID);

x_loc = [6 11 12 13 14 15 16 17];
datatime = data{1,1};
whm(:,1) = data{1,2}/100;
whm(:,2) = data{1,3}/100;
whm(:,3) = data{1,4}/100;
whm(:,4) = data{1,5}/100;
whm(:,5) = data{1,6}/100;
whm(:,6) = data{1,7}/100;
whm(:,7) = data{1,8}/100;
whm(:,8) = data{1,9}/100;




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

result_data = xb_analyse_nonh(datatime(idwanted), whm(idwanted,:),x_loc,(1/Tp)/2, true);

%%
figure('pos',[100 100 600 700]);
subplot(4,2,1)
plot(globalx, zb,'k','linewidth',2)
hold on
plot(globalx, squeeze( zs(100,1,:)),'b','linewidth',2)
grid on; title('domain')

subplot(4,2,2)
plot( x_loc, [result_data(:).Hm0_tot],'bo','linewidth',2);hold on
plot(pointx, [result(:).Hm0_tot],'r','linewidth',2)
grid on; title('H_{m0}')
legend('Observed','xb','location','southwest')
subplot(4,2,3)
plot(x_loc, [result_data(:).Hm0_LF],'bo','linewidth',2);hold on
plot(pointx, [result(:).Hm0_LF],'r','linewidth',2)
grid on; title('H_{m0,LF}')

subplot(4,2,4)
plot(x_loc, [result_data(:).eta],'bo','linewidth',2);hold on
plot(pointx, [result(:).eta]-mean([result(:).eta]),'r','linewidth',2)
grid on; title('\eta')

subplot(4,2,5)
plot(x_loc, [result_data(:).Tm_10],'bo','linewidth',2);hold on
plot(pointx, [result(:).Tm_10],'r','linewidth',2)
grid on; title('T_{m-1,0}')

subplot(4,2,6)
plot(x_loc, [result_data(:).Tm02],'bo','linewidth',2);hold on
plot(pointx, [result(:).Tm02],'r','linewidth',2)
grid on; title('T_{m02}')

subplot(4,2,7)
plot(x_loc, [result_data(:).skewness],'bo','linewidth',2);hold on
plot(pointx, [result(:).skewness],'r','linewidth',2)
grid on; title('skewness')

subplot(4,2,8)
plot(x_loc, [result_data(:).asymmetry],'bo','linewidth',2);hold on
plot(pointx, [result(:).asymmetry],'r','linewidth',2)
grid on; title('asymmetry')
xb_write_plot(gcf, dirs.output, 'hydro');
end