clear

load SECTION_year
bt_strmfn_yr=bt_strmfn;

figure(1)
fig_ps=[752 56 1537 1038];
set(gcf,'Position',fig_ps)

plotllc(bt_strmfn)
caxis([-60 160])
colormap(jet)
title('Mean barotropic streamfunction')

load SECTION_month
bt_strmfn_mn=mean(bt_strmfn,3);
bt_strmfn_sd=std(bt_strmfn,0,3);

figure(2)
set(gcf,'Position',fig_ps)
plotllc(bt_strmfn_mn-bt_strmfn_yr)

figure(3)
set(gcf,'Position',fig_ps)

subplot(211)
plotllc(bt_strmfn_mn)
caxis([-60 160])
colormap(jet)
title('Mean barotropic streamfunction')

subplot(212)
plotllc(bt_strmfn_sd)
%caxis([-60 160])
colormap(jet)
title('Std barotropic streamfunction')


figure(4)
set(gcf,'Position',fig_ps)
TT=3:14;
ii=[1 8 9 3 2 14 13 18]; @gcmfaces_lines_pairs.m
tt={'Bering Strait', 'Davis Strait', 'Denmark Strait', ...
    'Florida Strait', 'Gibraltar', 'Indonesia Throughflow', ...
    'Drake Passage', 'Australia Antarctica'};

for i=1:8
subplot(4,2,i)
a1=squeeze(fldTRANSPORTS(ii(i),:,:));
if i==6 %ITF
a1=-squeeze(nansum(fldTRANSPORTS(14:17,:,:),1));
end
a2=nansum(a1);
plot(TT,a2,'linew',2)
grid
xlim([2 15])
%xlabel('month from 2020/01')
a3=sprintf('%s (%7.2f Sv)' ,tt{i},mean(a2));
title(a3)
%title([tt{i} ' (' num2str(mean(a2)) 'Sv)'])
end
subplot(427)
xlabel('month from 2020/01')
subplot(428)
xlabel('month from 2020/01')

figure(3)
print -dpng lk_trsp_baro

figure(4)
print -dpng lk_trsp_sect