爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 53940|回复: 49

[其他] ERA5再分析资料下载详解

  [复制链接]

新浪微博达人勋

发表于 2019-5-17 11:19:58 | 显示全部楼层 |阅读模式

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

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

x
本帖最后由 普鲁斯特 于 2019-5-17 11:19 编辑

2019年3月ECMWF上线了最新一代再分析资料——ERA5。ERA5计划覆盖1950至今的时间段,目前只有第一部分,也就是从1979年至今,另外一部分将于今年晚些时候上线。
ERA5采用了4D-Var变分同化,相对于ERA-Interim来说具有更高的时空分辨率,时间上可以提供每天逐小时,垂直上有137个混合层,插值到37个等压面上。水平分辨率默认为0.25°×0.25°


文末有我写的脚本,可以循环时间以及设定下载区域,并逐日保存。
1. era5_by_zns.py比较全,需要的单独下载这个文件即可
2. era5.py是官网例子
3. pdf文档也就是这个帖子的内容了,排版好一点而已,没有更多的内容,可以不用下载


下载准备:
  1. 所需工具:
   a) win10 Linux 子系统(本文示例所用,以下简称WSL)
   b) python (WSL中默认有2.7版本)
   c) CDS API client
  2. CDS账户:
   a) 注册CDS账户
     如果之前没有注册过CDS账户,点击以下链接注册:
       CDS API key
     注册完成后登陆CDS.
     然后会得到你自己的url以及Key Snipaste_2019-05-17_10-20-25.png

  b) 安装CDS API client
   CDS API全称是Climate DataStore Application Program Interface,也就是访问CDS的数据接口。
   在WSL的用户根目录下(/home/xxxx/)新建一个文件.cdsapirc 注意文件以点号开头),并将上面获得的url和key写入文件,保存退出。
  然后输入pip installcdsapi 安装CDS API client.

 如果发现上图问题,根据提示安装pip完成后再输入上面命令。

05162057

05162057





05162054

05162054


  3. 编写python脚本下载数据
    脚本的格式为
  1. <font face="微软雅黑">import cdsapi
  2. c = cdsapi.Client()

  3. c.retrieve("dataset-short-name",
  4.     {... sub-selection request ...},
  5.     "target-file")
  6. </font>
复制代码
如果不知道怎么写,可以在下面网址上选择变量及时间后点击Show APIrequest 按钮查询(类似于下载ERA-Interim)

05162251

05162251

05162253

05162253





将上面的语句保存为python脚本,运行即可。
Snipaste_2019-05-16_22-59-33.png
注:
  上图中的语法错误是由于脚本中包含了非ASCII字符。
  解决方法:
  1)参看下面链接
  2)简单的方法,把报错那行删去重写。
闲话:
  关于WSL很多人安装了Linux子系统后,就只是用系统给出的终端,颜色丑陋,性能低下,这里推荐MobaXterm。MobaXterm相比于ConEmu来说,最大优点在于集成了X11,NCL出图的时候可以很方便查看,不需要另外安装Xming客户端。
  打开MobaXterm:Session-->WSL
  这样每次打开后可以在主页上看到WSL了,如果没有选择左侧Session。另外,WSL也可以设置字体,主题等等。
  Have fun!


Snipaste_2019-05-16_23-06-49.png

Snipaste_2019-05-16_23-07-48.png


Snipaste_2019-05-16_23-10-40.png



Snipaste_2019-05-16_23-05-39.png

era5_by_zns.py

2.31 KB, 下载次数: 153, 下载积分: 金钱 -5

ERA5资料下载-Nanshui.pdf

589.23 KB, 下载次数: 142, 下载积分: 金钱 -5

era5.py

368 Bytes, 下载次数: 46, 下载积分: 金钱 -5

评分

参与人数 4金钱 +50 贡献 +2 收起 理由
skullcole + 5
huangbicheng + 10 很给力!
sreveryone + 20 很给力!
尽头的尽头 + 15 + 2

查看全部评分

本帖被以下淘专辑推荐:

  • · 资源|主题: 56, 订阅: 8
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2020-9-14 19:07:34 | 显示全部楼层
发现一个小错误,era5_by_zns.py中第21行程序 if(imo=='6'): 中imo是循环的次数不是对应的月份,应该改成 if(mo[imo]=='06')
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-5-20 08:21:39 | 显示全部楼层
嘿嘿嘿嘿 发表于 2019-5-19 20:29
感谢楼主分享!请问楼主知道ERA5的数据说明在哪里找吗?下载了降雪数据但是数据都是负值

https://confluence.ecmwf.int/display/CKB/ERA5+data+documentation
你看看这个~
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

新浪微博达人勋

发表于 2019-5-17 13:45:15 | 显示全部楼层
感谢分享~~~
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2019-5-19 20:29:51 | 显示全部楼层
感谢楼主分享!请问楼主知道ERA5的数据说明在哪里找吗?下载了降雪数据但是数据都是负值
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-5-21 17:17:32 | 显示全部楼层
厉害了!!谢谢分享!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-5-22 14:54:00 | 显示全部楼层
楼主好,请问一下。我用web api下过ERA-Interim。我使用类似的手法创建完.cdsapirc文件,用系统自带的命令行窗口运行py文件的时候,报错了。
WechatIMG43899.png WechatIMG43898.png
可否贴一下您.cdsapirc的内容(只有两行?)api.py文件需要修改吗?

WechatIMG43899.png
QQ图片20170314195121.png
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-5-23 08:52:28 | 显示全部楼层
不知道楼主这样批量下载的速度怎么样?我之前写过一个python批量下载程序,速度太慢了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-5-23 09:02:21 | 显示全部楼层
写的很详细,谢谢分享
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-5-23 10:34:52 | 显示全部楼层
楼主的程序单单是6-8月的吧?下面涉及到循环里的dayEnd函数可以再优化一下。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-5-24 08:59:33 | 显示全部楼层
勿忘草 发表于 2019-5-23 08:52
不知道楼主这样批量下载的速度怎么样?我之前写过一个python批量下载程序,速度太慢了

这个是和网速有关的吧,反正我下的时候也不快
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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