lon=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/UWND/dat/uwnd.mon.mean.nc','lon');
lat=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/UWND/dat/uwnd.mon.mean.nc','lat');
year=1980:2015;
WIND=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/WIND/ncfile/WIND_NCEP.nc','WIND_DJF');
[~,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,length(year)-35:length(year)),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,length(year)-35:length(year)),2),1);