[经验]简单说一下如何用欧洲中心的资料驱动WRF

sunshinezhpb

最近看到有人在问，我就简单说一下我的用欧洲中心的资料驱动WRF的一些经验，算是抛砖引玉，有什么不对的地方还请大家指出。

首先是欧洲中心资料的下载，现在一般用的是ERA-Interim，下载地址为http://apps.ecmwf.int/datasets/data/interim_full_daily/ ，具体的下载步骤论坛上有相关帖子，非常详细，在此就不多赘述了。
除了直接在官网上下载，也可以用脚本下载。论坛上也有如何利用python脚本进行下载的帖子，也非常详细。在这个网站上https://software.ecmwf.int/wiki/ ... ta+servers+in+batch  也包括了Perl和Java脚本如何下载。

在下载下来的数据里，是包括高空（ERA-Int_pl_*）和地面（ERA-Int_sfc_*）两部分的。因此，我们在ungrib解码时也要解码两次。以1992年5月份的资料为例——ERA-Int_pl_19920501_19920531.grb和ERA-Int_sfc_19920501_19920531.grb：
链接Vtable:                       ln -sf ungrib/Variable_Tables/Vtable.ECMWF Vtable
链接高空数据：                ./link_grib.csh ./DATA/ERA-Interim/ERA-Int_pl_19920501_19920531.grb
编辑namelist.wps:             prefix = '3D',
运行ungrib:                      ./ungrib.exe
链接地面数据：               ./link_grib.csh ./DATA/ERA-Interim/ERA-Int_sfc_19920501_19920531.grb
编辑namelist.wps:            prefix = 'SFC',
运行ungrib:                     ./ungrib.exe
编辑namelist.wps里metgrid部分： fg_name = '3D','SFC'   并运行./metgrid.exe，后续部分跟用其余资料一致。

实际上，在使用ERA-Interim资料时，wrf自带的Vtable不止有Vtable.ECMWF，也有 Vtable.ERA-interim.pl 和 Vtable.ERA-interim.ml，其中Vtable.ERA-interim.pl对应的是气压层的ERA-Interim资料， Vtable.ERA-interim.ml对应的是模式层的ERA-Interim资料。
想更深入的了解可以看看这个网站：http://www.geomar.de/mitarbeiter/fb1/me/jharlass/wrf/erai/

Maggie89

其实如果Vtable 有把高空和地面数据都完整地列入，只需要执行一次ungrib.exe。不过前提是执行./link_csh.exe 时，需把所有数据都连接到执行文件中。

askira2011

谢谢楼主分享啊 ！刚好有急用

sun_shine_Xia

大赞一个！！

北风吹

谢谢分享，支持以下／／／

红豆夹薏仁

楼主，你好，我的步骤全部跟你一样，但是metgrid的时候出现了这样的问题：
WARNING: Field PRES has missing values at level 200100 at (i,j)=(1,53)
WARNING: Field PMSL has missing values at level 200100 at (i,j)=(1,53)
WARNING: Field PSFC has missing values at level 200100 at (i,j)=(1,53)
ERROR: Missing values encountered in interpolated fields. Stopping.
application called MPI_Abort(MPI_COMM_WORLD, 0) - process 0
应该是地面数据出了问题，但是又找不到原因，请楼主帮帮忙

sunshinezhpb

红豆夹薏仁 发表于 2014-11-6 20:18
楼主，你好，我的步骤全部跟你一样，但是metgrid的时候出现了这样的问题：
WARNING: Field PRES has missi ...

说的是你数据缺省...请仔细检查数据。
确定你是否在两次ungrib.exe时修改了namelist里的prefix 项并重新链接了数据。
还有，确定你下的地面数据和你在namelist里设置的所选区域是一致的，并确保你下载的数据完整（如果你是在欧洲中心网站直接下载数据有时候会断掉但提示下载完成）。

红豆夹薏仁

sunshinezhpb 发表于 2014-11-6 20:47
说的是你数据缺省...请仔细检查数据。
确定你是否在两次ungrib.exe时修改了namelist里的prefix 项并重新 ...

谢谢，高空和地面下的范围是一样的，我分开测试过，高空没问题，地面数据出了问题，可是我重新在欧洲中心下载资料后还是同样的问题，我确定每次都链接了新数据，修改了namelist，所以还是找不到原因

Maggie89

你可以尝试用WPS/util 里的 g1print.exe 来确保在欧洲中心下载的资料是否完整。

沿江西xi

欧洲中心的资料不是要钱的么？