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_ECHINA_10th=prctile(WIND_ECHINA,10);

WIND_ECHINA_weak=0;
WIND_year_ECHINA=zeros(150,1);
for yearID=1:length(year)
    if WIND_ECHINA(yearID)<WIND_ECHINA_10th
        WIND_year_ECHINA(yearID)=1;
        WIND_ECHINA_weak=WIND_ECHINA_weak+WIND_ECHINA(yearID)/15;
    end
end

WIND_ECHINA_aveg=mean(WIND_ECHINA);
WIND_ECHINA_diff=WIND_ECHINA_weak-WIND_ECHINA_aveg;



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

WIND_GOBI_10th=prctile(WIND_GOBI,10);

WIND_GOBI_weak=0;
WIND_year_GOBI=zeros(150,1);
for yearID=1:length(year)
    if WIND_GOBI(yearID)<WIND_GOBI_10th
        WIND_year_GOBI(yearID)=1;
        WIND_GOBI_weak=WIND_GOBI_weak+WIND_GOBI(yearID)/15;
    end
end

WIND_GOBI_aveg=mean(WIND_GOBI);
WIND_GOBI_diff=WIND_GOBI_weak-WIND_GOBI_aveg;


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;