Table of Contents
*************** This page is under construction! *****************
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 one hindcastDate
Info |
---|
The request below is for all members of the perturbed forecast, 10 meter U and V wind components, for all time-steps and for model version 2014-05-01 |
...
Code Block | ||
---|---|---|
| ||
#!/usr/bin/env python from ecmwfapi import ECMWFDataServer modelVersionDate = "2014-05-01" # This is the first model version for CMA (babj) hindcastDate = "2014-04-01" # The selected hindcast date server = ECMWFDataServer() 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": "data.pf.sfc", "time": "00", "number": "1/2/3", "type": "pf", }) |
...
A CMA reforecast request for
...
all the available hindcastDates
Info |
---|
|
...
|
v1.0
Code Block | ||
---|---|---|
| ||
#!/usr/bin/env python import calendar from ecmwfapi import ECMWFDataServer server = ECMWFDataServer() originhindcastYearStart = "babj"1994 modelVersionDatehindcastYearEnd = "2014-05-01"2014 hindcastMonthStart = 1 hindcastMonthEnd = 12 def retrieve_cma_datareforecast(hindcastDate,modelVersionDate,centre): targetfor = "target_s2s_%s_%s.grb" % centre % hindcastDatehindcastYear in list(range(hindcastYearStart, hindcastYearEnd + 1)): server.retrieve({ for hindcastMonth in list(range(hindcastMonthStart, hindcastMonthEnd + 1)): "class": "s2", numberOfDays = 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": centre, hindcastYear, hindcastMonth, hindcastDay + 1) "param": "165", "step": "0",cma_reforecast_request(hindcastDate) def cma_reforecast_request(hindcastDate): modelVersionDate = "2014-05-01" target = "stream": "enfh", data_s2s_%s.grb" % (hindcastDate) server.retrieve({ "targetclass": target"s2", "timedataset": "00s2s", "typedate": "cf"modelVersionDate, }) hindcastYearList = ["2012", "2013"] hindcastMonthList = ["08", "09"] hindcastDayList = ["01expver",: "02prod"], for hindcastYear in hindcastYearList: for hindcastMonth in hindcastMonthList:"hdate": hindcastDate, for hindcastDay in hindcastDayList: "levtype": "sfc", "origin": "babj", hindcastDate = hindcastYear + hindcastMonth + hindcastDay"param": "165/166", "step": "24/to/1440/by/24", retrieve_data(hindcastDate) |
Some further help for setting BoM hindcast dates
Code Block | ||||
---|---|---|---|---|
| ||||
for hindcastYear in ["2012stream",: "2013enfh"]:, for hindcastMonth in ["08", "09target"]: target, for hindcastDay in ["01" "06" "11" "16" "21" "26"]:"time": "00", "number": "1/2/3", "type": "pf", hindcastDate = hindcastYear+hindcastMonth+hindcastDay }) if __name__ == '__main__': retrieve_data(hindcastDate,"2014-01-01","ammc") |
...
cma_reforecast() |
Useful links
Info |
---|
*****************************************************************************************************