; Overplot
plot,years,t1,ystyle=8,/xs
axis,/yax,/save,yr=makerange(co21)
oplot,years,co21,col=3

runs=[101,201,301,401,501]
cols=[1,2,6,7,8]
th=[1,1,2,2,3]

axis,ystyle=4,yr=[-1,1],/save,/yax

for i=0,n_elements(runs)-1 do begin

  c=two(years,t1,co21,run=runs[i])

  oplot,years,c,th=th[i],col=cols(i)

  oplot,[0,inc*runs[i]]+max(years)*7/10.,-0.6-i*0.05*[1,1],th=th[i],col=cols[i]

endfor

gettwogifs,out='three'

end