clear
clc

%% slp
f='/Volumes/Backup Plus/backup/polar/ERA5/surface7922.nc';
for m=27:38
    mclim(:,:,m-26)=mean(squeeze(ncread(f,'msl',[1 1 1 m],[1440 361 1 30],[1 1 1 12])/100),3);
end
x0=ncread(f,'longitude');
y0=ncread(f,'latitude');
[X,Y]=meshgrid(x0,y0);
X=double(X);Y=double(Y);
ma(:,:,1)=permute(squeeze(ncread(f,'msl',[1 1 1 519],[1440 361 1 1],[1 1 1 1])/100)-mclim(:,:,1),[2 1]);
ma(:,:,2)=permute(squeeze(ncread(f,'msl',[1 1 1 520],[1440 361 1 1],[1 1 1 1])/100)-mclim(:,:,2),[2 1]);
ma(:,:,3)=permute(squeeze(ncread(f,'msl',[1 1 1 521],[1440 361 1 1],[1 1 1 1])/100)-mclim(:,:,3),[2 1]);
ma(:,:,4)=permute(squeeze(ncread(f,'msl',[1 1 1 522],[1440 361 1 1],[1 1 1 1])/100)-mclim(:,:,4),[2 1]);
ma(:,:,5)=permute(squeeze(ncread(f,'msl',[1 1 1 523],[1440 361 1 1],[1 1 1 1])/100)-mclim(:,:,5),[2 1]);
ma(:,:,6)=permute(squeeze(ncread(f,'msl',[1 1 1 524],[1440 361 1 1],[1 1 1 1])/100)-mclim(:,:,6),[2 1]);
ma(:,:,7)=permute(squeeze(ncread(f,'msl',[1 1 1 525],[1440 361 1 1],[1 1 1 1])/100)-mclim(:,:,7),[2 1]);
ma(:,:,8)=permute(squeeze(ncread(f,'msl',[1 1 2 526],[1440 361 1 1],[1 1 1 1])/100)-mclim(:,:,8),[2 1]);
ma(:,:,9)=permute(squeeze(ncread(f,'msl',[1 1 2 527],[1440 361 1 1],[1 1 1 1])/100)-mclim(:,:,9),[2 1]);
ma(:,:,10)=permute(squeeze(ncread(f,'msl',[1 1 2 528],[1440 361 1 1],[1 1 1 1])/100)-mclim(:,:,10),[2 1]);
f0='/Volumes/Backup Plus/backup/polar/ERA5/surface23.nc';
ma(:,:,11)=permute(squeeze(ncread(f0,'msl',[1 1 1 1],[1440 361 1 1],[1 1 1 1])/100)-mclim(:,:,11),[2 1]);
ma(:,:,12)=permute(squeeze(ncread(f0,'msl',[1 1 1 2],[1440 361 1 1],[1 1 1 1])/100)-mclim(:,:,12),[2 1]);
ano(:,:,1)=nanmean(ma(:,:,1:3),3);
ano(:,:,2)=nanmean(ma(:,:,4:6),3);
ano(:,:,3)=nanmean(ma(:,:,7:9),3);
ano(:,:,4)=nanmean(ma(:,:,10:12),3);

%% t u v
load lon
load lat
for m=27:38
    uclim(:,:,m-26)=mean(squeeze(ncread(f,'u10',[1 1 1 m],[1440 361 1 30],[1 1 1 12])),3);
    vclim(:,:,m-26)=mean(squeeze(ncread(f,'v10',[1 1 1 m],[1440 361 1 30],[1 1 1 12])),3);
end

