When accessing historic ecmwf forecasts downloaded from ecds, there is an error message “Index 85 out of bounds for length 85” (where the 85 varies with spatial domain specification). This prevents tools from Panoply from plotting. Access through xarray shows that the lat/lon in the file is very different from that specified in the download request. The issue does not show itself when downloading dwd (ensemble) data. Any idea? Would it be presumptuous for me to think that maybe there is a bug in the ecds somewhere?
Sample request (copied from the ecds web interface) is below:
import cdsapi
dataset = “tigge-forecasts”
request = {
“origin”: “ecmwf”,
“year”: “2026”,
“month”: “09”,
“day”: “01”,
“time”: “00:00”,
“level_type”: “single_level”,
“variable”: [“total_precipitation”],
“forecast_type”: “control_forecast”,
“leadtime_hour”: [
“0”,
“6”,
“12”,
“18”,
“24”,
“30”,
“36”,
“42”,
“48”,
“54”,
“60”,
“66”,
“72”,
“78”,
“84”,
“90”,
“96”,
“102”,
“108”,
“114”,
“120”,
“126”,
“132”,
“138”,
“144”,
“150”,
“156”,
“162”,
“168”,
“174”,
“180”,
“186”,
“192”,
“198”,
“204”,
“210”,
“216”,
“222”,
“228”,
“234”,
“240”,
“246”,
“252”,
“258”,
“264”,
“270”,
“276”,
“282”,
“288”,
“294”,
“300”,
“306”,
“312”,
“318”,
“324”,
“330”,
“336”,
“342”,
“348”,
“354”,
“360”
],
“data_format”: “grib”,
“area”: [53, 3, 51, 5]
}
client = cdsapi.Client()
client.retrieve(dataset, request).download()