http://blog.csdn.net/wei_guo_xd/article/details/73729581
http://blog.csdn.net/aimer_chen/article/details/78528000
basemap install
http://blog.csdn.net/aimer_chen/article/details/78528000
http://blog.sina.com.cn/s/blog_662dcb820102w21u.html
#coding=utf-8
__author__ = 'kong'
import matplotlib.pyplot as plt
import numpy as np
from mpl_toolkits.basemap import Basemap, cm
import netCDF4 as nc
filename="/home/cao/nc/wrfout_d01_2015-12-13_01_00_00"
data=nc.Dataset(filename)
for i in data.variables.keys():
print(i)
lat= data.variables['XLAT_U'][:].squeeze()
lon= data.variables['XLONG_U'][:].squeeze()
u=data.variables['U'][:].squeeze()
idate=data.variables['Times']
m=Basemap(resolution='l',area_thresh=10000,projection='cyl',llcrnrlon=70,urcrnrlon=150,llcrnrlat=10,urcrnrlat=60)
fig1 = plt.figure()
ax = fig1.add_axes([0.1,0.1,0.8,0.8])
#clevls_hgt=[-80,-60,-40,-20,0,20,40,60,80]
x, y = m(lon, lat)
m.drawparallels(np.arange(0,60,10.),labels=[1,1,0,0],fontsize=15)#longitude
m.drawmeridians(np.arange(70,150,10),labels=[0,0,0,1],fontsize=15)#latitude
#curve=m.contour(x,y,u[0,:,:],colors='k')
CS2=m.contourf(x,y,u[0,:,:])
m.colorbar(CS2)
m.drawcoastlines(linewidth=0.2)
plt.title('U',size=20)
plt.show()