plt.figure()
m = Basemap()
m.drawcoastlines()
m.imshow(dat)
plt.show()