function analyze_this(info, dirs)

% initialize variables
x_data 	= struct([]);
zb_data = struct([]);

% obtain required output times (dependant on measurement times)
bathy_files = dir(fullfile([dirs.data filesep info.run filesep '*.tek']));

for il = 1:length(bathy_files)
    bathy = readTekalFile([dirs.data filesep info.run filesep bathy_files(il).name]);
    
    x_data{il} = bathy.data(:,1);
    zb_data{il} = bathy.data(:,2);
end


% interpolate all measurements to a single grid
for in = 1:length(zb_data)
    zb_data{in} = interp1(x_data{in},zb_data{in},x_data{end});
end

% read output
xbo = xb_read_output(pwd,'vars',{'globalx', 'globaly', 'globaltime', 'zs', 'zb', 'H'});

xb_plot_profile(xbo,'measured',[x_data{end} zb_data{:}])
title(['profiles ' info.test ' ' info.run], 'interpreter', 'none')
xb_write_plot(gcf, dirs.output, 'profiles')

clear bathy x_data zb_data
end     %end of function