CDS API request on ERA5 on model levels return 1-dimensional grib data

Hi,

I am trying to retrieve some 3-D met vars from ERA5 on 137 model levels. However, using the CDS API script pulls data in 1-D as (410240,) using pygrib in Python, which should be (640,1280). 

I attached the cds script below. Any comments are appreciated. Thanks! Xia

#!/usr/bin/env python
import cdsapi
c = cdsapi.Client()
c.retrieve('reanalysis-era5-complete', {
    'class': 'ea',
    'date': '2020-02-04',
    'expver': '1',
    'levelist': '1/2/3/4/5/6/7/8/9/10/11/12/13/14/15/16/17/18/19/20/21/22/23/24/25/26/27/28/29/30/31/32/33/34/35/36/37/38/39/40/41/42/43/44/45/46/47/48/49/50/51/52/53/54/55/56/57/58/59/60/61/62/63/64/65/66/67/68/69/70/71/72/73/74/75/76/77/78/79/80/81/82/83/84/85/86/87/88/89/90/91/92/93/94/95/96/97/98/99/100/101/102/103/104/105/106/107/108/109/110/111/112/113/114/115/116/117/118/119/120/121/122/123/124/125/126/127/128/129/130/131/132/133/134/135/136/137',
    'levtype': 'ml',
    'param': '130',
    'stream': 'oper',
    'time': '00:00:00',
    'type': 'an',
}, 'output')