function analyze_this2(info, dirs)


if strcmp(info.run(1:2),'DP')
    test_name = info.run;
else
    test_name = info.run(1:3);
end
    


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





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

Sk_mean           = nc_varget(fname, 'Sk_mean'); 
As_mean           = nc_varget(fname, 'As_mean'); 
ue_mean           = nc_varget(fname, 'ue_mean'); 
u_var             = nc_varget(fname, 'u_var'); 
urms_mean         = nc_varget(fname, 'urms_mean'); 
ccg_mean          = nc_varget(fname, 'ccg_mean'); 

tsglobal         = squeeze(ncread(fname,'globaltime'));
%% read measurements
data = load([dirs.data filesep test_name filesep test_name '.txt'] );
x_data = data(:,1);
z_data = data(:,2:end);

% --- mat file
if ~strcmp(test_name,'T08')
    obs = load(fullfile(dirs.data, [test_name 'data.mat']));
end
%% analyse
zb  = squeeze(zb(:,1,:))';


if strcmp(test_name,'T01')
    xlim = [195 225];
    ylim = [3 6.5];
	zb = zb + 4.5;
    zs0 = 4.465;
elseif strcmp(test_name,'T02')
    xlim = [195 225];
    ylim = [3 6.5];
	zb = zb + 4.5;
	zs0 = 4.47;
elseif strcmp(test_name,'T03')
    xlim = [195 225];
    ylim = [3 6.5];
	zb = zb + 4.5;
	zs0 = 4.47;
elseif strcmp(test_name,'T08')
    xlim = [180 240];
    ylim = [3 6.5];
	zs0 = 4.5;
elseif strcmp(test_name,'T04')
    xlim = [180 240];
    ylim = [-1.5 2];
	zs0 = 0;
elseif strcmp(test_name,'DP01')
    xlim = [180 240];
    ylim = [-1.5 2];
	zs0 = 0;
elseif strcmp(test_name,'DP02')
    xlim = [180 240];
    ylim = [-1.5 2];
	zs0 = 0;
end


% --- compare profiles
[fig1 fig2] = xb_compare_profiles(globalx, zb,x_data, z_data, tsglobal,zs0, test_name, xlim ,ylim,[],true,true,true);
xb_write_plot(fig1, dirs.output, [info.run '_profiles'])
xb_write_plot(fig2, dirs.output, [info.run '_errorTime'])
if ~strcmp(test_name,'T08')
    % --- create morpho plot
    [fig1, fig2] = xb_morpho_analyse(globalx, zb,x_data, z_data, tsglobal,zs0);
    xb_write_plot(fig1, dirs.output, [test_name '_morpho'])
    xb_write_plot(fig2, dirs.output, [test_name '_volumes'])

    % -- hydro
    xb_hydro_analyse(globalx, zs_mean,zs_var,H_mean,zb,obs)
    xb_write_plot(gcf, dirs.output, [test_name '_hydro'])

    % -- velocity
    xb_hydro_velocity_analyse(globalx, ue_mean,urms_mean,u_var,zb,obs)
    xb_write_plot(gcf, dirs.output, [test_name '_velocity'])


    % -- shape
    xb_hydro_shape_analyse(globalx, Sk_mean,As_mean,zb,obs)
    xb_write_plot(gcf, dirs.output, [test_name '_shape'])

    % -- concentration
    xb_plot_concentration(globalx, ccg_mean,obs)
    xb_write_plot(gcf, dirs.output, [test_name '_concentration'])
end
end
