load Outputv2_2XCO2ramp25yr_nonoise;
t1slow = Ts1; t2slow = Ts2; tslow = Ts;
Ftslow = Ft;
CO2slow = CO2;

load Outputv2_2XCO2ramp25yr_nonoise_nowavafdk;
t1slow_nowa = Ts1; t2slow_nowa = Ts2; tslow_nowa = Ts;
Ftslow_nowa = Ft;
time_nowa = Time/YEAR;

load Outputv2_2XCO2step_nonoise;
t1fast = Ts1; t2fast = Ts2; tfast = Ts;
Ftfast = Ft;
time = Time/YEAR;
CO2fast = CO2;

load Outputv2_2XCO2step_nonoise_nowavafdk;
t1fast_nowa = Ts1; t2fast_nowa = Ts2; tfast_nowa = Ts;
Ftfast_nowa = Ft;
time_nowa = Time/YEAR;

load Output_CTRL_noise1;
t1n1 = Ts1; t2n1 = Ts2; tn1 = Ts;
Ftn1 = Ft;
timen = Time/YEAR;

figure(1), clf
[h]=plot(time,CO2fast,'b',time,CO2slow,'r');
set(h,'Linewidth',2);
set(gca,'Fontsize',15);
xlabel('Time (year)','Fontsize',15);
ylabel('CO2 concentration (ppm)','Fontsize',15);

figure(2), clf
[h]=plot(time,tfast-tfast(1),'b',time,tslow-tslow(1),'r',time_nowa,tfast_nowa-tfast_nowa(1),'b--',time_nowa,tslow_nowa-tslow_nowa(1),'r--');
set(h,'Linewidth',2);
set(gca,'Fontsize',15);
xlabel('Time (year)','Fontsize',15);
ylabel('Surface warming (K)','Fontsize',15);

figure(3), clf
[h]=plot(time,tfast-tfast(1),'b',time,tslow-tslow(1),'r',time_nowa,tfast_nowa-tfast_nowa(1),'b--',time_nowa,tslow_nowa-tslow_nowa(1),'r--');
set(h,'Linewidth',2);
set(gca,'Fontsize',15);
xlabel('Time (year)','Fontsize',15);
ylabel('Surface warming (K)','Fontsize',15);

figure(10), clf
subplot(3,1,1), plot(time,tfast,'b',time,tslow,'r',timen,tn1,'g',time_nowa,tfast_nowa,'b--',time_nowa,tslow_nowa,'r--');
subplot(3,1,2), plot(time,t1fast,'b',time,t1slow,'r',timen,t1n1,'g',time_nowa,t1fast_nowa,'b--',time_nowa,t1slow_nowa,'r--');
subplot(3,1,3), plot(time,t2fast,'b',time,t2slow,'r',timen,t2n1,'g',time_nowa,t2fast_nowa,'b--',time_nowa,t2slow_nowa,'r--');

figure(11), clf
plot(time(1:end-1),Ftfast,'b',time(1:end-1),Ftslow,'r',timen(1:end-1),Ftn1,'g',time_nowa(1:end-1),Ftfast_nowa,'b--',time_nowa(1:end-1),Ftslow_nowa,'r--');