load ~/NC/grid.mat
load ~/Aerosol_Variaiton/IRUN/SST/ncfile/SST.mat
load ~/NC/OCNFRAC.mat
year=1:150;

PDO=zeros(150,1);
PDO(:)=mean(mean(SST_ANN(47:99,59:80,:),2),1);


save('SST_PDO.mat','PDO')



nccreate('SST_PDO.nc','PDO','Dimensions',{'year',length(year)},'Format','classic')
ncwrite('SST_PDO.nc','PDO', PDO)
ncdisp('SST_PDO.nc')
