function combinedPlots(info,dirs)

files  = dir(dirs.output);
dirFlags = [files.isdir];
somnames = files(dirFlags);
somnames(ismember( {somnames.name}, {'.', '..'})) = [];  %remove . and ..
clear files

index_hydro = contains({somnames.name},'MS');
somnames_hydro = somnames(index_hydro);
somnames_morfo = somnames(~index_hydro);

%% Morfo table

VeroXB_all = [];
VeroMeas_all = [];
retreatXB_all = [];
retreatMeas_all = [];

for i=1:size(somnames_morfo,1)
    load([dirs.output filesep somnames_morfo(i).name filesep 'values.mat'])   
    try
        VeroXB_all(i,1) = VeroXB;
        VeroMeas_all(i,1) = VeroMeas;
        retreatXB_all(i,1) = retreatXB;
        retreatMeas_all(i,1) = retreatMeas;
    catch
        VeroXB_all(i,1) = NaN;
        VeroMeas_all(i,1) = NaN;
        retreatXB_all(i,1) = NaN;
        retreatMeas_all(i,1) = NaN;
    end
    clear VeroXB VeroMeas retreatXB retreatMeas 
end

save([dirs.output filesep 'skill.mat'], 'VeroXB_all','VeroMeas_all', 'retreatXB_all','retreatMeas_all');
dlm = {'VeroXB','VeroMeas', 'retreatXB','retreatMeas'};
matrix2latex([VeroXB_all, VeroMeas_all, retreatXB_all, retreatMeas_all], 'filename', 'skill.tex', ...
    'rowjustification','p{1.5cm}',...
    'rowlabel', {somnames_morfo.name}, ...
    'collabel', dlm, ...
    'format',	{'%s', '%4.2f', '%4.2f','%4.2f','%4.2f'});

end