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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 194849|回复: 361

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

  [复制链接]

新浪微博达人勋

发表于 2014-12-2 20:01:21 | 显示全部楼层 |阅读模式

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

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

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, 下载次数: 1198)

评分

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

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

0
早起挑战累计收入
发表于 2014-12-2 22:38:40 | 显示全部楼层
看了下楼主之前问过这个问题,看来是解决了来分享给大家,很值得支持!
密码修改失败请联系微信:mofangbao
回复 支持 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资料,刚运行一会就会有这个报错…是什么原因呢?求破!
密码修改失败请联系微信:mofangbao
回复 支持 2 反对 0

使用道具 举报

新浪微博达人勋

发表于 2017-7-26 15:18:50 | 显示全部楼层
If any masked field is ingested, then a LANDSEA field is recommended

If the WRF model is going to be run with the Noah LSM model, then at least 2 levels of Soil Temperature and Soil Moisture are required.

SST is needed if you plan on updating SST during the model run (often done for very long simulations).

Soil Height is recommended but not required. If Soil Height is supplied, then SOIL Temperatures and TSK can be adjusted to the WRF model terrain height.

Water equivalent snow depth (SNOW) is a nice field to have, but not required.

SEAICE is good to have for a high latitude winter case, but it is not required.
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 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.
求解答
密码修改失败请联系微信:mofangbao
回复 支持 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
     ...

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

使用道具 举报

新浪微博达人勋

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

使用道具 举报

新浪微博达人勋

发表于 2014-12-3 10:01:19 | 显示全部楼层
,感谢楼主,100个
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2014-12-3 12:44:42 | 显示全部楼层
楼主,这个绝对不错啊
密码修改失败请联系微信:mofangbao
回复 支持 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
密码修改失败请联系微信:mofangbao
回复 支持 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 就可以了。
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2014-12-2 22:24:15 | 显示全部楼层
感谢楼主的分享。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-12-2 22:25:32 | 显示全部楼层
能否麻烦楼主简单介绍一下如何运行*.py脚本文件?谢谢了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-12-3 08:44:43 | 显示全部楼层
32个赞~
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2014-12-3 12:28:45 | 显示全部楼层
感谢分享~~~
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

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

本版积分规则

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

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

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