利用python下载欧洲气象中心的数据（ECMWF）2018.3.31更

月夜猫妖

zuichu0128 发表于 2018-4-19 15:23
谢谢楼主，我按照这个里面的步骤来，第一步有点问题，要把自己账号的url和key以及email进行修改，这个步 ...

其实改改api中的代码也可以，修改aip.py内def _get_apikey_from_environ():下头的key、url、email就好了，把他们修改成自己的就ok了
key = "xxx"
url = "https://api.ecmwf.int/v1"
email = "xxx"

小莹子

你好，我想请教一下下载范围的问题。就是我下载的范围是50/60/15/110，然后出现错误WARNING: Field PSFC has missing values at level 200100 at (i,j)=(1,63)
ERROR: Missing values encountered in interpolated fields. Stopping.
查看之后发现是驱动范围和实际范围不符。在Google上查看是说“以为输入数据参考位置是它的中心，但实际上是在西南角”这个怎么理解哪？谢谢，期待您的回复

coolhoo

学习学习，谢谢分享

ShuyingChenUCAS

小莹子 发表于 2018-4-23 20:57
你好，我想请教一下下载范围的问题。就是我下载的范围是50/60/15/110，然后出现错误WARNING: Field PSFC ha ...

哈哈哈，发现你在两个帖子里都找到了我，回复过你啦：
我是这样做的，如果实际范围是50/60/20/110，下载ERA-Interim的时候选择东西南北各大1度的范围下载，也就是51/59/19/111，你试试行不行。还是有不行的话可能不是下载数据范围的问题。

小莹子

ShuyingChenUCAS 发表于 2018-4-24 19:11
哈哈哈，发现你在两个帖子里都找到了我，回复过你啦：
我是这样做的，如果实际范围是50/60/20/110，下载 ...

啊，不好意思，逛的帖子多了，就也不知道自己干了什么，谢谢您的回复。

ShuyingChenUCAS

小莹子 发表于 2018-4-25 09:27
啊，不好意思，逛的帖子多了，就也不知道自己干了什么，谢谢您的回复。

没事没事，完整的问答也许能帮到别人，我把你的回答粘过来了：
你好，我试了20度不行，30度就可以了，谢谢。但是现在我不小心把metgrid给删了，然后哎，剩下的。。

答：既然30度可以那么你的实际研究范围可能不是50/60/20/110。总之在下载数据前先仔细确定自己的研究范围，下载数据时比这个确定的范围大一些，大多少自己把握吧，当然越大越不容易出错但没必要，我每次大约是四边都向外扩展一度的样子。
耐心一下，多试几次总会转出来的，加油。

mockings

感谢楼主分享，，，，今天第一次尝试用这种方法，遇到了这种错误，又说文件超过30G，又说版本不符。。。现在一脸蒙蔽。。。请问这是因为python自身的问题吗

1. 
   
2.   File "C:\Users\mock1\Anaconda2\lib\site-packages\ecmwfapi\api.py", line 340, in call
   
3.     raise APIException("ecmwf.API error 1: %s" % (self.last["error"],))

复制代码

这是我运行的代码

1. #!/usr/bin/env python
   
2. from ecmwfapi import ECMWFDataServer
   
3.    
   
4. server = ECMWFDataServer()
   
5.    
   
6. server.retrieve({
   
7.     'stream'    : "oper",
   
8.     'levtype'   : "sfc",
   
9.     'param'     : "186.128/188.128",
   
10.     'dataset'   : "interim",
    
11.     'step'      : "0",
    
12.     'grid'      : "0.5/0.5",
    
13.     'time'      : "0600/1800",
    
14.     'date'      : "1979-01-031/to/2014-12-31",
    
15.     'type'      : "an",
    
16.     'class'     : "ei",
    
17.     'format'    : "netcdf",
    
18.     'area'      : "40/85/25/105",
    
19.     'domain'    : "G",
    
20.     'resol'     : "auto",
    
21.     'target'    : "interim_1979-01-01to2014-12-31_00061218.nc"
    
22. })

复制代码

谢谢楼主了

ShuyingChenUCAS

mockings 发表于 2018-4-29 19:35
感谢楼主分享，，，，今天第一次尝试用这种方法，遇到了这种错误，又说文件超过30G，又说版本不符。。。现 ...

看你代码应该是没什么问题，重新提交一次试试，你提取的时间有点长可能真是超过30G限制了。
版本是Python 2.7.x 或 Python 3 的python，如果这个没问题的话可能是平台的原因么？
不如，你可以先试试下一年的数据？

AcesLin

楼主，想请教一个问题，我下载了北半球的某一年的风速数据，然后文件格式是grid的，在arcgis打开为什么是一片黑的？？请问是什么问题？
谢谢解答~~

ShuyingChenUCAS

AcesLin 发表于 2018-5-3 20:06
楼主，想请教一个问题，我下载了北半球的某一年的风速数据，然后文件格式是grid的，在arcgis打开为什么是一 ...

grib数据？我用arcmap打开我下载的u和v分量没有问题，是不是你打开方式有问题或者数据没有下对？