year=1986:2006;
lon=-180:2.5:177.5;
lat=-90:2:90;

WIND=ncread('/Users/yangyang/Extreme_Aerosol/Meteorology/WIND850SPEED/ncfile/WIND850SPEED_daily.nc','wind850speed');
PM25=ncread('/Users/yangyang/Extreme_Aerosol/MET/ncfile/Aerosol_daily.nc','pm25');

WIND_ECHINA=zeros(21,1);
WIND_ECHINA(:)=mean(mean(mean(WIND(117:122,57:68,:,:),3),2),1);

WIND_ECHINA(13)=(WIND_ECHINA(11)+WIND_ECHINA(12))/2;

PM25_ECHINA=zeros(21,1);
PM25_ECHINA(:)=mean(mean(mean(PM25(117:122,57:68,:,:),3),2),1);
PM25_ECHINA(13)=(PM25_ECHINA(11)+PM25_ECHINA(12))/2;


Y=PM25_ECHINA;
X=[ones(length(year),1),WIND_ECHINA];
[b,~,~,~,stats]=regress(Y,X);
Slope=b(2);
R2=stats(1);