
%
CNAM='NBF20TROC_f19_tn11_01';
YEAR1=1850;
YEARN=2005;

fff=15;
%
load(['../Data/timeseries_' CNAM '_' sprintf('%4.4i-%4.4i',YEAR1,YEARN) '.mat']);
icearea_NSIDC;

%
ny=YEARN-YEAR1+1;

%
VARS={'N_mar','N_sep','S_mar','S_sep'};
VARS1={'N_mar (NSIDC)','N_sep (NSIDC)','S_mar (NSIDC)','S_sep (NSIDC)'};

yaxisss=[0 20];

figure;
h(1)=plot([1:ny]+YEAR1-1,ts_yr.aice_N_mar/1e6,'color','k','linewidth',1.5);hold on;
h(2)=plot([1:ny]+YEAR1-1,ts_yr.aice_N_sep/1e6,'color',[1 1 1]*0.7,'linewidth',1.5);hold on
h(3)=plot([1:ny]+YEAR1-1,ts_yr.aice_S_mar/1e6,'--','color','k','linewidth',1.5);hold on;
h(4)=plot([1:ny]+YEAR1-1,ts_yr.aice_S_sep/1e6,'--','color',[1 1 1]*0.7,'linewidth',1.5);hold on;

h0(1)=plot(1979:2005,ice_N_mar(1:27),'color',[204 0 0]/255,'linewidth',1.5);hold on;
h0(2)=plot(1979:2005,ice_N_sep(1:27),'color',[255 153 153]/255,'linewidth',1.5);hold on;
h0(3)=plot(1979:2005,ice_S_mar(1:27),'--','color',[204 0 0]/255,'linewidth',1.5);hold on;
h0(4)=plot(1979:2005,ice_S_sep(1:27),'--','color',[255 153 153]/255,'linewidth',1.5);hold on;

set(gca,'fontsize',fff);
xlabel('years','fontsize',fff);
ylabel('sea ice area (10^6 km^2)','fontsize',fff);

set(gca,'xlim',[0 ny]+YEAR1-1,'fontweight','bold');
set(gca,'xtick',1850:25:2000,'Xminortick','on','Yminortick','on');
legend(h,VARS,'fontweight','bold','interpreter','none','Location','west');
ah=axes('position',get(gca,'position'),'visible','off');
legend(ah,h0,VARS1,'fontweight','bold','interpreter','none','Location','east');
ylim(yaxisss);

set(gca,'fontsize',fff);
print(gcf,'-dpng','-r600',['plot_timeseries_ice_' CNAM '_' sprintf('%4.4i-%4.4i',YEAR1,YEARN) '.png'])
print(gcf,'-depsc',['plot_timeseries_ice_' CNAM '_' sprintf('%4.4i-%4.4i',YEAR1,YEARN) '.eps'])
