% show header ncdisp('/data/icdc/reanalyses/ncep_reanalysis1/DATA/2m_airtemp_monthly/air2m.mon.mean.nc') % get temperature variable at time step 1 air=ncread('/data/icdc/reanalyses/ncep_reanalysis1/DATA/2m_airtemp_monthly/air2m.mon.mean.nc', 'air',[1 1 804], [192 94 1]); % get lat and lon lat=ncread('/data/icdc/reanalyses/ncep_reanalysis1/DATA/2m_airtemp_monthly/air2m.mon.mean.nc', 'lat',[1], [94]); lon=ncread('/data/icdc/reanalyses/ncep_reanalysis1/DATA/2m_airtemp_monthly/air2m.mon.mean.nc', 'lon',[1], [192]); % set missing value to NaN air(air<-1000)=NaN % plot map imagesc(lon,lat,air') axis xy