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 data organised in MARS?
...
In general it is organised, as a huge tree, with the indentation below, showing different levels down that tree:
- centre (CMA, ECMWF, NCEP, JMA, ...)
- realtime or reforecast
- type of data (control forecast or perturbed forecast)
- type of level (single level or pressure level or potential temperature)
- model version date (2014-05-01 or ...)
- hindcast dates (2014-01-01 or 2014-01-02 or 2014-01-03, ...)
- time-steps
- members (for perturbed forecast)
- levels (for pl or pt)
- parameters
- levels (for pl or pt)
- members (for perturbed forecast)
- time-steps
- hindcast dates (2014-01-01 or 2014-01-02 or 2014-01-03, ...)
- model version date (2014-05-01 or ...)
- type of level (single level or pressure level or potential temperature)
- type of data (control forecast or perturbed forecast)
- realtime or reforecast
...
What would be the natural way to group requests?
Info |
---|
Following the previous paragraph, the natural way to group requests would be:
|
Best practices:
To iterate over several hindcastDates for a CMA request
...
title | The main idea in brief: |
---|
...
Best practise to iterate over all hindcastDates of
...
several hindcastYears
...
for
...
CMA
...
Info | |
---|---|
The main idea in brief: |
|
Info | title | The main idea in brief:
---|
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", A function to demonstrate how to iterate efficiently over all hindcastYears, hindcastMonths etc for a particular cma_reforecast_request. Change the variables below to adapt the iteration to your needs "dataset": "s2s", "date": modelVersionDate, hindcastYearStart = 1994 "expver": "prod",hindcastYearEnd = 2014 "hdate": hindcastDate, hindcastMonthStart = 1 hindcastMonthEnd = 12 "levtype": "sfc", "origin": "babj",for hindcastYear in list(range(hindcastYearStart, hindcastYearEnd + 1)): "param": "165", "step": "0", for hindcastMonth in list(range(hindcastMonthStart, hindcastMonthEnd + 1)): "stream": "enfh", numberOfDays "target": "data.cf.sfc", = calendar.monthrange(hindcastYear, hindcastMonth)[1] "time": "00", "type": "cf", }) |
A CMA reforecast request for several hindcastDates (iterate efficiently)
Info | ||
---|---|---|
Code Block | ||
| ||
#!/usr/bin/env python from ecmwfapi import ECMWFDataServer server = ECMWFDataServer() def retrieve_data(hindcastDate): target = "target_s2s_%s.grb" % hindcastDatefor hindcastDay in list(range(numberOfDays)): hindcastDate = '%04d%02d%02d' % ( hindcastYear, hindcastMonth, hindcastDay + 1) target = "%s%s.grb" % (origin, hindcastDate) cma_reforecast_request(hindcastDate, target) def cma_reforecast_request(hindcastDate, target): """ A CMA reforecast request. Change the keywords below to adapt it to your needs. """ server.retrieve({ "class": "s2", "dataset": "s2s", "date": "2014-05-01"modelVersionDate, "expver": "prod", "hdate": hindcastDate, "levtype": "sfc", "origin": "babj"origin, "param": "165/166", "step": "024/to/1440/by/24", "stream": "enfh", "target": target, "time": "00", "typenumber": "cf1/2/3", }) for hindcastYear in ["2012", "2013"]: for hindcastMonth in ["08", "09"]: for hindcastDay in ["01", "02"]:"type": "pf", }) if hindcastDate = hindcastYear+hindcastMonth+hindcastDay __name__ == '__main__': retrieve_datacma_reforecast(hindcastDate)) |
Useful links
Info |
---|