- 积分
- 54
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-5-22
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
这是脚本,用来下载ERA-20C近一百多年的500hPa位势高度场数据
#!/usr/bin/env python
import calendar
from ecmwfapi import ECMWFDataServer
server = ECMWFDataServer()
yearStart = 1900 # 按你的需求来改
yearEnd = 2010 # 可改
months = [1,2,3,4,5,6,7,8,9,10,11,12] # 可改
years = range(yearStart, yearEnd+1)
requestDates=" "
for y in years:
for m in months:
requestDates = requestDates+str(y)+(str(m)).zfill(2)+'01/'
requestDates = requestDates[:-1]
print (requestDates)
server.retrieve({
"class": "e2",
"dataset": "era20c",
"date" :requestDates,#输出类似于19790101/19790201/19790301/19790401/19790501/19790601/19790701/197908等等。。。''
"expver" :"1",
"levtype": "pl",
"levelisit" :"500",
"param": "129.128",#注意不同的数据,这里不同,要改
"stream" :"moda",
"type" :"an",
"grid" :"2/2",#分辨率
"area" :"90/-180/-90/180",#area 5S-15N;10W-20E
"format" :"netcdf", #nc文件格式
"target" :"D:/ecmwf-api-client-python/500geo_new2_325.nc"#输出路径一般在该脚本所在目录
})
但是脚本运行到一半,就出现报错,想问下大家是什么问题?
Traceback (most recent call last):
File "D:\ecmwf-api-client-python\ecmwf4.py", line 34, in <module>
"target" :"D:/ecmwf-api-client-python/500geo_new2_325.nc"#输出路径一般在该脚本所在目录
File "D:\ecmwf-api-client-python\ecmwfapi\api.py", line 520, in retrieve
c.execute(req, target)
File "D:\ecmwf-api-client-python\ecmwfapi\api.py", line 488, in execute
assert size == result["size"]
AssertionError
|
|