Esperiment with this:
T1 = readtable('dummy_model_data.txt', 'HeaderLines',1);
T1.Properties.VariableNames = {'Model','Min','Max','Mean'};
BlueBars = [reshape([reshape(T1.Min, 7, []); nan(1,3)], 1,[]); reshape([reshape(T1.Mean, 7, []); nan(1,3)], 1,[])];
RedBars = [reshape([reshape(T1.Max, 7, []); nan(1,3)], 1,[]); reshape([reshape(T1.Mean, 7, []); nan(1,3)],1,[])];
Models = reshape([reshape(T1.Model, 7, []); cell(1,3)], [], 1);
figure
plot([1;1]*(1:size(BlueBars,2)), BlueBars, 'b', 'LineWidth',15)
hold on
plot([1;1]*(1:size(RedBars,2)), RedBars, 'r', 'LineWidth',15)
hold off
set(gca, 'XTick',(1:size(BlueBars,2)), 'XTickLabel',Models, 'XTickLabelRotation',90)
pos = get(gcf,'Position');
text(4:8:24, ones(1,3)*-80, compose('%s 2007', ['JAN';'MAR';'MAY']))
set(gcf, 'Position', pos+[0 -250 500 250])
producing:
.