u0(:,:,1)=permute(squeeze(ncread(f,'u10',[1 1 1 519],[1440 361 1 1],[1 1 1 1]))-uclim(:,:,1),[2 1]);
u0(:,:,2)=permute(squeeze(ncread(f,'u10',[1 1 1 520],[1440 361 1 1],[1 1 1 1]))-uclim(:,:,2),[2 1]);
u0(:,:,3)=permute(squeeze(ncread(f,'u10',[1 1 1 521],[1440 361 1 1],[1 1 1 1]))-uclim(:,:,3),[2 1]);
u0(:,:,4)=permute(squeeze(ncread(f,'u10',[1 1 1 522],[1440 361 1 1],[1 1 1 1]))-uclim(:,:,4),[2 1]);
u0(:,:,5)=permute(squeeze(ncread(f,'u10',[1 1 1 523],[1440 361 1 1],[1 1 1 1]))-uclim(:,:,5),[2 1]);
u0(:,:,6)=permute(squeeze(ncread(f,'u10',[1 1 1 524],[1440 361 1 1],[1 1 1 1]))-uclim(:,:,6),[2 1]);
u0(:,:,7)=permute(squeeze(ncread(f,'u10',[1 1 1 525],[1440 361 1 1],[1 1 1 1]))-uclim(:,:,7),[2 1]);
u0(:,:,8)=permute(squeeze(ncread(f,'u10',[1 1 2 526],[1440 361 1 1],[1 1 1 1]))-uclim(:,:,8),[2 1]);
u0(:,:,9)=permute(squeeze(ncread(f,'u10',[1 1 2 527],[1440 361 1 1],[1 1 1 1]))-uclim(:,:,9),[2 1]);
u0(:,:,10)=permute(squeeze(ncread(f,'u10',[1 1 2 528],[1440 361 1 1],[1 1 1 1]))-uclim(:,:,10),[2 1]);
u0(:,:,11)=permute(squeeze(ncread(f0,'u10',[1 1 1 1],[1440 361 1 1],[1 1 1 1]))-uclim(:,:,11),[2 1]);
u0(:,:,12)=permute(squeeze(ncread(f0,'u10',[1 1 1 2],[1440 361 1 1],[1 1 1 1]))-uclim(:,:,12),[2 1]);

v0(:,:,1)=permute(squeeze(ncread(f,'v10',[1 1 1 519],[1440 361 1 1],[1 1 1 1]))-vclim(:,:,1),[2 1]);
v0(:,:,2)=permute(squeeze(ncread(f,'v10',[1 1 1 520],[1440 361 1 1],[1 1 1 1]))-vclim(:,:,2),[2 1]);
v0(:,:,3)=permute(squeeze(ncread(f,'v10',[1 1 1 521],[1440 361 1 1],[1 1 1 1]))-vclim(:,:,3),[2 1]);
v0(:,:,4)=permute(squeeze(ncread(f,'v10',[1 1 1 522],[1440 361 1 1],[1 1 1 1]))-vclim(:,:,4),[2 1]);
v0(:,:,5)=permute(squeeze(ncread(f,'v10',[1 1 1 523],[1440 361 1 1],[1 1 1 1]))-vclim(:,:,5),[2 1]);
v0(:,:,6)=permute(squeeze(ncread(f,'v10',[1 1 1 524],[1440 361 1 1],[1 1 1 1]))-vclim(:,:,6),[2 1]);
v0(:,:,7)=permute(squeeze(ncread(f,'v10',[1 1 1 525],[1440 361 1 1],[1 1 1 1]))-vclim(:,:,7),[2 1]);
v0(:,:,8)=permute(squeeze(ncread(f,'v10',[1 1 2 526],[1440 361 1 1],[1 1 1 1]))-vclim(:,:,8),[2 1]);
v0(:,:,9)=permute(squeeze(ncread(f,'v10',[1 1 2 527],[1440 361 1 1],[1 1 1 1]))-vclim(:,:,9),[2 1]);
v0(:,:,10)=permute(squeeze(ncread(f,'v10',[1 1 2 528],[1440 361 1 1],[1 1 1 1]))-vclim(:,:,10),[2 1]);
v0(:,:,11)=permute(squeeze(ncread(f0,'v10',[1 1 1 1],[1440 361 1 1],[1 1 1 1]))-vclim(:,:,11),[2 1]);
v0(:,:,12)=permute(squeeze(ncread(f0,'v10',[1 1 1 2],[1440 361 1 1],[1 1 1 1]))-vclim(:,:,12),[2 1]);

for i=1:12
ua(:,:,i)=griddata(X,Y,u0(:,:,i),lon,lat);
va(:,:,i)=griddata(X,Y,v0(:,:,i),lon,lat);
end

for m=27:38
    tclim(:,:,m-26)=mean(squeeze(ncread(f,'t2m',[1 1 1 m],[1440 361 1 30],[1 1 1 12])),3);
