year=1979:2014;
lon=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/RH/dat/RHum.mon.mean.nc','lon');
lat=ncread('/Users/yangyang/Aerosol_Variaiton/IRUN/RH/dat/RHum.mon.mean.nc','lat');
load ~/Aerosol_Variaiton/IRUN/RH/ncfile/RH_NCEP.mat

RH_EN_temp=RH_ANN(:,:,[4 9 13 14 15 16 19 24 26]);
RH_LL_temp=RH_ANN(:,:,[10 18 20 22 30]);

RH_EN=mean(RH_EN_temp,3);
RH_LL=mean(RH_LL_temp,3);

Diff_RH=RH_EN-RH_LL;

save('Diff_RH_NCEP_ENSO2.mat','RH_EN','RH_LL','Diff_RH')

nccreate('Diff_RH_NCEP_ENSO2.nc','lat','Dimensions',{'lat',73},'Format','classic')
ncwrite('Diff_RH_NCEP_ENSO2.nc','lat', lat)
ncwriteatt('Diff_RH_NCEP_ENSO2.nc','lat','units','degrees_north')
ncwriteatt('Diff_RH_NCEP_ENSO2.nc','lat','long_name','Latitude')
nccreate('Diff_RH_NCEP_ENSO2.nc','lon','Dimensions',{'lon',144},'Format','classic')
ncwrite('Diff_RH_NCEP_ENSO2.nc','lon', lon)
ncwriteatt('Diff_RH_NCEP_ENSO2.nc','lon','units','degrees_east')
ncwriteatt('Diff_RH_NCEP_ENSO2.nc','lon','long_name','Longitude')
nccreate('Diff_RH_NCEP_ENSO2.nc','RH_EN','Dimensions',{'lon',144,'lat',73},'Format','classic')
ncwrite('Diff_RH_NCEP_ENSO2.nc','RH_EN', RH_EN)
nccreate('Diff_RH_NCEP_ENSO2.nc','RH_LL','Dimensions',{'lon',144,'lat',73},'Format','classic')
ncwrite('Diff_RH_NCEP_ENSO2.nc','RH_LL', RH_LL)
nccreate('Diff_RH_NCEP_ENSO2.nc','Diff_RH','Dimensions',{'lon',144,'lat',73},'Format','classic')
ncwrite('Diff_RH_NCEP_ENSO2.nc','Diff_RH', Diff_RH)
ncdisp('Diff_RH_NCEP_ENSO2.nc')