month=1:168;
lon=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/SST/dat/air.mon.mean.nc','lon');
lat=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/SST/dat/air.mon.mean.nc','lat');
load ~/Aerosol_Variaiton/IRUN/SST/ncfile/NINO34_monthly_NCEP.mat

windowSize = 5;
x=NINO34;
y=filter(ones(1,windowSize)/windowSize,1,x);

EN_Month=zeros(length(month),1);
LL_Month=zeros(length(month),1);

for monthID=1:length(month);
    if y(monthID)>0.8
        EN_Month(monthID)=1;
    end
    if y(monthID)<-0.8
        LL_Month(monthID)=1;
    end
end


save('ENSO_Month_NCEP_extreme.mat','EN_Month','LL_Month')

nccreate('ENSO_Month_NCEP_extreme.nc','Month','Dimensions',{'time',Inf},'Format','classic')
ncwrite('ENSO_Month_NCEP_extreme.nc','Month', month)
ncwriteatt('ENSO_Month_NCEP_extreme.nc','Month','long_name','month')
nccreate('ENSO_Month_NCEP_extreme.nc','EN_Month','Dimensions',{'month',length(month)},'Format','classic')
ncwrite('ENSO_Month_NCEP_extreme.nc','EN_Month', EN_Month)
nccreate('ENSO_Month_NCEP_extreme.nc','LL_Month','Dimensions',{'month',length(month)},'Format','classic')
ncwrite('ENSO_Month_NCEP_extreme.nc','LL_Month', LL_Month)
ncdisp('ENSO_Month_NCEP_extreme.nc')