load ~/NC/grid
year=1:150;
WIND=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/WIND/ncfile/WIND.nc','WIND_DJF');

WIND_ECHINA=zeros(length(year),1);
WIND_ECHINA(:)=mean(mean(WIND(45:50,60:72,:),2),1);


WIND_GOBI=zeros(length(year),1);
WIND_GOBI(:)=mean(mean(WIND(41:45,67:73,:),2),1);


DUSTSF=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/DUSTSF/ncfile/DUSTSF.nc','DUSTSF_DJF');
DUSTSF_GOBI=mean(mean(DUSTSF(41:45,67:73,:),2),1)*60*60*24*90*1000/1e12;