end
ta(:,:,1)=permute(squeeze(ncread(f,'t2m',[1 1 1 519],[1440 361 1 1],[1 1 1 1]))-tclim(:,:,1),[2 1]);
ta(:,:,2)=permute(squeeze(ncread(f,'t2m',[1 1 1 520],[1440 361 1 1],[1 1 1 1]))-tclim(:,:,2),[2 1]);
ta(:,:,3)=permute(squeeze(ncread(f,'t2m',[1 1 1 521],[1440 361 1 1],[1 1 1 1]))-tclim(:,:,3),[2 1]);
ta(:,:,4)=permute(squeeze(ncread(f,'t2m',[1 1 1 522],[1440 361 1 1],[1 1 1 1]))-tclim(:,:,4),[2 1]);
ta(:,:,5)=permute(squeeze(ncread(f,'t2m',[1 1 1 523],[1440 361 1 1],[1 1 1 1]))-tclim(:,:,5),[2 1]);
ta(:,:,6)=permute(squeeze(ncread(f,'t2m',[1 1 1 524],[1440 361 1 1],[1 1 1 1]))-tclim(:,:,6),[2 1]);
ta(:,:,7)=permute(squeeze(ncread(f,'t2m',[1 1 1 525],[1440 361 1 1],[1 1 1 1]))-tclim(:,:,7),[2 1]);
ta(:,:,8)=permute(squeeze(ncread(f,'t2m',[1 1 2 526],[1440 361 1 1],[1 1 1 1]))-tclim(:,:,8),[2 1]);
ta(:,:,9)=permute(squeeze(ncread(f,'t2m',[1 1 2 527],[1440 361 1 1],[1 1 1 1]))-tclim(:,:,9),[2 1]);
ta(:,:,10)=permute(squeeze(ncread(f,'t2m',[1 1 2 528],[1440 361 1 1],[1 1 1 1]))-tclim(:,:,10),[2 1]);
ta(:,:,11)=permute(squeeze(ncread(f0,'t2m',[1 1 1 1],[1440 361 1 1],[1 1 1 1]))-tclim(:,:,11),[2 1]);
ta(:,:,12)=permute(squeeze(ncread(f0,'t2m',[1 1 1 2],[1440 361 1 1],[1 1 1 1]))-tclim(:,:,12),[2 1]);

ano(:,:,5)=nanmean(ta(:,:,1:3),3);
ano(:,:,6)=nanmean(ta(:,:,4:6),3);
ano(:,:,7)=nanmean(ta(:,:,7:9),3);
ano(:,:,8)=nanmean(ta(:,:,10:12),3);

uano(:,:,1)=nanmean(ua(:,:,1:3),3);
uano(:,:,2)=nanmean(ua(:,:,4:6),3);
uano(:,:,3)=nanmean(ua(:,:,7:9),3);
uano(:,:,4)=nanmean(ua(:,:,10:12),3);
vano(:,:,1)=nanmean(va(:,:,1:3),3);
vano(:,:,2)=nanmean(va(:,:,4:6),3);
vano(:,:,3)=nanmean(va(:,:,7:9),3);
vano(:,:,4)=nanmean(va(:,:,10:12),3);

%% ad
f='/Volumes/Backup Plus/backup/polar/ERA5/tuv1000.nc';
PI=3.1415926;
R0=6371004;
for i=27:38
    t=permute(squeeze(ncread(f,'t',[1 1 1 i],[1440 361 1 30],[1 1 1 12])),[2 1 3]);
    u=permute(squeeze(ncread(f,'u',[1 1 1 i],[1440 361 1 30],[1 1 1 12])),[2 1 3]);
    v=permute(squeeze(ncread(f,'v',[1 1 1 i],[1440 361 1 30],[1 1 1 12])),[2 1 3]);
    for n=1:30
        [dtx,dty]=gradient(t(:,:,n));
        dx=gradient(X)*PI/180;
        dy=dx;
        ad(:,:,n)=-((u(:,:,n).*dtx)./(dx.*cos(Y*PI/180))+v(:,:,n).*(dty./dy))/R0;
    end
    adclim(:,:,i-26)=nanmean(ad,3);
