求助！grads里不同预报时间的提取

txdjltao

我主要是想转换grib2格式数据的格式。
我的数据本身情况是这样的，6hr降水集合预报（50个集合），每天进行接下来9天的预报，如此进行6个月。

我现在成功在ctl里添加了EDEF，也就是实现了集合信息的读取。
但是我就不太明白时间轴的问题。
tdef 773 linear 12Z30apr2007 6hr
这样不就只有一个轴吗？
那么set t 2 是12Z30apr2007起前6小时预报？
那么set t 5 是12Z30apr2007起前24小时预报吗？
从数值上来看貌似也不是这样……
以及假若是这样，我怎么查看12Z01may2007起的预报呢？
不知道是不是要修改ctl？还是怎么读取呢？

附上wgrib2解码的情况：
1:0:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=+1
2:206:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=-2
3:412:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=+3
4:618:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=-4
5:824:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=+5
6:1030:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=-6
7:1236:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=+7
8:1442:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=-8
9:1648:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=+9
10:1854:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=-10
11:2060:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=+11
12:2266:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=-12
13:2472:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=+13
14:2678:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=-14
15:2884:d=2007043012:TPRATE:surface:0-0 day acc fcst:ENS=+15
.......
51:10300:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=+1
52:11640:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=-2
53:12980:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=+3
54:14320:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=-4
55:15660:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=+5
56:17000:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=-6
57:18340:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=+7
58:19680:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=-8
59:21020:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=+9
60:22360:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=-10
61:23700:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=+11
62:25040:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=-12
63:26380:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=+13
64:27720:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=-14
65:29060:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=+15
66:30400:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=-16
67:31740:d=2007043012:TPRATE:surface:0-6 hour acc fcst:ENS=+17
.......
239:262220:d=2007043012:TPRATE:surface:0-1 day acc fcst:ENS=+39
240:263560:d=2007043012:TPRATE:surface:0-1 day acc fcst:ENS=-40
241:264900:d=2007043012:TPRATE:surface:0-1 day acc fcst:ENS=+41
242:266240:d=2007043012:TPRATE:surface:0-1 day acc fcst:ENS=-42
243:267580:d=2007043012:TPRATE:surface:0-1 day acc fcst:ENS=+43
244:268920:d=2007043012:TPRATE:surface:0-1 day acc fcst:ENS=-44
245:270260:d=2007043012:TPRATE:surface:0-1 day acc fcst:ENS=+45
246:271600:d=2007043012:TPRATE:surface:0-1 day acc fcst:ENS=-46
247:272940:d=2007043012:TPRATE:surface:0-1 day acc fcst:ENS=+47
248:274280:d=2007043012:TPRATE:surface:0-1 day acc fcst:ENS=-48
249:275620:d=2007043012:TPRATE:surface:0-1 day acc fcst:ENS=+49
250:276960:d=2007043012:TPRATE:surface:0-1 day acc fcst:ENS=-50
251:278300:d=2007043012:TPRATE:surface:0-30 hour acc fcst:ENS=+1
252:279640:d=2007043012:TPRATE:surface:0-30 hour acc fcst:ENS=-2
253:280980:d=2007043012:TPRATE:surface:0-30 hour acc fcst:ENS=+3
254:282320:d=2007043012:TPRATE:surface:0-30 hour acc fcst:ENS=-4
255:283660:d=2007043012:TPRATE:surface:0-30 hour acc fcst:ENS=+5
256:285000:d=2007043012:TPRATE:surface:0-30 hour acc fcst:ENS=-6
257:286340:d=2007043012:TPRATE:surface:0-30 hour acc fcst:ENS=+7
258:287680:d=2007043012:TPRATE:surface:0-30 hour acc fcst:ENS=-8
259:289020:d=2007043012:TPRATE:surface:0-30 hour acc fcst:ENS=+9
.......
842:2410240:d=2007043012:TPRATE:surface:0-9 day acc fcst:ENS=-42
1843:2411580:d=2007043012:TPRATE:surface:0-9 day acc fcst:ENS=+43
1844:2412920:d=2007043012:TPRATE:surface:0-9 day acc fcst:ENS=-44
1845:2414260:d=2007043012:TPRATE:surface:0-9 day acc fcst:ENS=+45
1846:2415600:d=2007043012:TPRATE:surface:0-9 day acc fcst:ENS=-46
1847:2416940:d=2007043012:TPRATE:surface:0-9 day acc fcst:ENS=+47
1848:2418280:d=2007043012:TPRATE:surface:0-9 day acc fcst:ENS=-48
1849:2419620:d=2007043012:TPRATE:surface:0-9 day acc fcst:ENS=+49
1850:2420960:d=2007043012:TPRATE:surface:0-9 day acc fcst:ENS=-50
1851:2422300:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=+1
1852:2422506:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=-2
1853:2422712:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=+3
1854:2422918:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=-4
1855:2423124:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=+5
1856:2423330:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=-6
1857:2423536:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=+7
1858:2423742:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=-8
1859:2423948:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=+9
1860:2424154:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=-10
1861:2424360:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=+11
1862:2424566:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=-12
1863:2424772:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=+13
1864:2424978:d=2007050112:TPRATE:surface:0-0 day acc fcst:ENS=-14
......

