请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

新浪微博登陆

只需一步, 快速开始

QQ登录

只需一步,快速开始

搜索
查看: 46769|回复: 343

使用ECMWF资料作为WRF的初始场,以ERA-Interim资料为例

  [复制链接] |关注本帖

新浪微博达人勋

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 夏朗的芒果 于 2015-6-27 07:49 编辑

通过论坛各位的帮助,我总结了一些关于如何使用ECMWF资料运行WRF的经验,肯定有不足之处,还望各位批评指正。
首先需要获得正确的数据

论坛里很多帖子讲如何下载ERA-Interim数据的步骤:
一般情况下,需要分别下载地表变量和高空变量。高空变量又可以分为气压层和模式层。由于ECMWF数据分辨率较高,如果下载全部变量的话,文件一般较大,增加下载时间。如果只下载WRF运行所必需的变量,将大大缩减下载时间。
根据WRF官网教程翻译,http://bbs.06climate.com/forum.php?mod=viewthread&tid=30940,WRF需要一些必要的变量才能运行起来。我根据Table 128加入了相应ECMWF所对应的Grib code。

以下部分摘自翻译。
-----------------------------------------------------------------------------------------------
三维数据(例如,气压层的数据)
· 温度(130.128)
· 风场u,v分量(131.128, 132.128)
· 位势高度(129.128对应位势)
· 相对湿度(157.128/比湿 (133.128)
二维数据
· 地面气压 (134.128)
· 平均海平面气压(151.128)
· 地表温度(235.128)
· 2米温度 (167.128, 2m 露点温度168.128)
· 2米相对湿度或比湿
· 10风场u,v分量(165.128, 166.128)  
       如果需要用到masked,建议提供LANDSEA场。  
      如果在进行物理过程参数化时需要用到Noah LSM陆面模式,需要提供两层Soil Temperature (139,170,183,236.128)和Soil       Moisture(39,40,41,42.128)。
      如果模式运行时需要更新SST,需要提供SST(34.128)。
      建议提供Soil Height(可选),这样就可以调整土壤温度达到与WRF模式地形高度向适应。
      建议提供SNOW(141.128场。
      建议对高纬度冬季个例,提供SEAICE (31.128) 场。
-------------------------------------------------------------------------------------------------

根据WRF输入变量的要求,我改写了批量下载气压层pl和地面层sfc的python脚本,大家下载后只需修改相应的时间即可。
lupit_pl_14_20.py (1.16 KB, 下载次数: 839)

评分

参与人数 6金钱 +49 贡献 +5 收起 理由
XinxuZhao + 2
花纹云彩 + 2 很给力!
balfulosa + 26 谢谢楼主!
mofangbao + 15 + 5
zhxjanet + 2 很给力!
林子木 + 2 很给力!

查看全部评分

已有16人关注本帖

liusinuoWillnlk467回到未来BLUEANGEL2012
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2014-12-2 22:38:40 | 显示全部楼层 |取消关注该作者的回复
看了下楼主之前问过这个问题,看来是解决了来分享给大家,很值得支持!
密码修改失败请联系qq:937062711
回复 支持 2 反对 0

使用道具 举报

新浪微博达人勋

发表于 2015-1-22 16:21:27 | 显示全部楼层 |取消关注该作者的回复
Traceback (most recent call last):
  File "D:\tools\Python27\lupit_pl_14_20.py", line 33, in <module>
    'target'  : "ERA-Int_pl_19890531-19891001.grib"
  File "D:\tools\Python27\ecmwfapi\api.py", line 414, in retrieve
    c.execute(req, target)
  File "D:\tools\Python27\ecmwfapi\api.py", line 377, in execute
    size = self._transfer(result["href"], target, result["size"])
  File "D:\tools\Python27\ecmwfapi\api.py", line 76, in wrapped
    return func(*args,**kwargs)
  File "D:\tools\Python27\ecmwfapi\api.py", line 340, in _transfer
    assert total == size
AssertionError
我用版主的脚本下载interim资料,刚运行一会就会有这个报错…是什么原因呢?求破!
密码修改失败请联系qq:937062711
回复 支持 2 反对 0

使用道具 举报

新浪微博达人勋

发表于 2017-10-27 11:27:18 | 显示全部楼层 |取消关注该作者的回复
楼住,数据下载下来是grib格式,但是你的教程里 链接的时候链接的是grb格式,我按照教程,会提示 no match,改成grib后缀,./ungrib.exe会提示错误
Start_date =  1979-07-22_00:00:00 ,      End_date = 1979-07-24_18:00:00
output format is WPS
Path to intermediate files is ./
ERROR: Grib2 file or date problem, stopping in edition_num.
求解答
密码修改失败请联系qq:937062711
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2015-8-28 10:57:20 | 显示全部楼层 |取消关注该作者的回复
微微一笑 发表于 2015-1-22 16:21
Traceback (most recent call last):
  File "D:\tools\Python27\lupit_pl_14_20.py", line 33, in
     ...

时隔半年,我也出现了这个问题。解决方法是运行脚本成功下载时,会提供一个下载地址,利用这个地址用迅雷下载,可以避免这个问题。
密码修改失败请联系qq:937062711
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2014-12-3 09:53:02 | 显示全部楼层 |取消关注该作者的回复
我也一直有这个疑问,楼主好人啊!!!太感谢了~
密码修改失败请联系qq:937062711
回复 支持 0 反对 1

使用道具 举报

新浪微博达人勋

发表于 2014-12-3 10:01:19 | 显示全部楼层 |取消关注该作者的回复
,感谢楼主,100个
密码修改失败请联系qq:937062711
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2014-12-3 12:44:42 | 显示全部楼层 |取消关注该作者的回复
楼主,这个绝对不错啊{:soso_e179:}
密码修改失败请联系qq:937062711
回复 支持 0 反对 1

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2014-12-3 16:24:24 | 显示全部楼层 |取消关注该作者的回复
hxyj 发表于 2014-12-2 22:25
能否麻烦楼主简单介绍一下如何运行*.py脚本文件?谢谢了

这个是python的脚本文件,运行详情见http://bbs.06climate.com/forum.php?mod=viewthread&tid=24435
密码修改失败请联系qq:937062711
回复 支持 0 反对 1

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2014-12-20 08:53:02 | 显示全部楼层 |取消关注该作者的回复
陈颖 发表于 2014-12-17 21:31
非常感谢分享!很有帮助!
请问楼主 那两个时间无关的文件  Geopotential.grib 和 Land-sea mask.grib    ...

这两个是可选的,geopotential.grib 别被它名字误导了,其实就是地形高度,这个在进行插值土壤温度有一点用。land-sea mask 就不用我说了吧?
这两者下载下来后和之前的没多大区别。
先修改namelist.input里面的时间项,起止时间都是1989-01-01_12
前者的prefix = z ,后者的prefix =lsm
然后选择 vtable.ecmwfvtable。对两个文件分别进行 ungrib.exe
这样就有 sfc、pl、z、lsm 这四类解压后的文件了。
在进行metgrid 前,修改
constants_name = 'LSM:1989-01-01_12', 'Z:1989-01-01_12'
fg_name = 'sfc',‘pl’
最后进行 metgrid.exe 就可以了。
密码修改失败请联系qq:937062711
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2015-6-29 21:53:54 | 显示全部楼层 |取消关注该作者的回复
楼主  运行WRF时候 下载这个初始场数据必须是全球的么?可以把范围缩小到 比最大嵌套稍大一点的区域吗?
密码修改失败请联系qq:937062711
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2014-12-2 22:24:15 | 显示全部楼层 |取消关注该作者的回复
感谢楼主的分享。
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2014-12-2 22:25:32 | 显示全部楼层 |取消关注该作者的回复
能否麻烦楼主简单介绍一下如何运行*.py脚本文件?谢谢了
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2014-12-3 08:44:43 | 显示全部楼层 |取消关注该作者的回复
32个赞~
密码修改失败请联系qq:937062711
回复

使用道具 举报

新浪微博达人勋

发表于 2014-12-3 12:28:45 | 显示全部楼层 |取消关注该作者的回复
感谢分享~~~
密码修改失败请联系qq:937062711
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表