关于NCEP/NCAR再分析资料学习

在NECP/NCAR官网上下载再分析资料https://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.pressure.html

这里插一句，我突然发现NECP有一个在线绘图的网站，可以让我这样完全不懂的小白校对自己的程序和资料，也可以节省想用其资料分析的研究者的时间，可以先在上面点一点在线绘图，然后再决定是否要下载，确实非常实用。

我这个链接是绘制每日平均的资料，在网站右侧还有许多可以选的资料选项。

https://www.esrl.noaa.gov/psd/data/composites/day/

我选了每日的位势高度进行了下载，点那个see list。

因为我要每日的，所以继续点这个Daily Mean的See list。

有一堆年份可以选，我选了其中一个年份下载，大概80几M。

看下这个nc文件的属性，可以看到：一共17层，平面上是144*73的，下载的是每日的，所以time是365个时刻。

;本程序用于读取hgt500_1961_2016.nc  ;先抽查随机年份的夏季、冬季500hPa环流图 begin     fin=addfile("/home/honghx/date/NECPhgt/hgt500_1961_2016.nc","r")     hgt=fin->hgt     hgt=hgt*0.1     size=dimsizes(hgt)     ntime=size(0)     nlat=size(1);nlat=73     nlon=size(2);nlon=144     year=ispan(1961,2016,1)     nyear=dimsizes(year)     ; 做count用来记录每一年的第一天。     count=new(nyear+1,"integer")     count(0)=0     do i=1,nyear         if(year(i-1)%4.eq.0)then             count(i)=count(i-1)+366;下一年的起始值         else             count(i)=count(i-1)+365         end if     picture=new((/nlat,nlon/),"float")     picture=hgt(count(54)+115,:,:);2015年4月26日     ; wks=gsn_open_wks("png","January 1st 1961 500hPa")     wks=gsn_open_wks("png","April 26th 2015 500hPa")     gsn_define_colormap(wks,"")     res=True     res@cnLevelSelectionMode="ManualLevels"     res@cnMinLevelValF=500.0     res@cnMaxLevelValF=590.0     res@cnLevelSpacingF=5.0     res@mpCenterLonF=180;检查的时候才发现NECP文件是0-357.5的，差点全错！还好自己检查了！     res@cnFillOn=True     res@cnLinesOn=False     res@cnLineLabelsOn=False     res@mpOutlineOn=True     res@mpGeophysicalLineThicknessF=2.0 ;这两行是为了加粗边界和国界线     res@mpNationalLineThicknessF=2.0     res@tiMainString=("April 26th 500hPa")     plot=gsn_csm_contour_map(wks,picture,res)

end

我画了2015年4月26日的500hPa等压面的位势高度场。

而在线绘图网站上的图片如下：（一开始我画的和他不一样，我才发现我NECP/NCAR资料应该是0-357.5°的，而NCL自动画的是-180~0~+180的，是错的，还好检查了- -）

至此，NECP/NCAR再分析资料算是刚刚弄明白怎么下载及操作。