GRIB2数据 Cannot contour grid - all undefined values

雨竹晓风

在ftp://ftp.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.2015051812/批量下载了gfs.t12z.pgrb2.0p25.f000......gfs.t12z.pgrb2.0p25.f048的文件，给g2ctl生成ctl文件，改写ctl后，再用gribmap生成idx，编写gs文件，运行grads后，出现

Cannot contour grid - all undefined values错误，无法出图。

给出具体操作，求大家帮助！
一、批量下载gfs.t12z.pgrb2.0p25.f000......gfs.t12z.pgrb2.0p25.f048文件。
二、g2ctl.exe d:\gfs\gfs.t00z.pgrbf00.grib2  > d:\gfs\gfs.ctl 生成ctl文件；
三、改写gfs1.ctl文件：
dset D:\gfs\2015051820\gfs.t12z.pgrb2.0p25.f%h3
options template
index D:\gfs\2015051820\gfs.idx
undef 9.999E+20
title D:\gfs\2015051820\gfs.t12z.pgrb2.0p25.f000
*  produced by g2ctl v0.0.3d
* griddef=1:0:(1440 x 721):grid_template=0: lat-lon grid:(1440 x 721) units 1e-06 input WE:NS output WE:SN res 48 lat 90.000000 to -90.000000 by 0.250000 lon 0.000000 to 359.750000 by 0.250000 #points=1038240
dtype grib2
ydef 721 linear -90.000000 0.25
xdef 1440 linear 0.000000 0.250000
tdef 17 linear 12Z18may2015 3hr
*  z has 26 levels -prs
zdef 26 levels 100000 97500 95000 92500 90000 85000 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 7000 5000 3000 2000 1000
options pascals
vars 121
变量略去......
四、 生成idx文件：gribmap  -v -i d:\gfs\gfs.ctl -0
五、编写gs文件：
      因为是多时次的，所以做单站要素时序图。
      'reinit'
'open D:\gfs\2015051820\gfs.ctl'
'set grads off'
'set grid on'
'set parea 1 10 1.5 7'
'set ylevs 1000 975 950 925 850 700 600 500 200 '
t=1
'q time'
yy=substr(result,16,4)
mm=substr(result,13,3)
dd=substr(result,11,2)
hh=substr(result,8,2)
'set lon 105'
'set lat 35'
'set lev 925 200'
'set t 2 17'
'set gxout shaded'
'set cint 10'
'set cmin 60'
'set clevs 60 70 80 90 '
'set rbcols -1 34 35 37 39 '
'd rhprs'
'cbarn 1 0'
'set gxout barb'
'set ccolor 1'
'set cthick 6'
'd ugrdprs*2.5;vgrdprs*2.5'
'set gxout contour'
'set cmax -0.1'
'set ccolor 2'
'd vvelprs'
'draw title RH(%) & omega(Pa/s) & wind'
filen='d:\gfs\cross-rh-omega-'%mm%dd%hh%'.png'
'draw string 9.1 0.7 forcasted by GFS'
'draw string 9.1 0.4 Initial time:'yy'.'mm'.'dd'.'hh''
'printim 'filen' white'
'disable print'
'c'
;
六、运行gs，结果就出现

Cannot contour grid - all undefined values错误。

      但是如果不是多时次处理，而是单时次处理运行的话，可以出图。

     求教大家给出帮助，

     dset D:\gfs\2015051820\gfs.t12z.pgrb2.0p25.f%h3里的gfs.t12z.pgrb2.0p25.f000..f048，应该用什么表示是多时次的？参考过别人的关于修改对应的文件名，采用通配符方式，文件名用％表示。其中：   
      ％y2   代表两位数年   
      ％y4   代表四位数年   
      ％m1   代表一位或者两位数的月   
      ％m2   代表两位数月（用0补齐1位数）  
      ％mc   3个字符月份的缩写   
      ％d1   1或2位天   
      ％d2   两位天   
      ％h1   1或者2位时   
      ％h2   2位时。

    麻烦大家帮助！！！谢谢！

鲸鱼便子

有两种方法：
1.把grads换成1.9的
2.把风场、湿度场、和位势高度场变量后面的解释缩短，就可以了；
具体啥原因还真不知道

雨竹晓风

@鲸鱼便子
RHprs   26,100  0,1,1 ** (1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10) none Relative Humidity [%]  这个是不是正常？

鲸鱼便子

你试试，我记不太清楚了，不过风场和位势高度都不行
你把后面的气压删掉一些试一下 ，变成这样
RHprs   26,100  0,1,1 ** (1000 -10) none Relative Humidity [%]

雨竹晓风

改了还是不行！