end
for i=[519:525 527:528]
    t=permute(squeeze(ncread(f,'t',[1 1 1 i],[1440 361 1 1],[1 1 1 1])),[2 1 3]);
    u=permute(squeeze(ncread(f,'u',[1 1 1 i],[1440 361 1 1],[1 1 1 1])),[2 1 3]);
    v=permute(squeeze(ncread(f,'v',[1 1 1 i],[1440 361 1 1],[1 1 1 1])),[2 1 3]);
    [dtx,dty]=gradient(t);
    dx=gradient(X)*PI/180;
    dy=dx;
    ad2206=-((u.*dtx)./(dx.*cos(Y*PI/180))+v.*(dty./dy))/R0;
    adano(:,:,i-518)=ad2206-adclim(:,:,i-518);
end
t=permute(squeeze(ncread(f,'t',[1 1 2 526],[1440 361 1 1],[1 1 1 1])),[2 1 3]);
u=permute(squeeze(ncread(f,'u',[1 1 2 526],[1440 361 1 1],[1 1 1 1])),[2 1 3]);
v=permute(squeeze(ncread(f,'v',[1 1 2 526],[1440 361 1 1],[1 1 1 1])),[2 1 3]);
[dtx,dty]=gradient(t);
dx=gradient(X)*PI/180;
dy=dx;
ad2206=-((u.*dtx)./(dx.*cos(Y*PI/180))+v.*(dty./dy))/R0;
adano(:,:,8)=ad2206-adclim(:,:,8);
f0='/Volumes/Backup Plus/backup/polar/ERA5/tuv23.nc';
for i=1:2
    t=permute(ncread(f0,'t',[1 1 i],[1440 361 1],[1 1 1]),[2 1]);
    u=permute(ncread(f0,'u',[1 1 i],[1440 361 1],[1 1 1]),[2 1]);
    v=permute(ncread(f0,'v',[1 1 i],[1440 361 1],[1 1 1]),[2 1]);
    [dtx,dty]=gradient(t);
    dx=gradient(X)*PI/180;
    dy=dx;
    ad2206=-((u.*dtx)./(dx.*cos(Y*PI/180))+v.*(dty./dy))/R0;
    adano(:,:,i+10)=ad2206-adclim(:,:,i+10);
end
h=fspecial('average',7);
for i=1:12
    adanof(:,:,i)=imfilter(adano(:,:,i),h,'conv');
end
ano(:,:,9)=nanmean(adanof(:,:,1:3),3);
ano(:,:,10)=nanmean(adanof(:,:,4:6),3);
ano(:,:,11)=nanmean(adanof(:,:,7:9),3);
ano(:,:,12)=nanmean(adanof(:,:,10:12),3);

%%
load ntSIC
for i=27:38
    Cclim(:,:,i-26)=nanmean(C(:,:,i:12:386),3);
end
filepath='/Volumes/Backup Plus/backup/polar/dailySICv2_nrt/';
filelist=dir([filepath,'N*.nc']);
k=length(filelist);
for i=1:167
  filename=[filepath,filelist(i).name];
  SIC=ncread(filename,'F18_ICECON')*100;
  SIC(SIC>100)=nan;
  C0(:,:,i)=SIC;
end

for i=519:528
    Ca(:,:,i-518)=C(:,:,i)-Cclim(:,:,i-518);
end
Ca(:,:,11)=nanmean(C0(:,:,1:31),3)-Cclim(:,:,11);
Ca(:,:,12)=nanmean(C0(:,:,32:59),3)-Cclim(:,:,12);

anoC(:,:,1)=nanmean(Ca(:,:,1:3),3);
anoC(:,:,2)=nanmean(Ca(:,:,4:6),3);
anoC(:,:,3)=nanmean(Ca(:,:,7:9),3);
anoC(:,:,4)=nanmean(Ca(:,:,10:12),3);

edgeC(:,:,1)=nanmean(Cclim(:,:,1:3),3);
edgeC(:,:,2)=nanmean(Cclim(:,:,4:6),3);
edgeC(:,:,3)=nanmean(Cclim(:,:,7:9),3);
edgeC(:,:,4)=nanmean(Cclim(:,:,10:12),3);

