clear
data=forcingdata;
aero=data.Sul+data.Soot+data.Carbon+data.Biomas;
volc=data.Volcanic;
bar(data.year,volc,1.9,'r');
hold on
plot(data.year,aero);
ylabel('Wymuszanie radiacyjne [W/m^2]');
xlabel('Rok')
axis([1850 2008 -3 0.5])
legend('Aerozol Antropogeniczny.','Wulkany')
gridnew
print -dpng FORCING