clear
clc

load eke_anns_ofes2.mat
datestr(tim);
poslat=lat>=-5&lat<=8;
poslon=lon>=200&lon<=260;
postim=tim>=1993&tim<=2021;
data0=squeeze(nanmean(eke(postim,:,poslat,poslon),4))*1024;
lon=lat(poslat);
lat=dep;
tim=tim(postim);
nt=length(tim);
ny=length(lat);
nx=length(lon);
data1=nan(ny,nx,nt);
run=4;
data2=nan(ny,nx,nt-2*run);
for t=1:nt
    data1(:,:,t)=data0(t,:,:);
end
for j=1:ny 
for i=1:nx
    s=squeeze(data1(j,i,:));
    s1=movmean(s,2*run+1);
    s1(isnan(s))=nan;
    data2(j,i,:)=s1(1+run:end-run);
end
end
data3=trend(data2,1)*10;
data4=mann_kendall(data2,0.05);
data20=data2;
data30=data3;

[x,y]=meshgrid(lon,lat);

% plot
figure(1)
set(gcf,'unit','centimeters','position',[5 5 9 5])
set(gca,'position',[0.06 0.35 0.73 0.6])
cax1=-5;
cax2=5;  
ax1=-5;
ax2=8;
ax3=-150;
ax4=0;
nl=40;
data2t=data3;
data2t=movmean(data2t,1,2,'omitnan');
data2t(isnan(data3))=nan;
data2t(data2t>=cax2)=cax2;
data2t(data2t<=cax1)=cax1;

contourf(x,y,data2t,nl,'linestyle','none');
hold on
% m_contour(x,y,data2t,[0. 0.02],'-','color',[0.1 0.1 0.1]*6,'linewidth',0.3);
plot([0 0],[ax3 ax4],'k:','linewidth',.5)
stipple(x,y,data4,'markersize',1,'density',100,'color',[0.1 0.1 0.1]*5)

load u_anns_ofes2
datestr(tim);
poslat=lat>=-5&lat<=8;
poslon=lon>=200&lon<=260;
postim=tim>=1993&tim<=2021;
bb=2.28e-11;
data0=squeeze(nanmean(u(postim,:,poslat,poslon),4))*1e2;
lon=lat(poslat);
lat=dep;
tim=tim(postim);
nt=length(tim);
ny=length(lat);
nx=length(lon);
data1=nan(ny,nx,nt);
run=4;
data2=nan(ny,nx,nt-2*run);
for t=1:nt
    data1(:,:,t)=data0(t,:,:);
end
for j=1:ny 
for i=1:nx
    s=squeeze(data1(j,i,:));
    s1=movmean(s,2*run+1);
    s1(isnan(s))=nan;
    data2(j,i,:)=s1(1+run:end-run);
end
end
data3=trend(data2,1)*10;
sstm=data3;
sstm=movmean(sstm,1,2,'omitnan');
sstm(isnan(data3))=nan;
sstm=movmean(sstm,1,1,'omitnan');
sstm(isnan(data3))=nan;
level=1:2:50;
[c,h]=contour(x,y,sstm,level,'color',[.1 .1 .1]*5,'linewidth',.5);
clabel(c,h,'labelspacing',720,'color',[0.1 0.1 0.1]*2,'fontsize',7)
[c,h]=contour(x,y,sstm,-level,'--','color',[.1 .1 .1]*5,'linewidth',1);
clabel(c,h,'labelspacing',720,'color',[0.1 0.1 0.1]*2,'fontsize',7)
% [c,h]=contour(x,y,sstm,-[0 0],':','color',[.1 .1 .1]*0,'linewidth',1);
% clabel(c,h,'labelspacing',720,'color',[0.1 0.1 0.1]*0,'fontsize',7)

hold off

h1=colorbar("eastoutside");
set(h1,'position',[0.82 0.46 0.02 0.35])
ticks=[cax1:cax2:cax2];
set(h1,'ticks',ticks,'fontsize',7,'TickLength',0.03,'tickdirection','out')
caxis([cax1 cax2])
% h2=get(h1,'title');
% set(h2,'string',['J m-3',newline,'dec'],'rotation',0)
hold on
hm=text(10.3,-78,['J m-3 / dec'],...
    'rotation',-90,'fontsize',7,'HorizontalAlignment','center');
% hm=m_text(178,9,['Satellite'],...
%     'rotation',0,'fontsize',7,'HorizontalAlignment','center');
% hm=m_text(220.,-1,['TAO'],...
%     'rotation',0,'fontsize',7,'HorizontalAlignment','center');
hold off
load bwr_jet
colormap(bwr_jet)
% colormap(othercolor('BuOr_12'))
hc=cbarrow;
uistack(hc,'top')

xtick=-4:4:10;
ytick=-150:50:0;

yticklabel={'150','100','50','0'};
xticklabel={'4°S','0°','4°N','8°N'};

set(gca,'tickdir','out')
set(gca,'xtick',xtick)
set(gca,'xticklabel',xticklabel)
set(gca,'ytick',ytick)
set(gca,'yticklabel',yticklabel)
% set(gca,'yminortick','on')
set(gca,'ticklength',[0.01 0.025])

set(gca,'fontsize',7)

axis([ax1 ax2 ax3 ax4])