GR_
10/21
本文最后更新于
2023年10月21日，已超过
569天没有更新。如果文章内容或图片资源失效，请留言反馈，我会及时处理，谢谢！
问题
有一个grib2数据，使用wgrib2查看元信息：
1 | $ wgrib2 in.grb |
其中SPFH场垂直层信息为reserved，未能显示正确的等压面层次。
使用eccodes中grib_ls工具查看
1 | $ grib_ls in.grb |
输出如下：
其中typeOfLevel和level两个key信息对于比湿q（即wgrib2中显示的SPFH）有问题。
层次类型值为unknown，而不是像正常场一样为isobaricInhPa，其level单位是Pa，而其它正常显示的场单位是hPa，应该就是这个问题导致在wgrib2查看时出现异常。
修改key信息
通过eccodes提供的grib_set工具修改key信息使grib文件正常显示元信息。
1 | $ grib_set -s typeOfLevel=isobaricInhPa -w typeOfLevel=unknown in.grib out.grib |
-w只对匹配该选项条件的GRIB messages执行set操作，其他不符合的messages原封不动拷贝到新的grib文件。-s设置key=value
修改后wgrib2查看
附:eccodes简要安装说明
最新的eccodes基于CMake构建。
1 2 3 4 5 6 7 8 9 10 11 | $ tar xzvf eccodes-2.30.0-Source.tar.gz$ cd eccodes-2.30.0-Source $ cmake -DCMAKE_INSTALL_PREFIX=/path/to/install/eccodes -DENABLE_AEC=OFF -B build$ cmake --build build --target install# 设置环境变量$ vi ~/.bashrcexport PATH=/path/to/install/eccodes/bin:$PATHexport ECCODES_DIR=/path/to/install/eccodes$ source ~/.bashrc |
eccodes一些功能需要依赖其它库，比如NETCDF/JPG/PNG等，具体cmake开关见官网说明。
版权属于：
MetMan's Blog
本文链接：
https://blog.metman.top/index.php/archives/31/（转载时请注明本文出处及文章链接）
作品采用：
《署名-非商业性使用-相同方式共享 4.0 国际 》许可协议授权
