Page History

...

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) according to your needs,
In this way you can extend this request to download the whole S2S CMA reforecast, but don't forget to check CMA availability

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

Code Block

language	py

#!/usr/bin/env python
from ecmwfapi import ECMWFDataServer
server = ECMWFDataServer()
def retrieve_data(hindcastDate,modelVersionDate,centre):
    target = "target_s2s_%s_%s.grb" % centre % hindcastDate
    server.retrieve({
       "class": "s2",
       "dataset": "s2s",
       "date": modelVersionDate,
       "expver": "prod",
       "hdate": hindcastDate,
       "levtype": "sfc",
       "origin": centre,
       "param": "165",
       "step": "0",
       "stream": "enfh",
       "target": target,
       "time": "00",
       "type": "cf",
})

for hindcastYear in ["2012", "2013"]:
    for hindcastMonth in ["08", "09"]:
		 for hindcastDay in ["01" "02" ]:
            hindcastDate = hindcastYear+hindcastMonth+hindcastDay
            retrieve_data(hindcastDate,"2014-05-01","babj")

Code Block

language	py

#!/usr/bin/env python
from ecmwfapi import ECMWFDataServer
server = ECMWFDataServer()

origin = "babj"
modelVersionDate = "2014-05-01" 
hindcastDate = "2014-04-01"      
hindcastYearList = ["2012", "2013"] 
hindcastMonthList = ["08", "09"]
hindcastDayList = ["01" "02" ]

def retrieve_data(hindcastDate, modelVersionDate, origin):
    target = "target_s2s_%s_%s.grb" % centre % 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",
})


def iterate():
for hindcastYear in hindcastYearList :
    for hindcastMonth in hindcastMonthList :
         for hindcastDay in hindcastDayList:
            hindcastDate = hindcastYear+hindcastMonth+hindcastDay
            retrieve_data(hindcastDate,"2014-05-01", "babj")

if __name__ == '__main__':
    print "The current settings server is:"
    print os.environ['WREP_SETTINGS_URL']
    iterate()

Code Block

language	py
title	CMA

for hindcastYear in ["2012", "2013"]:
    for hindcastMonth in ["08", "09"]:
         if hindcastMonth in ["01", "03", "05", "07", "08", "10", "12"]: hindcastDateEnd = "31"
         if hindcastMonth in [ "04", "06", "09", "11"]: hindcastDateEnd = "30"
         for hindcastDay in ["01" ... hindcastDateEnd ]:
            hindcastDate = hindcastYear+hindcastMonth+hindcastDay
            retrieve_data(hindcastDate,"2014-05-01","babj")

...

Space shortcuts

Page tree

Versions Compared

Old Version 83

New Version 84

Key