year=1981:2015;
year=1:67;

%VIS=ncread('/Users/yangyang/Aerosol_Variaiton/VIS/ncfile/VIS_1981_2015.nc','vis_mean_average');
%HAZE=ncread('/Users/yangyang/Aerosol_Variaiton/VIS/ncfile/VIS_1981_2015.nc','vis_ndays_average');

%WIND=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/WIND/ncfile/WIND_ECMWF_1979_2014.nc','WIND_DJF');
%lon=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/UWND/dat/wind850_ECMWF.nc','longitude');
%lat=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/UWND/dat/wind850_ECMWF.nc','latitude');
%[~,lonmin]=min(abs(lon-110));
%[~,lonmax]=min(abs(lon-122.5));
%[~,latmax]=min(abs(lat-20));
%[~,latmin]=min(abs(lat-45));
%WIND_ECHINA=zeros(length(year),1);
%WIND_ECHINA(:)=mean(mean(WIND(lonmin:lonmax,latmin:latmax,2:36),2),1);


%[~,lonmin]=min(abs(lon-100));
%[~,lonmax]=min(abs(lon-110));
%[~,latmax]=min(abs(lat-35));
%[~,latmin]=min(abs(lat-45));
%WIND_GOBI=zeros(length(year),1);
%WIND_GOBI(:)=mean(mean(WIND(lonmin:lonmax,latmin:latmax,2:36),2),1);


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)