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

% load 
input=input*4;

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

% load 

olivInf_con_input=input(:,1);
oliv10_con_input=input(:,2);
oliv100_con_input=input(:,3);
oliv1000_con_input=input(:,4);
olivInf_con_t_input=input(:,5);
oliv10_con_t_input=input(:,6);
oliv100_con_t_input=input(:,7);
oliv1000_con_t_input=input(:,8);
olivInf_input30x2E4_input=input(:,9);
oliv10_input30x2E4_input=input(:,10);
oliv100_input30x2E4_input=input(:,11);
oliv1000_input30x2E4_input=input(:,12);
olivInf_input9_input=input(:,13);
oliv10_input9_input=input(:,14);
oliv100_input9_input=input(:,15);
oliv1000_input9_input=input(:,16);
control_input=input(:,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_input,'b--','linewidth',2);
h2=plot(t(1:4:end),oliv1000_con_input(1:4:end),'b-x','linewidth',2,'markersize',10);
h3=plot(t,oliv100_con_t_input,'g--','linewidth',2);
h4=plot(t(1:4:end),oliv100_con_input(1:4:end),'g-x','linewidth',2,'markersize',10);
h5=plot(t,oliv10_con_t_input,'r--','linewidth',2);
h6=plot(t(1:4:end),oliv10_con_input(1:4:end),'r-x','linewidth',2,'markersize',10);
h7=plot(t,olivInf_con_t_input,'--','color',[0.5 0.5 0.5],'linewidth',2);
h8=plot(t(1:4:end),olivInf_con_input(1:4:end),'-x','color',[0.5 0.5 0.5],'linewidth',2,'markersize',10);
h9=plot(t(1:4:end),control_input(1:4:end),'k-x','linewidth',2,'markersize',10);

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

h5=plot(t,olivInf_input9_input,'-','color',[0.5 0.5 0.5],'linewidth',1.5);
h5=plot(t(3:4:end),olivInf_input9_input(3:4:end),'o','color',[0.5 0.5 0.5],'linewidth',1.5,'markersize',8);
h6=plot(t,oliv10_input9_input,'-','color','r','linewidth',1.5);
h6=plot(t(1:4:end),oliv10_input9_input(1:4:end),'o','color','r','linewidth',1.5,'markersize',8);
h7=plot(t,oliv100_input9_input,'-','color','g','linewidth',1.5);
h7=plot(t(3:4:end),oliv100_input9_input(3:4:end),'o','color','g','linewidth',1.5,'markersize',8);
h8=plot(t,oliv1000_input9_input,'-','color','b','linewidth',1.5);
h8=plot(t(1:4:end),oliv1000_input9_input(1:4:end),'o','color','b','linewidth',1.5,'markersize',8);
h9=plot(t,control_input,'k','linewidth',2);
set(gca,'Yticklabel',{' '})



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


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