%%
str=["a","b","c","d","e","f","g","h","i","j","k","l","m","n","o","p"];
ti=["Autumn (MAM)","Winter (JJA)","Spring (SON)","Summer (DJF)"];
set(gcf,'Position',[300 400 1000 800]);
ha = tight_subplot(4,4,[.03 0.01],[.02 .04],[.05 0.06]);
for i=1:4
    axes(ha(i))
    m_proj('azimuthal equal-area','latitude',-90,'radius',55,'rectbox','on');
    hold on;
    m_pcolor(X,Y,ano(:,:,i));
    colortable = textread('/Users/jinfei/Documents/MATLAB/colorbar/cmp_flux.txt');
    colormap(ha(i),colortable);
    m_gshhs_i('color','k');
    m_grid('xaxislocation','top','xtick',[-120 -60 0 60 120 180],'ytick',[-90 -80 -70 -60],'fontsize',11);
    caxis([-15 15]);
    text(-0.6,0.57,str(i),'FontSize',16,'FontWeight','bold','fontname','Arial')
    tit=title(ti(i),'FontSize',14);
    set(tit,'position',get(tit,'Position') + [0 0.06 0])
end
hcb1 = colorbar;
title(hcb1,'hPa','FontSize',14)
set(hcb1,'position',[0.94 0.78 0.02 0.15],'YTick',[-15 0 15],'FontSize',12)
text(-6.1,-0.1,'SLP','FontSize',16,'rotation',90)

for i=5:8
    axes(ha(i))
    m_proj('azimuthal equal-area','latitude',-90,'radius',55,'rectbox','on');
    hold on;
    m_pcolor(X,Y,ano(:,:,i));
    colortable = textread('/Users/jinfei/Documents/MATLAB/colorbar/nrl_sirkes.txt');
    colormap(ha(i),colortable);
    m_gshhs_i('color','k');
    m_grid('xaxislocation','top','xtick',[-120 -60 0 60 120 180],'ytick',[-90 -80 -70 -60],'fontsize',11);
    caxis([-10 10]);
    text(-0.6,0.57,str(i),'FontSize',16,'FontWeight','bold','fontname','Arial')
end
hcb2 = colorbar;
title(hcb2,'^oC','FontSize',14)
set(hcb2,'position',[0.94 0.54 0.02 0.15],'YTick',[-10 0 10],'FontSize',12)
text(-6.1,-0.1,'T2m','FontSize',16,'rotation',90)

d = 18;dd = 0.6;
for i=9:12
    axes(ha(i))
    m_proj('azimuthal equal-area','latitude',-90,'radius',55,'rectbox','on');
    hold on;
    m_pcolor(X,Y,-ano(:,:,i));
    h1=m_quiver(lon(1:d:end,1:d:end),lat(1:d:end,1:d:end),uano(1:d:end,1:d:end,i-8)/dd,vano(1:d:end,1:d:end,i-8)/dd,0,'MaxHeadSize',3);
    set(h1,'color','k','linewidth',0.8);
    colortable = textread('/Users/jinfei/Documents/MATLAB/colorbar/ncl_default.txt');
    colormap(ha(i),colortable);
    m_gshhs_i('patch',[.8 .8 .8],'edgecolor',[.8 .8 .8],'MarkerSize',0.1);
    m_grid('xaxislocation','top','xtick',[-120 -60 0 60 120 180],'ytick',[-90 -80 -70 -60],'fontsize',11);
    caxis([-5e-5 5e-5]);
    text(-0.6,0.57,str(i),'FontSize',16,'FontWeight','bold','fontname','Arial')
end
hcb2 = colorbar;
title(hcb2,'K s^-^1','FontSize',14)
set(hcb2,'position',[0.94 0.3 0.02 0.15],'FontSize',12)
text(-6.1,-0.5,'Advection & {\itv}10','FontSize',16,'rotation',90)
h=m_quiver(0,-85,5/dd,0,0,'MaxHeadSize',3);
set(h,'color','k','linewidth',0.2);
m_text(-15,-80,'5 m/s','fontsize',14);

axes(ha(13))
m_proj('azimuthal equal-area','latitude',-90,'radius',55,'rectbox','on');
hold on;
m_gshhs_i('patch',[.8 .8 .8],'edgecolor',[.8 .8 .8],'MarkerSize',0.1);
m_pcolor(lon,lat,anoC(:,:,1));
m_contour(lon,lat,edgeC(:,:,1),[15 15],'color',sixteen2ten('#4E23D1')/255,'LineWidth',1.5)
colortable = textread('/Users/jinfei/Documents/MATLAB/colorbar/sunshine_diff_12lev.txt');
colormap(ha(13),colortable);
m_grid('xaxislocation','top','xtick',[-120 -60 0 60 120 180],'ytick',[-90 -80 -70 -60],'fontsize',11);
caxis([-50 50]);
text(-0.6,0.57,str(13),'FontSize',16,'FontWeight','bold','fontname','Arial')
    
