使用欧洲数值天气预报EC run WRF的一点经验探索（二）

andrewsoong

前段时间是在是太忙，累成狗！真没时间做这些了。
但有一些坛友咨询这方面的问题，尤其是山东省局的朱同学，故而昨晚弄到很晚，现在写一下，欢迎大家来喷，无妨，但要提出建议，否则我可不待见。哈哈哈，言归正传：


再上一个帖子里面http://bbs.06climate.com/forum.p ... ad&tid=35277&extra=写到了将EC数据转换成grib1数据，接下来就是运行WPS了。我是用的是WRF3.7-bugfix，将处理好的EC数据link.csh到WPS，Vtable使用ECMWF的，即可UNGRIB。需要指出的是应朱同学要求，这里使用T639资料作为土壤温湿度的来源。EC也缺少地面气压，可用T639的补上。


T639资料的使用：WPS/ungrib/src
首先将rd_grib2.F替换，否则读不出T639的土壤温湿度数据，替换程序WPS/ungrib/src rd_grib2.F (39.3 KB, 下载次数: 93)

2018-4-11 09:51 上传

点击文件名下载附件


然后重新编译WPS即可！Vtabele使用这个： vtable.T639-soil (1.46 KB, 下载次数: 88)

2015-9-23 10:25 上传

点击文件名下载附件 然后就可以ungrib T639的资料。

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

然后将ungrib后的T639拷贝到EC的WPS下，即可运行met，得到metXXXX.nc文件。
接下来连接到WRFV3/test/em_real下，修改namelist： num_metgrid_levels                  = 18,
num_metgrid_soil_levels             = 4,
即可运行real。

写的比较粗，单流程基本就这样，本人已经跑通。结果没时间看了。以后再说。欢迎各位给出意见和测评结果。QQ讨论群：气象家园-WRF模式应用 546530151
20150923于呼和浩特。

pandasp1213

感谢楼主分享的经验，我最近在用18年的ec数据做，grib_api用的1.24.0，到ungrib时出现了一些问题，
End-of-record mark (7777) not found
stopping in findgrib in gribcode
I could not find the GRIB string in the input file after testing the first 100,000 bytes.


这个文件是2018年10月10日15时的，文件处理成grib1后大小224M，另外试了一个12时的文件，大小142M，12时的文件ungrib没有问题。
12时的文件相当于再分析场，15时是3小时的预报场。

不同时次的预报场，大小也不一致。
不知楼主碰到过类似问题没有，感谢

tornadoarps

我用楼主的方法将混合格式转换为grib1格式会出错，/usr/local/grib_api/bin/grib_set -f -s editionNumber=1 W_NAFP_C_ECMF_20170110055402_P_C1D01100000011006001 W_NAFP_C_ECMF_20170110055402_P_C1D01100000011006001.grib1
出现以下错误：
GRIB_API ERROR   :  unable to set NV=276 as long (Encoding invalid)
GRIB_API ERROR   :  grib_set_values[0] editionNumber (1) failed: Encoding invalid
GRIB_API ERROR   :  Key "numberOfVerticalCoordinateValues": Trying to encode value of 276 but the maximum allowable value is 255 (number of bits=8)
改用以下命令则运行正常：/usr/local/grib_api/bin/grib_copy -v -w edition=1
希望大家注意！

qq469015280

请问我应经下载好了fnl数据，是grib1格式
WPS读取时，出现错误。
它为什么无法识别我的数据呢？

xo_shi

好帖子，感谢楼主的分享，按照楼主思路，经过这段时间断断续续的折腾，终于弄出来了。感觉最重要的是grib-api的安装使用。我的grib-api反反复复装了3天，仍然不能实现转码（grib_set），幸好grib_get_data可以使用，所以就参照以前用wgrib解码给模式应用的经验，类似来用grib_get_data。
流程是：通过grib_get_data来输出各要素数据，然后在脚本里面用sed删除说明行（Lat，LON，value），用awk输出第三列（value），这样就得到基本的变量，接着用fortran来读写，并转成中间格式。
需要说明的是：（1）细网格有17层（我们这边cmacast只解了10层（最开始我尝试的是10层的数据），可能有些单位在cmacast使用上有改进，层次要更全）
                        （2）土壤变量和地面气压是没有的，需要向楼主这样借用T639（对于基层业务单位，T639资料相对算是好的，当然如果有互联网用GFS来补充缺省量更好）
                        （3）细网格的地面层变量0.125的分辨率，而且是全球的，用grib_get_data导出来数据是很大的，因此要用sed去掉0.125/0.375/0.625/0.875的行，再用awk输出第三列，这样处理速度要快一些。
                        （4）最后回到grib-api的安装使用，我觉得和机器环境有很大关系，jasper和grib-api的安装要注意enable-shared和-fpic这两个编译选项，有些同行反映jasper安装好了但grib-api不认，应该是grib-api解析不了jasper的静态库进而出现“bad value”

西西

扬子一出关就有大手笔啊！

andrewsoong

西西 发表于 2015-9-23 10:37
扬子一出关就有大手笔啊！

姐姐多多指点~~~

平凡的幸福

真心感谢@andrewsoong，不耐其烦的一步步指导，这个过程遇到很多问题，占用了@andrewsoong大量的时间。当然借此感谢气象家园这个平台，让我遇到很多优秀的同学，从中学到很多。

andrewsoong

平凡的幸福 发表于 2015-9-23 16:03
真心感谢@andrewsoong，不耐其烦的一步步指导，这个过程遇到很多问题，占用了@andrewsoong大量的时间。当然 ...

andrewsoong

xo_shi 发表于 2015-11-14 22:00
好帖子，感谢楼主的分享，按照楼主思路，经过这段时间断断续续的折腾，终于弄出来了。感觉最重要的是grib-a ...

超赞

461449061

谢谢楼主分享，非常感谢

平凡的幸福

xo_shi 发表于 2015-11-14 22:00
好帖子，感谢楼主的分享，按照楼主思路，经过这段时间断断续续的折腾，终于弄出来了。感觉最重要的是grib-a ...

土壤温度和湿度没有，地面气压也没有吗/

xo_shi

平凡的幸福 发表于 2015-11-20 16:20
土壤温度和湿度没有，地面气压也没有吗/

是的，没有地面气压