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');

WIND_daily=zeros(90,21);
WIND_seasonal=zeros(21,1);

WIND_seasonal(:)=mean(mean(mean(WIND(117:122,57:68,:,:),3),2),1);
WIND_seasonal(13)=(WIND_seasonal(11)+WIND_seasonal(12))/2;

WIND_daily(:,:)=mean(mean(WIND(117:122,57:68,:,:),2),1);