使用欧洲数值天气预报EC run WRF的一点经验探索（三）关于2016下半年国家局修改数...

andrewsoong

在前面两个帖子关于使用国家局下发的EC预报场驱动WRF进行预报，有不少业务人员反馈。也得到了很多帮助。
从2016年8月31日开始，为进一步支持现代气象业务发展，根据ECMWF业务调整要求（停止使用复杂压缩格式，改为bzip2压缩文件），2016年8月31日12时次（开始接收和下发时间为北京时间9月1日2:40）起，调整目前通过CMACast广播下发、以及通过中国气象局骨干网、全国宽带网（MPLS VPN）共享的ECMWF高分辨率数值预报产品。那么以前的两个帖子中的部分流程就无法使用了，一直想再写个帖子，苦于没时间（主要是瞎忙）。现在，在原帖的基础上改进吧，先奉上前面的两个帖子地址以便于对比哈：http://bbs.06climate.com/forum.p ... 77&page=1#pid471423
和http://bbs.06climate.com/forum.php?mod=viewthread&tid=38302

主要说一下需要改进的地方：其他参考前面两个帖子：
1，数据的解压：2016年9月1日开始，根据ECMWF业务调整要求，全部ECMWF高分辨率数值预报产品，停止使用复杂（二阶）压缩格式，恢复为原简单压缩（grib_simple）格式，改为bzip2格式压缩文件（需自行解压缩），文件扩展名改为.bz2。注意：此次更新，增加了地面气压！共19层。

Linux/AIX/UNIX等操作系统上解压缩说明：
Linux/AIX/UNIX等操作系统提供bzip2工具包。
解压缩参考命令：

游客，如果您要查看本帖隐藏内容请回复

命令执行完成后，该bz2文件自动删除，当前目录下出现解压缩后的文件，文件名为C1D08040000080400001。
注：
1）解压缩之前，请确认该目录（文件夹）有写权限。
2）若执行命令时出现“command not found（命令未找到）”等错误信息，请联系该服务器管理员安装bzip2工具包。转换后的C1D08040000080400001为GRIB1和GRIB2混合的，其中grib2格式的数据为预报降水相态等。


2，提取WPS所需数据命令如下：

游客，如果您要查看本帖隐藏内容请回复

PS：不要喷，不要说隐藏回复可见浪费你的时间什么的，为了杜绝伸手党吧，谁都很忙的。你如果觉得回复可见浪费了你的时间，那你大可不必看了，自己研究去吧。我相信你会更多的时间to waste。。。。。。解释权归清风@mofangbao
如此得到的20170109-1200new.grib1就能被WPS ungrib了，剩下的工作参考帖子（二）。
欢迎讨论，QQ群：气象家园-WRF模式应用546530151。

2017-01-11于内蒙古气象局-呼和浩特。

andrewsoong

又是那隻貓 发表于 2017-1-11 17:14
兹泚宋首席

谢谢猫

cbcbc

W_NAFP_C_ECMF_20180424054725_P_C1D04240000042406001的bz2数据，
用grib_copy -v -w edition=1 filename1 filename2指令提取的grib1要素成功，
但是用WPS下的ungrib.exe数据时，报下面的错

End-of-record mark (7777) not found 806629396
Sec0(1) =  8462100 0
164  34 SSTK     1    0    0  2018-04-24_00:00 + 03
*** stopping in findgrib in gribcode ***
        I could not find the GRIB string in the input file
        after testing the first 100,000 bytes.
        The file may be corrupt or it is not a GRIB file.
        Perhaps a gzipped GRIB file or netcdf?
STOP findgrib
什么意思？是说在806629396字节处应该是7777？
我打开看该文件末尾是的啊！。。。。。。。。
怎么回事呢？有人知道吗？请指教！谢谢！

平凡的幸福

zqvvcl 发表于 2017-7-26 10:22
关于数据提取，我用grib_set -f -s editionNumber=1 会出错，用grib_copy -v -w edition=1没问题
但是到了 ...

你好，我也是解码第一个数据正常，第二个数据就出问题了
请问你是怎么解决的

zqvvcl

关于数据提取，我用grib_set -f -s editionNumber=1 会出错，用grib_copy -v -w edition=1没问题
但是到了ungrib这一步，每次预报场的第一个数据（是分析场吗？）可以顺利通过，后面预报时次的数据就出错了：
===========================
*** Starting program ungrib.exe ***
Start_date =  2017-07-27_03:00:00 ,      End_date = 2017-07-27_03:00:00
output format is WPS
Path to intermediate files is ./
End-of-record mark (7777) not found1863675669
Sec0(1) =  8462100 0
*** stopping in findgrib in gribcode ***\n
\tI could not find the GRIB string in the input file
\tafter testing the first 100,000 bytes.
\tThe file may be corrupt or it is not a GRIB file.
\tPerhaps a gzipped GRIB file or netcdf?\n
findgrib
=================================
我用./util/g1print.exe 072700copy.grib查看好像是正常的。大家遇到过这个问题吗？求教

又是那隻貓

兹泚宋首席

平凡的幸福

先顶！d=====(￣▽￣*)上去

topmad

你就是砖家
呲呲soong首席

这么牛逼的东东

andrewsoong

topmad 发表于 2017-1-11 17:41
你就是砖家
呲呲soong首席

我是砖家，老5说的没错。我就是搬砖的砖家

tornadoarps

回复一下

andrewsoong

tornadoarps 发表于 2017-1-11 17:53
回复一下

你给的那个我测试了，ungrib的时候出问题。你是怎么解决的

tornadoarps

andrewsoong 发表于 2017-1-11 18:04
你给的那个我测试了，ungrib的时候出问题。你是怎么解决的

没有解决。我用grib_set命令时出错，用grib_copy正常，生成的grib1文件用wgrib解码也正常，但ungrib处理时出错。真是奇怪了，难道是grib api版本过低？我的grib_api 版本是1.12.3

andrewsoong

tornadoarps 发表于 2017-1-11 18:19
没有解决。我用grib_set命令时出错，用grib_copy正常，生成的grib1文件用wgrib解码也正常，但ungrib处理 ...

这个问题我测试过很多版本，只有1.10.4可以。其他版本都出错。这个我测试过，ungrib都没问题。