% fluxtoa
load('/work/cgu025/noresm/micom_diag/micom/toa/Data/toaflux_NBF1850OC_f19_tn11_02_1001-2000.mat');

% fluxo
load('/work/cgu025/noresm/micom_diag/micom/toa/Data/toasrf_NBF1850OC_f19_tn11_02_1001-2000.mat');
fluxo=fluxo/5.1010e+14;


% ts_yr.sst_yr & ts_yr.sss_yr
load('/work/cgu025/noresm/micom_diag/micom/timeseries/Data/timeseries_NBF1850OC_f19_tn11_02_1001-2000.mat');

% amoc265
amoc265=ncread('/work/cgu025/noresm/micom_diag/NBF1850OC_f19_tn11_01/moc/Figs/NBF1850OC_f19_tn11_02.hy_maxmoc_1001-2000.nc','amoc265');




yr=1001:2000;

nn=5
figure;
%orient tall;

subplot(5,1,1);
plot(yr,fluxtoa,'color',0.7*ones(1,3));hold on;
plot(yr,rm(fluxtoa,10),'color','k');
set(gca,'xticklabel',[],'xminortick','on','yminortick','on');
ylabel('Net TOA (W m^{-2})')

%subplot(5,1,2);
%plot(yr,fluxo,'color',0.7*ones(1,3));hold on;
%plot(yr,rm(fluxo,10),'color','k');
%set(gca,'xticklabel',[],'xminortick','on','yminortick','on');
%ylabel('Ocn srf h f (W m^{-2})')

subplot(5,1,2);
plot(yr,ts_yr.t2m_yr,'color',0.7*ones(1,3));hold on;
plot(yr,rm(ts_yr.t2m_yr,10),'color','k');
set(gca,'xticklabel',[],'xminortick','on','yminortick','on');
ylabel('T_{2m} (^{\circ}C)')

subplot(5,1,3);
plot(yr,ts_yr.sst_yr,'color',0.7*ones(1,3));hold on;
plot(yr,rm(ts_yr.sst_yr,10),'color','k');
set(gca,'xticklabel',[],'xminortick','on','yminortick','on');
ylabel('SST (^{\circ}C)')

subplot(5,1,4);
plot(yr,ts_yr.sss_yr,'color',0.7*ones(1,3));hold on;
plot(yr,rm(ts_yr.sss_yr,10),'color','k');
ylim([33.9 34.1]);
set(gca,'xticklabel',[],'xminortick','on','yminortick','on');
ylabel('SSS (g kg^{-1})');

subplot(5,1,5);
plot(yr,amoc265,'color',0.7*ones(1,3));hold on;
plot(yr,rm(amoc265,10),'color','k');
set(gca,'xminortick','on','yminortick','on');
ylabel('AMOC 26.5^{\circ}N (Sv)');

print('-depsc','plot_timeseries.eps');
print('-dpng','-r300','plot_timeseries.png');
