python下载EC-Interim再分析资料

xunmengfei8

该程序用于批量下载EC-Interim再分析资料，每次重新启动程序后，具备自检跳过已存在数据，下载不存在数据的功能。

xunmengfei8

Masterpiece 发表于 2019-5-7 20:43
？？？？
这个玩意还要5个贡献？
官网上边能自动生成的脚本，这么做有点不够厚道吧？

第一次发帖子，还不会弄呢，望大家见谅。直接贴出来，大家看看，有用大家用用，不好勿喷。

import datetime
import os
import shutil
from ecmwfapi import ECMWFDataServer
server = ECMWFDataServer()
start_time = '19810101'
end_time = '20101231'
path = 'H:/ERA_Interim_daily0.5/'
pathout = 'D:/ERA-Interim-0.5/'

start = datetime.datetime.strptime(start_time,'%Y%m%d')
end = datetime.datetime.strptime(end_time,'%Y%m%d')
# print start
time = start
while time <= end:
    time = datetime.datetime.strftime(time,'%Y%m%d')
    dyear = time[:4]
    dmonth = time[4:6]
    dday = time[6:8]
    filepath = path + dyear +dmonth +dday+'.nc'
    filepathout = pathout + dyear + dmonth + dday + '.nc'
    print(filepath)
    if os.path.exists(filepath):
        print(os.path.getsize(filepath))
        if os.path.getsize(filepath)==399244228:
            time = datetime.datetime.strptime(time,'%Y%m%d')
        else:
            server.retrieve({
                'stream': 'oper',
                'levtype': 'pl',
                'param': '129.128/130.128/131.128/132.128/133.128/157.128',
                'dataset': 'interim',
                'step': '0',
                'grid': '0.50/0.50',
                'time': '00/06/12/18',
                'date': dyear +dmonth +dday,
                'type': 'an',
                'class': 'ei',
                'expver': '0001',
                'repres': 'SH',
                'levelist': '10/20/30/50/70/100/125/150/175/200/225/250/300/350/400/450/500/550/600/650/700/750/775/800/825/850/875/900/925/950/975/1000',
                'domain': 'G',
                'padding': '0',
                'format': 'netcdf',
                'target': filepath
                 })
            time = datetime.datetime.strptime(time,'%Y%m%d')
    else:
        server.retrieve({
            'stream': 'oper',
            'levtype': 'pl',
            'param': '129.128/130.128/131.128/132.128/133.128/157.128',
            'dataset': 'interim',
            'step': '0',
            'grid': '0.50/0.50',
            'time': '00/06/12/18',
            'date': dyear +dmonth +dday,
            'type': 'an',
            'class': 'ei',
            'expver': '0001',
            'repres': 'SH',
            'levelist': '10/20/30/50/70/100/125/150/175/200/225/250/300/350/400/450/500/550/600/650/700/750/775/800/825/850/875/900/925/950/975/1000',
            'domain': 'G',
            'padding': '0',
            'format': 'netcdf',
            'target': filepathout
             })
        shutil.move(filepathout,path)
        time = datetime.datetime.strptime(time,'%Y%m%d')
    time = time + datetime.timedelta(1)

Masterpiece

？？？？
这个玩意还要5个贡献？
官网上边能自动生成的脚本，这么做有点不够厚道吧？

1649518749

Masterpiece 发表于 2019-5-7 20:43
？？？？
这个玩意还要5个贡献？
官网上边能自动生成的脚本，这么做有点不够厚道吧？

气象家园早就有人发过贴，如何用py批下载EC数据。

xunmengfei8

Masterpiece 发表于 2019-5-7 20:43
？？？？
这个玩意还要5个贡献？
官网上边能自动生成的脚本，这么做有点不够厚道吧？

第一次发帖子，还不会弄

werewolf

不妨写个ERA5的吧，最新的高质量再分析数据。

小哥8号

期待ERA5的脚本

snoopy-_-99

嗯嗯，期待更新的脚本。

寻路人

不明觉厉，