clear all
close all

pathi='/data1/home/yhtseng/samuel0302/OMIP2021/S01_timeseries/result/'
patho = '/data1/home/yhtseng/samuel0302/OMIP2021/S01_timeseries/figure/';
addpath '/data1/home/yhtseng/samuel0302/cbrewer/';
addpath('/data1/home/yhtseng/samuel0302/plot_tools/')

%[ha, hfigure]= tight_subplot_cm(1, 1, [.4 .1], [1.3 1], [2.5 .7], 6, 20);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
load ([pathi 'cal03AMOC_v2.2.mat']);
%load ([pathi 'cal03AMOC.mat']);% variables: 'AMOC_iaf','iafperiod','AMOC_jra','jraperiod'
load (['~/cal03OOBSAMOC.mat']);% variables: 'OBS_MOC_mon','OBS_MOC_year' 

%=======================
% annual mean
%=======================
for iyr = 1:length(AMOC_iaf)/12
  annual_AMOC_iaf(iyr) = mean(AMOC_iaf((iyr-1)*12+1:(iyr-1)*12+12));
end
for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra1(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end
for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra2(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end
for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra3(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end
for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra4(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end
for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra5(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end
for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra5(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end
for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra6(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end
for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra7(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end
for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra8(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end
for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra9(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end

for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra10(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end

for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra_CTL(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end
for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra_oldGM(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end

for iyr = 1:length(AMOC_jra)/12
  annual_AMOC_jra_oldGM_3(iyr) = mean(AMOC_jra((iyr-1)*12+1:(iyr-1)*12+12));
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%axes(ha(1));
%plot([1948:1:2009],annual_AMOC_iaf,'color','r','linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
%hold on;plot([1978:1:2018],annual_AMOC_jra_oldGM(1:41),'color','k','linestyle','o','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
%hold on;plot([2013:1:2018],annual_AMOC_jra_oldGM_3(1:6),'color','b','linestyle','--','marker','o','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
hold on;plot([1997:1:2018],annual_AMOC_CTL(1:22),'color','y','linestyle','--','marker','o','MarkerSize',7,'MarkerFaceColor',[228/255,200/255,200/255],'Linewidth',1);
hold on;plot([1997:1:2018],annual_AMOC_jra1(1:22),'color','b','linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
hold on;plot([1997:1:2018],annual_AMOC_jra2(1:22),'color','r','linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
hold on;plot([1997:1:2018],annual_AMOC_jra3(1:22),'color','[0.3 0.3 0.7]','linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
hold on;plot([1997:1:2018],annual_AMOC_jra4(1:22),'color','[0.9 0.5 0.9]','linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
%hold on;plot([1997:1:2018],annual_AMOC_jra5(1:22),'color','[0.7 0.1 0.5]','linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);%hold on;plot([1997:1:2018],annual_AMOC_jra6(1:22),'color','[0.1 0.7 0.9]','linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
hold on;plot([1997:1:2018],annual_AMOC_jra6(1:22),'color','[0.9 0.5 0.3]','linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
%hold on;plot([1997:1:2018],annual_AMOC_jra8(1:22),'color','[0.4 0.9 0.6]','linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
%hold on;plot([1997:1:2018],annual_AMOC_jra9(1:22),'color','[0.8 0.8 0.6]','linestyle','--','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
%hold on;plot([1997:1:2018],annual_AMOC_jra10(1:22),'color','[0.8 0.5 0.3]','linestyle','--','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
%hold on;plot([1997:1:2018],annual_AMOC_jra11(1:22),'color','[0.9 0.9 0.6]','linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);

hold on;plot([2004:2018],OBS_MOC_year,'color',[0/255,0/255,0/255],'linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[0/255,0/255,0/255],'Linewidth',2);
set(gca,'fontsize',12);
ylabel({'AMOC [SV]'},'Fontsize',14,'color','k')
legend({'v2.8.1o 1st ','v2.8.1o1 1st','v2.8.1o 1st ','v2.8.1o1 1st','RAPID observation'},'Fontsize',12)

legend({'v2.8.1j 6th ','v2.8.1k 6th','v2.8.1l 5th','v2.8.1o 4th','RAPID observation'},'Fontsize',12)
%legend({'v2.8.1i 1st','v2.8.1j 1st','v2.8.1k 1st','v2.8.1l 1st','v2.8.1m 1st','v2.8.1e 1st','RAPID observation'},'Fontsize',12)

legend({'v2.8.1f 4th','v2.8.1g 4th','v2.8.1h 4th','v2.8.1i 4th','v2.8.1j 4th','v2.8.1k 2nd','v2.8.1l 2nd','RAPID observation'},'Fontsize',12)

legend boxoff
title({'AMOC at 26.5^oN '},'Fontsize',14)
grid on;
set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'TickDir','out')
ax=gca
set(ax, {'XColor', 'YColor'}, {'k','k'});
set(gcf,'color','w'); 
%set(gca,'color','none'); % Sets axes background
%export_fig([patho 'Figure03_AMOC26'], '-pdf' ,'-png', '-transparent','-painters');
%export_fig([patho 'Figure03_AMOC26_old'], '-pdf' ,'-png', '-transparent','-painters');
%axis([2012 2020 14.5 18.5])










hold on;plot([1997:1:2018],annual_AMOC_jra1(1:22),'color','b','linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
hold on;plot([1997:1:2018],annual_AMOC_jra4(1:22),'color','[0.9 0.5 0.9]','linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[228/255,0/255,0/255],'Linewidth',2);
hold on;plot([2004:2018],OBS_MOC_year,'color',[0/255,0/255,0/255],'linestyle','-','marker','none','MarkerSize',35,'MarkerFaceColor',[0/255,0/255,0/255],'Linewidth',2);
set(gca,'fontsize',12);
ylabel({'AMOC [SV]'},'Fontsize',14,'color','k')
legend({'hi vdc (g)','lo vdc (j)','RAPID observation'},'Fontsize',12)

legend boxoff
title({'AMOC at 26.5^oN '},'Fontsize',14)
grid on;
set(gca,'XMinorTick','on','YMinorTick','on')
set(gca,'TickDir','out')
ax=gca
set(ax, {'XColor', 'YColor'}, {'k','k'});
set(gcf,'color','w');