load ~/NC/grid.mat
load ~/Aerosol_Variaiton/IRUN/SST/ncfile/SST_EOF.mat
load ~/Aerosol_Variaiton/IRUN/RH/ncfile/RH_region.mat
PC_SST=-PC_SST;
year=1:150;

[r, p_SPO]=corrcoef(RH_SPO,PC_SST);
corr_SPO=r(1,2);

[r, p_SAO]=corrcoef(RH_SAO,PC_SST);
corr_SAO=r(1,2);

[r, p_SIO]=corrcoef(RH_SIO,PC_SST);
corr_SIO=r(1,2);

[r, p_TPO]=corrcoef(RH_TPO,PC_SST);
corr_TPO=r(1,2);

[r, p_TAO]=corrcoef(RH_TAO,PC_SST);
corr_TAO=r(1,2);

[r, p_TIO]=corrcoef(RH_TIO,PC_SST);
corr_TIO=r(1,2);

[r, p_NPO]=corrcoef(RH_NPO,PC_SST);
corr_NPO=r(1,2);

[r, p_NAO]=corrcoef(RH_NAO,PC_SST);
corr_NAO=r(1,2);

[r, p_SAM]=corrcoef(RH_SAM,PC_SST);
corr_SAM=r(1,2);

[r, p_AFR]=corrcoef(RH_AFR,PC_SST);
corr_AFR=r(1,2);

[r, p_OCE]=corrcoef(RH_OCE,PC_SST);
corr_OCE=r(1,2);

[r, p_NAM]=corrcoef(RH_NAM,PC_SST);
corr_NAM=r(1,2);

[r, p_EUR]=corrcoef(RH_EUR,PC_SST);
corr_EUR=r(1,2);

[r, p_ASI]=corrcoef(RH_ASI,PC_SST);
corr_ASI=r(1,2);