Table of Contents
What is the objective of this page?
Info |
---|
To help users to improve S2S CMA MARS requests performance via the WebAPI.
|
How is the S2S reforecast data organised in MARS?
Info |
---|
In general it is organised, as a huge tree, with the indentation below, showing different levels down that tree:
|
What would be the natural way to group requests?
Info |
---|
The idea is to request as much data as possible from the same tape file. The natural way to group requests would be:
|
Best practise to iterate over all hindcastDates of several hindcastYears for CMA
Info |
---|
|
Info |
---|
for hindcastYear in hindcastYears |
Web-API examples:
A CMA reforecast request for
...
all the available hindcastDates
Info |
---|
Please note:
|
Code Block | ||
---|---|---|
| ||
#!/usr/bin/env python import calendar from ecmwfapi import ECMWFDataServer server = ECMWFDataServer() origin = "babj" modelVersionDate = "2014-05-01" hindcastDate = "2014-04-01 def retrieve_cma_reforecast(): """ # The selected hindcast date server = ECMWFDataServer() server.retrieve({ "class": "s2", "dataset": "s2s", "date": modelVersionDate, A function to demonstrate how to iterate efficiently over all hindcastYears, hindcastMonths etc "expver": "prod", for a "hdate": hindcastDate, particular cma_reforecast_request. "levtype": "sfc", "origin": "babj", Change the "param": "165/166", "step": "24/to/1440/by/24",variables below to adapt the iteration to your needs "stream": "enfh", "target": "data.pf.sfc", hindcastYearStart = 1994 hindcastYearEnd "time": "00", = 2014 hindcastMonthStart "number":= "1/2/3", hindcastMonthEnd "type": "pf", }) |
A CMA reforecast request for several hindcastDates
Info | ||||
---|---|---|---|---|
Code Block | ||||
| ||||
#!/usr/bin/env python from ecmwfapi import ECMWFDataServer server = ECMWFDataServer() def retrieve_data(hindcastDate,modelVersionDate): target = "target_s2s_%s.grb" % hindcastDate server.retrieve({ = 12 for hindcastYear in list(range(hindcastYearStart, hindcastYearEnd + 1)): for hindcastMonth in list(range(hindcastMonthStart, hindcastMonthEnd + 1)): numberOfDays "class": "s2",= calendar.monthrange(hindcastYear, hindcastMonth)[1] "dataset": "s2s", for hindcastDay "date": modelVersionDate,in list(range(numberOfDays)): "expver": "prod", "hdate": hindcastDate, = '%04d%02d%02d' % ( "levtype": "sfc", "origin": "babj", "param": "165", hindcastYear, hindcastMonth, hindcastDay + 1) "step": "0", "stream": "enfh", target = "target": target,%s%s.grb" % (origin, hindcastDate) "time": "00", "type": "cf", } cma_reforecast_request(hindcastDate, target) for hindcastYear in ["2012", "2013"]def cma_reforecast_request(hindcastDate, target): for hindcastMonth in ["08", "09"]: for hindcastDay in ["01" "02" ]:""" A CMA reforecast request. hindcastDateChange = hindcastYear+hindcastMonth+hindcastDay the keywords below to adapt it to your needs. retrieve_data(hindcastDate,"2014-05-01") | ||||
Code Block | ||||
| ||||
for hindcastYear in ["2012", "2013"]: for hindcastMonth in ["08", "09"]: """ server.retrieve({ if hindcastMonth in ["01class",: "03s2", "05", "07", "08", "10", "12"]: hindcastDateEnd = "31" "dataset": "s2s", "date": modelVersionDate, if hindcastMonth in [ "04expver",: "06prod", "09", "11"]: hindcastDateEnd = "30" "hdate": hindcastDate, for hindcastDay in ["01levtype" ... hindcastDateEnd ]:: "sfc", "origin": origin, hindcastDate = hindcastYear+hindcastMonth+hindcastDay "param": "165/166", retrieve_data(hindcastDate,"2014-05-01") | ||||
Code Block | ||||
| ||||
for hindcastYear in ["2012", "2013"]: "step": "24/to/1440/by/24", for hindcastMonth in ["08stream",: "09enfh"]:, for hindcastDay in ["01" "06" "11" "16" "21" "26"]:"target": target, "time": "00", "number": "1/2/3", hindcastDate = hindcastYear+hindcastMonth+hindcastDay "type": "pf", }) if __name__ == '__main__': retrieve_cma_data(hindcastDate,"2014-01-01") |
...
reforecast() |
Useful links
Info |
---|
...