这⾥插⼀句，我突然发现NECP有⼀个在线绘图的⽹站，可以让我这样完全不懂的⼩⽩校对⾃⼰的程序和资料，也可以节省想⽤其资料分析的研究者的时间，可以先在上⾯点⼀点在线绘图，然后再决定是否要下载，确实⾮常实⽤。 我这个链接是绘制每⽇平均的资料，在⽹站右侧还有许多可以选的资料选项。
我选了每⽇的位势⾼度进⾏了下载，点那个see list。 因为我要每⽇的，所以继续点这个Daily Mean的See list。
有⼀堆年份可以选，我选了其中⼀个年份下载，⼤概80⼏M。
看下这个nc⽂件的属性，可以看到：⼀共17层，平⾯上是144*73的，下载的是每⽇的，所以time是365个时刻。
；本程序⽤于读取hgt500_
；先抽查随机年份的夏季、冬季500hPa环流图
begin
fin=addfile("/home/honghx/date/NECPhgt/hgt500_","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)
workbench; 做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")dopodo
plot=gsn_csm_contour_map(wks,picture,res)
end
我画了2015年4⽉26⽇的500hPa等压⾯的位势⾼度场。
2014年sci影响因子⽽在线绘图⽹站上的图⽚如下：（⼀开始我画的和他不⼀样，我才发现我NECP/NCAR资料应该是0-357.5°的，⽽NCL⾃动画的是-180~0~+180的，是错的，还好检查了- -）
⾄此，NECP/NCAR再分析资料算是刚刚弄明⽩怎么下载及操作。
