如何正确编写TIGGE-EC资料的Vtable

一顾留白

目前想用TIGGE-EC 00时刻的资料跑WRF,卡在了WPS的ungrib上，因为没有合适对应的Vtable，所以就想自己写一个，之前有参考其他大佬运行欧洲中心集合数据的帖子（http://bbs.06climate.com/forum.p ... 62&highlight=Vtable）和（http://bbs.06climate.com/forum.php?mod=viewthread&tid=48433），但是发现其中的Vtable并不适用于TIGGE-EC资料，我用g2print.exe输出了一下变量，不知道哪里出了问题，请各位大佬给看一下，感谢

ljchen1989

要不你po一个数据上来看看？

一顾留白

抱歉啊数据有点大，我上传了下载数据的python脚本，如果不方便的话您留个我发过去也行~~麻烦啦
另外我看之前的帖子只用了TIGGE-EC的高空资料，地面资料用的fnl的，所以我也只用了TIGGE-EC的pl资料，不过sfc的脚本我也上传啦~·

一顾留白

抱歉，忘了脚本还要用金钱下，下面是TIGGE_EC-pl.py
#!/usr/bin/env python
from ecmwfapi import ECMWFDataServer
server = ECMWFDataServer()
server.retrieve({
    "class": "ti",
    "dataset": "tigge",
    "date": "2019-04-23/to/2019-04-25",
    "expver": "prod",
    "grid": "0.5/0.5",
    "levelist": "200/250/300/500/700/850/925/1000",
    "levtype": "pl",
    "origin": "ecmf",
    "param": "130/131/132/133/156",
    "step": "0",
    "time": "00:00:00/12:00:00",
    "type": "cf",
    "target": "D:/2019042325pl.grb2",
})
TIGGE_EC-sfc.py
#!/usr/bin/env python
from ecmwfapi import ECMWFDataServer
server = ECMWFDataServer()
server.retrieve({
    "class": "ti",
    "dataset": "tigge",
    "date": "2019-04-23/to/2019-04-25",
    "expver": "prod",
    "grid": "0.5/0.5",
    "levtype": "sfc",
    "origin": "ecmf",
    "param": "59/134/136/146/147/151/165/166/167/168/172/176/177/179/189/235/228001/228002/228039/228139/228141/228144/228164/228228",
    "step": "0",
    "time": "00:00:00/12:00:00",
    "type": "fc",
    "target": "D:/20190423-25sfc.grib2",
})

一顾留白

自己填下坑，vtable没有问题，是数据的问题，虽然有6个时次但是读的时候只会显示第一个时次，所以会导致数据缺失的错误，但是我分别每个时次ungrib的话，metgrid还是会出错，错误如下：
Processing domain 1 of 1
Processing 2019-04-23_00
    SFC
    PRESS
Processing 2019-04-23_12
    SFC
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image              PC                Routine            Line        Source            
metgrid.exe        0000000000456617  Unknown               Unknown  Unknown
metgrid.exe        000000000045279F  Unknown               Unknown  Unknown
metgrid.exe        0000000000418576  Unknown               Unknown  Unknown
metgrid.exe        0000000000410A5C  Unknown               Unknown  Unknown
libc.so.6          0000003BD661ECDD  Unknown               Unknown  Unknown
metgrid.exe        0000000000410959  Unknown               Unknown  Unknown
继续探索中~

ljchen1989

建议你要不在python或者shell里面写循环，把每个时次的grib文件单独写出来，要不用wgrib2把你下的文件拆开
话说，你看过ungrib的log吗？确认变量都ok嘛？

一顾留白

我现在是单个时次的下的，ungrib.log也看过的，不是Vtable的问题，是这种资料几个时次一块下只能显示初始时次，所以显示资料缺失，我分别下好各个时次的统一link就可以了~~·谢谢啦~~~

解脱

你下载的数据时间和你模式开始运行时间不同呀。模式开始时间是4月23日12时的，但是你下载的是00时的数据。提示里也说了没有找到4月23日12时的数据，是不是你配置参数有问题呀？！

                               
登录/注册后可看大图

一顾留白

是数据的问题，我的数据是23日00时到25日00时，但是需要从23日12时开始转，WPS只能识别初始时刻，也就是说WPS只能识别到这是23日00时的数据。我看ERA5好像也有这种问题，目前解决啦~~谢谢~~~