m=Basemap(projection='cyl', llcrnrlon=-180, \
urcrnrlon=180.,llcrnrlat=lat.min(),urcrnrlat=lat.max(), \
resolution='c')
x, y = m(lon,lat)
cs = m.pcolormesh(x,y,data,shading='flat',cmap=plt.cm.hot, norm=colors.LogNorm())
m.drawcoastlines()
m.drawmapboundary()
m.drawparallels(np.arange(-90.,120.,30.),labels=[1,0,0,0])
m.drawmeridians(np.arange(-180.,180.,60.),labels=[0,0,0,1])
plt.colorbar(cs,orientation='vertical', shrink=0.5)
plt.title('CAMS GFAS fire radiative power (W / m2)') # Set the name of the variable to plot
plt.savefig('gfas_05.png') # |