Базовая карта не отображает данные GFS

Question

У меня есть следующий код:

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
import pygrib

file = pygrib.open('./gfs_3_20180220_0000_000.grb2')

z1000Pa = file[36].data()[0]
lon = file[36].data()[1]
lat = file[36].data()[2]

m = Basemap(llcrnrlon = 180, llcrnrlat = 15., urcrnrlon = 360., 
urcrnrlat = 80., resolution='i', projection = 'merc')


m.drawcoastlines()
m.drawparallels(np.arange(15,80,15),labels=[1,0,0,0], linewidth=0.0)
m.drawmeridians(np.arange(180,360,30),labels=[0,0,0,1], linewidth=0.0)

X = np.meshgrid(lon)
Y = np.meshgrid(lat)

cs = m.contourf(Y, X, z1000Pa, np.arange(0,6000,40), colors='red')
plt.clabel(cs, fmt = '%1.f')

Итак, я пытаюсь нанести z1000Pa на карту базовой карты, но я получил следующую ошибку:

AttributeError: 'list' object has no attribute 'shape'

Я пробовал разные вещи, но это не работает.