严重急切求助啊！！！拜托了！！
新手琢磨了好久了……好不容易走到这一步……

txdjltao

不知道解释清楚没有，如果不明白请大家问哦！！
急切求助！！！

txdjltao

附上ctl如下：（这个是可以在grads里正常打开的，只是我不明白那个时间轴……）
dset ^rjtd_sl_tp_pf_2007.grib
index ^rjtd_sl_tp_pf_2007.grib.idx
undef 9.999E+20
title rjtd_sl_tp_pf_2007.grib
*  produced by g2ctl v0.0.3d
* griddef=1:0:(27 x 21):grid_template=0:winds(N/S): lat-lon grid:(27 x 21) units 1e-06 input WE:NS output WE:SN res 48 lat 38.000000 to 28.000000 by 0.500000 lon 108.000000 to 121.000000 by 0.500000 #points=567

dtype grib2
ydef 21 linear 28.000000 0.5
xdef 27 linear 108.000000 0.500000
tdef 773 linear 12Z30apr2007 6hr
edef 50
ens01 773 12Z30apr2007 3,1
ens02 773 12Z30apr2007 2,2
ens03 773 12Z30apr2007 3,3
ens04 773 12Z30apr2007 2,4
ens05 773 12Z30apr2007 3,5
ens06 773 12Z30apr2007 2,6
ens07 773 12Z30apr2007 3,7
ens08 773 12Z30apr2007 2,8
ens09 773 12Z30apr2007 3,9
ens10 773 12Z30apr2007 2,10
ens11 773 12Z30apr2007 3,11
ens12 773 12Z30apr2007 2,12
ens13 773 12Z30apr2007 3,13
......
ens49 773 12Z30apr2007 3,49
ens50 773 12Z30apr2007 2,50
ENDEDEF
zdef 1 linear 1 1
vars 1
TPRATEsfc 0,1   0,1,52,1 ** surface acc Total Precipitation Rate [kg/m^2/s]
ENDVARS

风行浪尖

虽然不懂，帮忙顶下

linikiku

我也想知道。。

覃小覃

ctl时用哪个命令生成的哦亲？g2ctl -ens ...?

rhythm9806

请问一下，你如何在ctl文件里面添加集合预报数据的成员？

rhythm9806

覃小覃 发表于 2014-11-16 13:30
ctl时用哪个命令生成的哦亲？g2ctl -ens ...?

-ens 使用的时候出错呢，你怎么解决的？

xiaoyudian

同问，ctl如何添加集合预报成员信息

chenshengjie

求教。我从TIGGE网站下载的EC降水集合预报，数据格式是grib，如何利用grads读取呢？非常感谢