year=1981:2007;

temp=ncread('/Users/yangyang/Aerosol_Variaiton/VIS/ncfile/VIS_1981_2015_NC.nc','vis_ndays_average');
Haze_NC=temp(1:27);
temp=ncread('/Users/yangyang/Aerosol_Variaiton/VIS/ncfile/VIS_1981_2015_SC.nc','vis_ndays_average');
Haze_SC=temp(1:27);
DSFI=[0.1,-0.05,0.5,-0.1,-0.1,-0.3,-0.25,-0.6,-0.8,-0.7,-1.1,-0.85,-0.95,-0.8,-0.9,-1.25,-1.4,-0.9,-1.05,-1.2,-0.7,-1.1,-1.5,-1.4,-1.7,-1.25,-1.2];



WIND_DJF=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/WIND/ncfile/WIND_NCEP.nc','WIND_DJF');
WIND_MAM=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/WIND/ncfile/WIND_NCEP.nc','WIND_MAM');

lon=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/WIND/ncfile/WIND_NCEP.nc','lon');
lat=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/WIND/ncfile/WIND_NCEP.nc','lat');
[~,lonmin]=min(abs(lon-100));
[~,lonmax]=min(abs(lon-110));
[~,latmax]=min(abs(lat-35));
[~,latmin]=min(abs(lat-45));
WIND_ECHINA_DJF=zeros(length(year),1);
WIND_ECHINA_DJF(:)=mean(mean(WIND_DJF(lonmin:lonmax,latmin:latmax,:),2),1);
WIND_ECHINA_MAM=zeros(length(year),1);
WIND_ECHINA_MAM(:)=mean(mean(WIND_MAM(lonmin:lonmax,latmin:latmax,:),2),1);

corrcoef(WIND_ECHINA_DJF,WIND_ECHINA_MAM)