关于自己写Vtable

林子木 · 发表于 2016-3-13 11:33:15

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

我最近在用欧洲的集合预报跑WRF，因为没有合适的Vtable所以自己按照wrf手册上写的，但是手册上只是介绍了关于GRIB1的资料的Vtable的写法，现在的资料都是GRIB2的，在解码上两者还是有差别的，我用g2print.exe 查看资料后给出信息如下：
rec Prod Cat Param  Lvl Lvl    Lvl    Prod Name          Time
   Fcst
  num Disc    num code one    two    Templ
   hour
-------------------------------------------------------------------------------
--------
1 0 0 0    100  100000    0    1    TMP    2013-08-12_00:00:00 06  PDT4.1
2 0 2 2    100  100000    0    1    UGRD    2013-08-12_00:00:00 06  PDT4.1
3 0 2 3    100  100000    0    1    VGRD    2013-08-12_00:00:00 06  PDT4.1
4 0 1 0    100  100000    0    1    SPFH    2013-08-12_00:00:00 06  PDT4.1
5 0 3 5    100  100000    0    1    HGT    2013-08-12_00:00:00 06  PDT4.1
6 0 0 0    100 92500    0    1    TMP    2013-08-12_00:00:00 06  PDT4.1
7 0 2 2    100 92500    0    1    UGRD    2013-08-12_00:00:00 06  PDT4.1
8 0 2 3    100 92500    0    1    VGRD    2013-08-12_00:00:00 06  PDT4.1
9 0 1 0    100 92500    0    1    SPFH    2013-08-12_00:00:00 06  PDT4.1
  10 0 3 5    100 92500    0    1    HGT    2013-08-12_00:00:00 06  PDT4.1
  11 0 0 0    100 85000    0    1    TMP    2013-08-12_00:00:00 06  PDT4.1
  12 0 2 2    100 85000    0    1    UGRD    2013-08-12_00:00:00 06  PDT4.1
  13 0 2 3    100 85000    0    1    VGRD    2013-08-12_00:00:00 06  PDT4.1
  14 0 1 0    100 85000    0    1    SPFH    2013-08-12_00:00:00 06  PDT4.1
  15 0 3 5    100 85000    0    1    HGT    2013-08-12_00:00:00 06  PDT4.1
  16 0 0 0    100 70000    0    1    TMP    2013-08-12_00:00:00 06  PDT4.1
  17 0 2 2    100 70000    0    1    UGRD    2013-08-12_00:00:00 06  PDT4.1
  18 0 2 3    100 70000    0    1    VGRD    2013-08-12_00:00:00 06  PDT4.1
  19 0 1 0    100 70000    0    1    SPFH    2013-08-12_00:00:00 06  PDT4.1
  20 0 3 5    100 70000    0    1    HGT    2013-08-12_00:00:00 06  PDT4.1
  21 0 0 0    100 50000    0    1    TMP    2013-08-12_00:00:00 06  PDT4.1
  22 0 2 2    100 50000    0    1    UGRD    2013-08-12_00:00:00 06  PDT4.1
  23 0 2 3    100 50000    0    1    VGRD    2013-08-12_00:00:00 06  PDT4.1
  24 0 1 0    100 50000    0    1    SPFH    2013-08-12_00:00:00 06  PDT4.1
  25 0 3 5    100 50000    0    1    HGT    2013-08-12_00:00:00 06  PDT4.1
  26 0 0 0    100 30000    0    1    TMP    2013-08-12_00:00:00 06  PDT4.1
  27 0 2 2    100 30000    0    1    UGRD    2013-08-12_00:00:00 06  PDT4.1
  28 0 2 3    100 30000    0    1    VGRD    2013-08-12_00:00:00 06  PDT4.1
  29 0 1 0    100 30000    0    1    SPFH    2013-08-12_00:00:00 06  PDT4.1
  30 0 3 5    100 30000    0    1    HGT    2013-08-12_00:00:00 06  PDT4.1
  31 0 0 0    100 25000    0    1    TMP    2013-08-12_00:00:00 06  PDT4.1
  32 0 2 2    100 25000    0    1    UGRD    2013-08-12_00:00:00 06  PDT4.1
  33 0 2 3    100 25000    0    1    VGRD    2013-08-12_00:00:00 06  PDT4.1
  34 0 1 0    100 25000    0    1    SPFH    2013-08-12_00:00:00 06  PDT4.1
  35 0 3 5    100 25000    0    1    HGT    2013-08-12_00:00:00 06  PDT4.1
  36 0 0 0    100 20000    0    1    TMP    2013-08-12_00:00:00 06  PDT4.1
  37 0 2 2    100 20000    0    1    UGRD    2013-08-12_00:00:00 06  PDT4.1
  38 0 2 3    100 20000    0    1    VGRD    2013-08-12_00:00:00 06  PDT4.1
  39 0 1 0    100 20000    0    1    SPFH    2013-08-12_00:00:00 06  PDT4.1
  40 0 3 5    100 20000    0    1    HGT    2013-08-12_00:00:00 06  PDT4.1
  41 0 3 5    100 5000    0    1    HGT    2013-08-12_00:00:00 06  PDT4.1

