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'); 
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'); 
globalx         = globalx+1
ue_mean           = nc_varget(fname, 'ue_mean'); 
u_var             = nc_varget(fname, 'u_var'); 

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

obs       = struct();
f       = dir(fullfile(dirs.data, '*.asc'));
for j = 1:length(f)
    var     = f(j).name(1:end-4);
    obs.(var) = load(fullfile(dirs.data, f(j).name));

end
obs.Hrms_hf = obs.HRMS_H_I(:,2);
obs.Hrms_lf = obs.HRMS_L_O(:,2);
obs.xhrms = obs.HRMS_L_O(:,1);
obs.setup = obs.ETA(:,2)-mean(obs.ETA(:,2));
obs.xsetup = obs.ETA(:,1);
else
load([dirs.data filesep test_name '.mat']);


obs.Hrms_hf = s.Hrms_hf(:,2);
obs.Hrms_lf = s.Hrms_lf(:,2);
obs.xhrms = s.xps;
obs.setup = s.setup;
obs.xsetup = s.xps;

z_data = [s.z{1}' ];
x_data = [s.x{1}' ];
tsglobal = [tsglobal(end)];
end


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

zs0 = 4.56; % checken?


[fig1 fig2] = xb_compare_profiles(globalx, zb,x_data, z_data, tsglobal,zs0, test_name, [160 200] ,[3 7],[],true,true,true  );
xb_write_plot(fig1, dirs.output, [info.run '_profiles'])
xb_write_plot(fig2, dirs.output, [info.run '_errorTime'])

% --- 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-mean(zs_mean),zs_var,H_mean,zb,obs)
xb_write_plot(gcf, dirs.output, [test_name '_hydro'])


end
