retrieve_seas5.py

import cdsapi

c = cdsapi.Client()

r = c.retrieve(
    'seasonal-original-single-levels',
    {
        'variable': 'total_precipitation',
        'originating_centre': 'ecmwf',
        'system': '5',
        'year': [str(y) for y in range(1981, 2022)],
        'month': [
            '05', '06', '07',
            '08', '09',
        ],
        'day': '01',
        'area': [
            39, 66, 6,
            100,
        ],
        'leadtime_hour': [
            '24', '48', '72',
            '96', '120', '144',
            '168', '192', '216',
            '240', '264', '288',
            '312', '336', '360',
            '384', '408', '432',
            '456', '480', '504',
            '528', '552', '576',
            '600', '624', '648',
            '672', '696', '720',
            '744', '768', '792',
            '816', '840', '864',
            '888', '912', '936',
            '960', '984', '1008',
            '1032', '1056', '1080',
            '1104', '1128', '1152',
            '1176', '1200', '1224',
            '1248', '1272', '1296',
            '1320', '1344', '1368',
            '1392', '1416', '1440',
        ],
        'format': 'netcdf'
    })
r.download('data/seas5_forecasts.nc')