Successful completion of g2print
然后自己写的Vtable这样
GRIB | Level| Level| Level| metgrid  |  metgrid | metgrid                               |GRIB2|GRIB2|GRIB2|GRIB2|GB2|
Code | Code | 1  | 2  | Name    |  Units | Description                            |Discp|Catgy|Param|Level|PDT|
-----+------+------+------+----------+----------+------------------------------------------+---------------------------+
  5  | 100  | *  |    | HGT    | m       | Height                                  |  0  |  2  |  5  | 100 | 1 |
  0  | 100  | *  |    | TT    | K       | Temperature                            |  0  |  2  |  0  | 100 | 1 |
  2  | 100  | *  |    | UU    | m s-1 | U                                     |  0  |  1  |  2  | 100 | 1 |
  3  | 100  | *  |    | VV    | m s-1 | V                                     |  0  |  3  |  3  | 100 | 1 |
  1  | 100  | *  |    | SPECHUMD | kg kg-1  | Specific humidity                      |  0  |  1  |  1  | 100 | 1 |
-----+------+------+------+----------+----------+------------------------------------------+-----------------------+---+

然后。。。。。。。就一个变量都解不出来，我觉得主要是后面关于GRIB2数据描述的那一部分不准确，但也不知道怎么写，有大神指点一二吗

lukisy · 发表于 2016-3-13 14:45:58

虽然不是很清楚,还是支持下哈 ,同求解

ruilier2014 · 发表于 2017-1-6 12:03:23

请问后来解决了吗？

ljchen1989 · 发表于 2017-1-6 12:09:36

楼楼都能编译出g2print.exe了，为啥不直接用wgrib2呢？
rewrite grib2 data to binary format
and modify the wps intermediate file
it should be easier than creating vtable

林子木 · 发表于 2017-2-24 10:19:16

ljchen1989 发表于 2017-1-6 12:09
楼楼都能编译出g2print.exe了，为啥不直接用wgrib2呢？
rewrite grib2 data to binary format
and modify ...

受教了，我不知道还有这种方法，以后请大神多赐教

林子木 · 发表于 2017-2-24 10:20:09

ruilier2014 发表于 2017-1-6 12:03
请问后来解决了吗？

解决了，又写了一个帖子，把Vtable放上去了。

小其其格 · 发表于 2017-8-14 06:28:00

厉害！

kylin · 发表于 2018-5-15 16:44:53

楼主怎么解决的，能否附下链接，还有wgrib2这个是怎么用的呢，求指教呀

qimukakaxi · 发表于 2018-6-12 12:47:37

学习学习，了解一下

梁小梁 · 发表于 2019-3-24 10:38:00

学习了，感谢楼主！！！

		自动登录	找回密码
密码			立即注册