axes(ha(14))
m_proj('azimuthal equal-area','latitude',-90,'radius',55,'rectbox','on');
hold on;
m_gshhs_i('patch',[.8 .8 .8],'edgecolor',[.8 .8 .8],'MarkerSize',0.1);
m_pcolor(lon,lat,anoC(:,:,2));
m_contour(lon,lat,edgeC(:,:,2),[15 15],'color',sixteen2ten('#4E23D1')/255,'LineWidth',1.5)
colortable = textread('/Users/jinfei/Documents/MATLAB/colorbar/sunshine_diff_12lev.txt');
colormap(ha(14),colortable);
m_grid('xaxislocation','top','xtick',[-120 -60 0 60 120 180],'ytick',[-90 -80 -70 -60],'fontsize',11);
caxis([-50 50]);
text(-0.6,0.57,str(14),'FontSize',16,'FontWeight','bold','fontname','Arial')
hcb3 = colorbar;
title(hcb3,'%','FontSize',14)
set(hcb3,'position',[0.94 0.06 0.02 0.15],'YTick',[-50 0 50],'FontSize',12)
% m_line([-20,-20],[-60,-55],'color','r','linewi',2);
% m_line([20,20],[-60,-55],'color','r','linewi',2);
% m_line(-20:20,-55:-55,'color','r','linewi',2);
% m_line(-20:20,-60:-60,'color','r','linewi',2);

axes(ha(15))
m_proj('azimuthal equal-area','latitude',-90,'radius',55,'rectbox','on');
hold on;
m_gshhs_i('patch',[.8 .8 .8],'edgecolor',[.8 .8 .8],'MarkerSize',0.1);
m_pcolor(lon,lat,anoC(:,:,3));
m_contour(lon,lat,edgeC(:,:,3),[15 15],'color',sixteen2ten('#4E23D1')/255,'LineWidth',1.5)
colortable = textread('/Users/jinfei/Documents/MATLAB/colorbar/sunshine_diff_12lev.txt');
colormap(ha(15),colortable);
m_grid('xaxislocation','top','xtick',[-120 -60 0 60 120 180],'ytick',[-90 -80 -70 -60],'fontsize',11);
caxis([-50 50]);
text(-0.6,0.57,str(15),'FontSize',16,'FontWeight','bold','fontname','Arial')
hcb3 = colorbar;
title(hcb3,'%','FontSize',14)
set(hcb3,'position',[0.94 0.06 0.02 0.15],'YTick',[-50 0 50],'FontSize',12)

axes(ha(16))
m_proj('azimuthal equal-area','latitude',-90,'radius',55,'rectbox','on');
hold on;
m_gshhs_i('patch',[.8 .8 .8],'edgecolor',[.8 .8 .8],'MarkerSize',0.1);
m_pcolor(lon,lat,anoC(:,:,4));
m_contour(lon,lat,edgeC(:,:,4),[15 15],'color',sixteen2ten('#4E23D1')/255,'LineWidth',1.5)
colortable = textread('/Users/jinfei/Documents/MATLAB/colorbar/sunshine_diff_12lev.txt');
colormap(ha(16),colortable);
m_grid('xaxislocation','top','xtick',[-120 -60 0 60 120 180],'ytick',[-90 -80 -70 -60],'fontsize',11);
caxis([-50 50]);
text(-0.6,0.57,str(16),'FontSize',16,'FontWeight','bold','fontname','Arial')
m_line([-120,-120],[-80,-70],'color','r','linewi',2);
m_line([-180,-180],[-80,-70],'color','r','linewi',2);
m_line(-180:-120,-70:-70,'color','r','linewi',2);
m_line(-180:-120,-80:-80,'color','r','linewi',2);
m_line([-120,-120],[-70,-60],'color','b','linewi',2);
m_line([-180,-180],[-70,-60],'color','b','linewi',2);
% m_line(-180:-120,-70:-70,'color','b','linewi',2);
m_line(-180:-120,-60:-60,'color','b','linewi',2);

text(-6.1,-0.1,'SIC','FontSize',16,'rotation',90)




