Versions Compared

Old Version 100

changes.mady.by.user Matthew Manoussakis

Saved on

compared with

New Version 101

changes.mady.by.user Matthew Manoussakis

Saved on

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Info
  • The objective of this example is to demonstrate how to make a request efficient by iterating properly over several hindcastYears, hindcastMonths and hindcastDays
  • It can be used as a starting point however you need to adapt it to your needseg to set the keyword values (eg hindcastYear).
  • In this way you can extend this request to download the whole S2S CMA reforecast,  but don't forget to check CMA availability (warning)

(warning) Please note: use the variable "target"  to write each hindcastDate on a separate file .

 

v1.0

Code Block
languagepy
#!/usr/bin/env python
import calendar
from ecmwfapi import ECMWFDataServer
server = ECMWFDataServer()

origin = "babj"
modelVersionDate = "2014-05-01"

def retrieve_data(hindcastDate):
    target = "target_s2s_%s_%s.grb" % (origin, hindcastDate)
    server.retrieve({
       "class": "s2",
       "dataset": "s2s",
       "date": modelVersionDate,
       "expver": "prod",
       "hdate": hindcastDate,
       "levtype": "sfc",
       "origin": origin,
       "param": "165",
       "step": "0",
       "stream": "enfh",
       "target": target,
       "time": "00",
       "type": "cf",
})

hindcastYearStart = 1994
hindcastYearEnd = 2014
hindcastMonthStart = 1
hindcastMonthrEndhindcastMonthEnd = 12

for hindcastYear in list(range(hindcastYearStart, hindcastYearEnd + 1)):
    for hindcastMonth in list(range(hindcastMonthStart, hindcastMonthrEnd + 1)):
        numberOfDays = calendar.monthrange(hindcastYear, hindcastMonth)[1]
        for hindcastDay in list(range(1, numberOfDays + 1)):
            hindcastDate = '%04d%02d%02d' % (hindcastYear, hindcastMonth, hindcastDay)
            retrieve_data(hindcastDate) 

v1.1

Code Block
languagepy
#!/usr/bin/env python
import calendar
from ecmwfapi import ECMWFDataServer
server = ECMWFDataServer()

origin = "babj"
modelVersionDate = "2014-05-01"

hindcastYearStart = 1994
hindcastYearEnd = 2014
hindcastMonthStart = 1
hindcastMonthrEnd = 12

def retrieve_cma_reforecast():
    for hindcastYear in list(range(hindcastYearStart, hindcastYearEnd + 1)):
        for hindcastMonth in list(range(hindcastMonthStart, hindcastMonthrEndhindcastMonthEnd + 1)):
            numberOfDays = calendar.monthrange(hindcastYear, hindcastMonth)[1]
            for hindcastDay in list(range(1, numberOfDays + 1)):
                hindcastDate = '%04d%02d%02d' % (
                    hindcastYear, hindcastMonth, hindcastDay)
                cma_reforecast_request(hindcastDate)

def cma_reforecast_request(hindcastDate):
    target = "target_s2s_%s_%s.grb" % (origin, hindcastDate)
    server.retrieve({
        "class": "s2",
        "dataset": "s2s",
        "date": modelVersionDate,
        "expver": "prod",
        "hdate": hindcastDate,
        "levtype": "sfc",
        "origin": origin,
        "param": "165",
        "step": "0",
        "stream": "enfh",
        "target": target,
        "time": "00",
        "type": "cf",
    })

if __name__ == '__main__':
    retrieve_cma_reforecast()

...