function analyze_this2(info, dirs)

test_name = info.run(1:2);

% ---
dummy = info;
dummy.run = test_name;
try
	analyze_this(dummy, dirs)
catch
end





%% read xbeach
fname            = [dirs.output filesep 'xboutput.nc'];
zb               = nc_varget(fname, 'zb'); 
globalx          = nc_varget(fname, 'globalx'); 

tsglobal         = squeeze(ncread(fname,'globaltime'));
%% read measurements

data_path_exp = [dirs.data filesep test_name];
filez = dir([data_path_exp filesep '*.tek']);


for ifile = 1:length(filez)
    tmp(ifile) = readTekalFile([data_path_exp filesep filez(ifile).name]);
    plot(tmp(ifile).data(:,1),tmp(ifile).data(:,2)); hold on
    if ifile==1
        X = tmp(ifile).data(:,1);
    end
    
end

figure;
bed         = [];
z0 = tmp(1).data(:,2);
for ifile = 1:length(filez)
    x = tmp(ifile).data(:,1);
    z = tmp(ifile).data(:,2);

    z = interp1(x,z,X);
    z(isnan(z)) = z0(isnan(z));
    plot(X,z); hold on
    
    bed = [bed z];
end
x_data = X;
z_data = bed;

%% analyse
zb  = squeeze(zb(:,1,:))';

xlim = [200 300];
ylim = [3 7];
zs0 = 5;



[fig1 fig2] = xb_compare_profiles(globalx, zb,x_data, z_data, tsglobal,zs0, test_name, xlim ,ylim,[],true,true,false  );
xb_write_plot(fig1, dirs.output, [info.run '_profiles'])
xb_write_plot(fig2, dirs.output, [info.run '_errorTime'])

end
