plotting bars with respect to three different y axis

imp_2d_E=[1 0;-1 -3;-2 2;-1 0;5 2;-1 0;3 -1;2 2;1 -2;0 2;2 4;1 2;-3 -2;1 1;-2 -7;-3 -1;1 1];
bar(imp_2d_E)
xtv = get(gca, 'XTick');
xtl = {'ALGO','ALIC','ALRT','ANKR','BRAZ','CAS1','CCJ2','CKIS','MAS1','MAW1','MCM4','NOVM','NYA1','ROTH','SCOR','TIXI','ZAMB'};
xtnew = linspace(1, numel(xtl), numel(xtl));
set(gca, 'XTick',xtnew, 'XTickLabel',xtl, 'XTickLabelRotation',90);
title('Changes','FontWeight','bold','FontSize',22);
xlabel('STATIONS', 'FontWeight', 'bold','FontSize',20);
ylabel('millimeter ', 'FontWeight','bold','FontSize',20)
legend('24h','12h');
set(gca,'FontWeight','bold')
set(gca,'fontsize',20)

The above codes plot the changes in 17 stations with 24-12 h legends. I have two more data sets with 24-12 h legends and I need to plot all three data sets in one figure using the vertical shift between each data set. How can I modify the above codes to do it?