关于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再分析资料算是刚刚弄明白怎么下载及操作。