GFS GRIB2数据处理后没有数值？

雨竹晓风

在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

写成gfs.t12z.pgrb2.0p25.f0%h2试试

elegisite

楼主的问题解决了么？我到第四步生成idx文件的时候就出错了，显示打不开描述文件，不知道楼主之前有没有遇到过这个问题？

QQtangxian

请问楼主最后怎么解决的呢改写成f0%h2？

雨竹晓风

最后改成F%f3就行了。

QQtangxian

我也遇到相同问题请问怎么解决的呢

QQtangxian

elegisite 发表于 2015-7-30 18:26
楼主的问题解决了么？我到第四步生成idx文件的时候就出错了，显示打不开描述文件，不知道楼主之前有没有遇 ...

请问您的问题解决了么

雨竹晓风

gribmap的时候用参数 gribmap -v -i  

风卷云

请问楼主问题解决了没

飘渺居

请问你知道生成idx文件时后面的 -0 是啥意思吗？