GRIB 1 到 2 的转换

grib_set -s edition=2 in.grib1 out.grib2

上述命令会转换如下内容：

• time
• geography
• vertical
• parameter
• local
• data

因为两种编码标准各有不同，再加上一些气象中心会使用本地码表，所以参数转换比较复杂。

转换会基于变量参数的唯一标识符 paramId。

基于 paramId 的转换

如何生成 2 米温度的 GRIB 文件：

grib_set -s paramId=165 in.grib1 out.grib1
grib_set -s paramId=165 in.grib2 out.grib2

如何将 GRIB 1 文件转码为 GRIB 2 文件：

grib_set -s edition=2 in.grib1 out.grib2

在转换到版本 2 的过程中，ecCodes 从 GRIB 1 拷贝 paramId 到 GRIB 2：

1. 从 GRIB 1 中获取 paramId（=165）
2. 修改版本号为 2，生成 GRIB 2 文件
3. 设置 GRIB 2 的 paramId（=165）

转换仅在 paramId 定义了两个版本的情况下可行。

访问参数数据库网站查看转换是否可行。

练习

使用上一节生成的文件。

1. 将 p.grib1 转换为 GRIB 版本 2 文件 p_converted.grib2

    grib_set -s edition=2 p.grib1 p_converted.grib2
   

2. 查询 p.grib2 的参数命名空间。

    $ grib_ls -n parameter p.grib2
    p.grib2
    centre      paramId     shortName   units       name        
    ecmf        165         10u         m s**-1     10 metre U wind component 
    1 of 1 messages in p.grib2
   
    1 of 1 total messages in 1 file
   

3. 查询 p_converted.grib2 的参数命名空间。

    $ grib_ls -n parameter p_converted.grib2 
    p_converted.grib2
    centre      paramId     shortName   units       name        
    ecmf        165         10u         m s**-1     10 metre U wind component 
    1 of 1 messages in p_converted.grib2
   
    1 of 1 total messages in 1 files
   

   可以看到，不同版本的 GRIB 消息的变量信息相同。

4. 准备一个当前定义配置无法识别的 GRIB 消息：

    $ grib_ls -n parameter undefined.grib2 
    undefined.grib2
    centre      paramId     shortName   units       name        
    babj        0           unknown     unknown     unknown    
    1 of 1 messages in undefined.grib2
   
    1 of 1 total messages in 1 files
   

5. 将上面的文件转换成 GRIB 1 文件。

    $ grib_set -s edition=1 undefined.grib2 undefined.grib1
    ECCODES ERROR   :  concept: no match for paramId=0
    ECCODES ERROR   :  grib_set_values[0] edition (1) failed: Concept no match