clear

%grid
lats=-89:89; nz=90;

load FWHT_year0.mat %from U^bar*T^bar(at T point)
aFW_yr0=aFW; gFW_yr0=gFW; pFW_yr0=pFW;
aHT_yr0=aHT; gHT_yr0=gHT; pHT_yr0=pHT;
kk=find(lats<-35|lats>70); aFW_yr0(kk,:)=nan; aHT_yr0(kk,:)=nan;
kk=find(lats<-35|lats>65); pFW_yr0(kk,:)=nan; pHT_yr0(kk,:)=nan;

load FWHT_year.mat %from U^bar*T^bar(at U point)
aFW_yr=aFW; gFW_yr=gFW; pFW_yr=pFW;
aHT_yr=aHT; gHT_yr=gHT; pHT_yr=pHT;
kk=find(lats<-35|lats>70); aFW_yr(kk,:)=nan; aHT_yr(kk,:)=nan;
kk=find(lats<-35|lats>65); pFW_yr(kk,:)=nan; pHT_yr(kk,:)=nan;

load FWHT_year2.mat %from (U*T)^bar(at T point)
aFW_yr2=aFW; gFW_yr2=gFW; pFW_yr2=pFW;
aHT_yr2=aHT; gHT_yr2=gHT; pHT_yr2=pHT;
kk=find(lats<-35|lats>70); aFW_yr2(kk,:)=nan; aHT_yr2(kk,:)=nan;
kk=find(lats<-35|lats>65); pFW_yr2(kk,:)=nan; pHT_yr2(kk,:)=nan;


figure(1)
subplot(211)
plot(lats,[gFW_yr0 aFW_yr0 pFW_yr0],'linewidth',2)
grid on
xlabel('Latitude')
xlim([-90 90])
legend('Global','Atlantic','Pacific+Indian');
title('Meridional seawater transport (in Sv)');

subplot(212)
plot(lats,[gHT_yr0 aHT_yr0 pHT_yr0],'linewidth',2)
grid on
xlabel('Latitude')
xlim([-90 90])
legend('Global','Atlantic','Pacific+Indian');
title('Meridional Heat Transport (in PW)');

figure(2)
subplot(211)
plot(lats,[gFW_yr0-gFW_yr aFW_yr0-aFW_yr pFW_yr0-pFW_yr],'linewidth',2)
grid on
xlabel('Latitude')
xlim([-90 90])
legend('Global','Atlantic','Pacific+Indian');
title('Meridional seawater transport (in Sv)');

subplot(212)
plot(lats,[gHT_yr0-gHT_yr aHT_yr0-aHT_yr pHT_yr0-pHT_yr],'linewidth',2)
grid on
xlabel('Latitude')
xlim([-90 90])
legend('Global','Atlantic','Pacific+Indian');
title('Meridional Heat Transport (in PW)');

figure(3)
subplot(211)
plot(lats,[gFW_yr2-gFW_yr aFW_yr2-aFW_yr pFW_yr2-pFW_yr],'linewidth',2)
grid on
xlabel('Latitude')
xlim([-90 90])
legend('Global','Atlantic','Pacific+Indian');
title('Meridional seawater transport (in Sv)');

subplot(212)
plot(lats,[gHT_yr2-gHT_yr aHT_yr2-aHT_yr pHT_yr2-pHT_yr],'linewidth',2)
grid on
xlabel('Latitude')
xlim([-90 90])
legend('Global','Atlantic','Pacific+Indian');
title('Meridional Heat Transport (in PW)');
% 1% difference between gHT_yr0 and gHT_y
% no bother go mk_statUTVT_new.m

load FWHT_month.mat
kk=find(lats<-35|lats>70); aFW(kk,:)=nan; aHT(kk,:)=nan;
kk=find(lats<-35|lats>65); pFW(kk,:)=nan; pHT(kk,:)=nan;
aFW_yrr=mean(aFW,2); gFW_yrr=mean(gFW,2); pFW_yrr=mean(pFW,2);
aHT_yrr=mean(aHT,2); gHT_yrr=mean(gHT,2); pHT_yrr=mean(pHT,2);
figure(4)
subplot(211)
plot(lats,[gFW_yrr-gFW_yr aFW_yrr-aFW_yr pFW_yrr-pFW_yr],'linewidth',2)
grid on
xlabel('Latitude')
xlim([-90 90])
legend('Global','Atlantic','Pacific+Indian');
title('Meridional seawater transport (in Sv)');

subplot(212)
plot(lats,[gHT_yrr-gHT_yr aHT_yrr-aHT_yr pHT_yrr-pHT_yr],'linewidth',2)
grid on
xlabel('Latitude')
xlim([-90 90])
legend('Global','Atlantic','Pacific+Indian');
title('Meridional Heat Transport (in PW)');

% 10% difference between gHT_yr and gHT_yrr
% need (U*T)^bar as opposed to U^bar*T^bar