#!/usr/bin/env python
import calendar
from ecmwfapi import ECMWFDataServer
server = ECMWFDataServer()
hindcastYearStart = 1994
hindcastYearEnd = 2014
hindcastMonthStart = 1
hindcastMonthEnd = 12
def retrieve_cma_reforecast():
for hindcastYear in list(range(hindcastYearStart, hindcastYearEnd + 1)):
for hindcastMonth in list(range(hindcastMonthStart, hindcastMonthEnd + 1)):
numberOfDays = calendar.monthrange(hindcastYear, hindcastMonth)[1]
for hindcastDay in list(range(numberOfDays)):
hindcastDate = '%04d%02d%02d' % (
hindcastYear, hindcastMonth, hindcastDay + 1)
cma_reforecast_request(hindcastDate)
def cma_reforecast_request(hindcastDate):
modelVersionDate = "2014-05-01"
target = "data_s2s_%s.grb" % (hindcastDate)
server.retrieve({
"class": "s2",
"dataset": "s2s",
"date": modelVersionDate,
"expver": "prod",
"hdate": hindcastDate,
"levtype": "sfc",
"origin": "babj",
"param": "165/166",
"step": "24/to/1440/by/24",
"stream": "enfh",
"target": target,
"time": "00",
"number": "1/2/3",
"type": "pf",
})
if __name__ == '__main__':
retrieve_cma_reforecast() |