%Better do it on laptop because the dashed curve will look better on that display
% This goes to generate the time series 
% First I need to copy the time_series.exl into matlab

% load 

olivInf_con_co2=co2(:,1);
oliv10_con_co2=co2(:,2);
oliv100_con_co2=co2(:,3);
oliv1000_con_co2=co2(:,4);
olivInf_con_t_co2=co2(:,5);
oliv10_con_t_co2=co2(:,6);
oliv100_con_t_co2=co2(:,7);
oliv1000_con_t_co2=co2(:,8);
olivInf_co230x2E4_co2=co2(:,9);
oliv10_co230x2E4_co2=co2(:,10);
oliv100_co230x2E4_co2=co2(:,11);
oliv1000_co230x2E4_co2=co2(:,12);
olivInf_co29_co2=co2(:,13);
oliv10_co29_co2=co2(:,14);
oliv100_co29_co2=co2(:,15);
oliv1000_co29_co2=co2(:,16);
control_co2=co2(:,17);

% colors: 
% black
% '0.5 0.5 0.5','c','m','r','g','b' and its dashed
% [0.6 0.65 0.1],[0.5 0.75 0.5][0.7 0.2 0.2][0.5 0.75 1]

% figure 1

t=2020:2099

figure(1)
hold on
h1=plot(t,oliv1000_con_t_co2,'b--','linewidth',2);
h2=plot(t,oliv1000_con_co2,'b','linewidth',5);
h3=plot(t,oliv100_con_t_co2,'g--','linewidth',2);
h4=plot(t,oliv100_con_co2,'g','linewidth',5);
h5=plot(t,oliv10_con_t_co2,'r--','linewidth',2);
h6=plot(t,oliv10_con_co2,'r','linewidth',5);
h7=plot(t,olivInf_con_t_co2,'--','color',[0.5 0.5 0.5],'linewidth',2);
h8=plot(t,olivInf_con_co2,'color',[0.5 0.5 0.5],'linewidth',5);
h9=plot(t,control_co2,'k','linewidth',2);
set(gca,'Yticklabel',{' '})

set(gca,'fontsize',24)
grid on
set(figure(1),'position',[100,100,450,600]);
xlim=[2020 2100];

set(gcf,'PaperPositionMode','auto')
set(gcf,'color','w')


print('group1_co2','-depsc')

clf

hold on
h1=plot(t,olivInf_co230x2E4_co2,'-','color',[0.5 0.5 0.5],'linewidth',1.5);
h1=plot(t(1:4:end),olivInf_co230x2E4_co2(1:4:end),'^','color',[0.5 0.5 0.5],'linewidth',1.5,'markersize',8);
h2=plot(t,oliv10_co230x2E4_co2,'-','color','r','linewidth',1.5);
h2=plot(t(3:4:end),oliv10_co230x2E4_co2(3:4:end),'^','color','r','linewidth',1.5,'markersize',8);
h3=plot(t,oliv100_co230x2E4_co2,'-','color','g','linewidth',1.5);
h3=plot(t(1:4:end),oliv100_co230x2E4_co2(1:4:end),'^','color','g','linewidth',1.5,'markersize',8);
h4=plot(t,oliv1000_co230x2E4_co2,'-','color','b','linewidth',1.5);
h4=plot(t(3:4:end),oliv1000_co230x2E4_co2(3:4:end),'^','color','b','linewidth',1.5,'markersize',8);

h5=plot(t,olivInf_co29_co2,'-','color',[0.5 0.5 0.5],'linewidth',1.5);
h5=plot(t(3:4:end),olivInf_co29_co2(3:4:end),'o','color',[0.5 0.5 0.5],'linewidth',1.5,'markersize',8);
h6=plot(t,oliv10_co29_co2,'-','color','r','linewidth',1.5);
h6=plot(t(1:4:end),oliv10_co29_co2(1:4:end),'o','color','r','linewidth',1.5,'markersize',8);
h7=plot(t,oliv100_co29_co2,'-','color','g','linewidth',1.5);
h7=plot(t(3:4:end),oliv100_co29_co2(3:4:end),'o','color','g','linewidth',1.5,'markersize',8);
h8=plot(t,oliv1000_co29_co2,'-','color','b','linewidth',1.5);
h8=plot(t(1:4:end),oliv1000_co29_co2(1:4:end),'o','color','b','linewidth',1.5,'markersize',8);
h9=plot(t,control_co2,'k','linewidth',2);
set(gca,'Yticklabel',{' '})




set(gca,'fontsize',24)
grid on
set(figure(1),'position',[100,100,450,600]);
xlim=[2020 2100];

set(gcf,'PaperPositionMode','auto')
set(gcf,'color','w')
print('group2_co2','-depsc')
clf


