clc
clear
close all;

load('I:\OHC_trend_by_IAP\LFK\SLAT\LFK_ctrl_SLAT_F50.mat')
load('I:\OHC_trend_by_IAP\LFK\SLAT\LFK_NAEU_SLAT_F50.mat')
S_NAEU_F=S_NAEU-S_ctrl;

load('I:\OHC_trend_by_IAP\LFK\SLAT\LFK_EA_SLAT_F50.mat')
S_EA_F=S_EA-S_ctrl;
S_1920_1980=S_NAEU_F+S_EA_F./1.2;
S_1980_2005=S_EA_F./1.2-S_NAEU_F;
%鎻掑??
x=0:1:360;%鎻掑?间箣鍚庣殑lon,lat
y=-90:1:90;
[X,Y]=meshgrid(y,x);
SS_1980_2005=griddata(latData,lonData,S_1980_2005,X,Y);
SS_1920_1980=griddata(latData,lonData,S_1920_1980,X,Y);


figure('Units','centimeters','Position',[8,2,21,15]);
subplot('position',[0.08 0.55 0.7 0.4])
% [long,latt]=meshgrid(double(LON),double(LAT));
contourf(SS_1920_1980',200,'linestyle','none');
hold on
grid on
load coast;
lat=lat+91;
long=long+(long<0)*360+1;
long(long>359)=NaN;
plot(long,lat,'color','k')
ylim([-10+90 30+90])
xlim([90 180])
% colorbar;
set(gca,'xtick',90:20:180)
set(gca,'Xticklabel',{'90E','110E','130E','150E','180',...
    '150W','120W','90W'})
set(gca,'Ytick',80:10:60+90,'Yticklabel',{'10S','EQ',...
    '10N','20N','30N'},'fontname','Times New Roman','fontsize',12,'fontweight','normal')
caxis([-0.61 0.61]);
title('(c)S0-55m NAEU-F + ASIA-F ','fontsize',15)

subplot('position',[0.08 0.05 0.7 0.4])
% [long,latt]=meshgrid(double(LON),double(LAT));
contourf(SS_1980_2005',200,'linestyle','none');
hold on
grid on
load coast;
lat=lat+91;
long=long+(long<0)*360+1;
long(long>359)=NaN;
plot(long,lat,'color','k')
ylim([-10+90 30+90])
xlim([90 180])
% colorbar;
set(gca,'xtick',90:20:180)
set(gca,'Xticklabel',{'90E','110E','130E','150E','180',...
    '150W','120W','90W'})
set(gca,'Ytick',80:10:60+90,'Yticklabel',{'10S','EQ',...
    '10N','20N','30N'},'fontname','Times New Roman','fontsize',12,'fontweight','normal')
caxis([-0.61 0.61]);
title('(d)S0-55m ASIA-F - NAEU-F','fontsize',15)

load('J:\盐度\年代际变化\color\b_r3.mat')
L=length(colors);
L2=floor(L/2);
a=colors(1:L2,:);b=colors(L2:end,:);
c=[a;b];
colormap(c);
hcb=colorbar('position',[0.8,0.08,0.025,0.85]);
% print(gcf,'-dpng','-r600','LFK_Ctrl_SLAT_F50');
% print -depsc2 -painters LFK_1980qh.eps
% print(gcf,'-dpng','-r600','LFK_1980qh');