我使用python 3.6绘制CMIP5的降水数据,我下载的文件是netCDF4文件。我已经在另一个类似的文件上使用了这个代码,它运行得很好,所以我不确定是什么问题。我没有收到任何有关此代码的错误消息,它只显示一个世界地图,它应该是各种颜色的所有颜色。在此文件中找到的变量是time,time_bnds,lat,lat_bnds,lon,lon_bnds和prc。 prc是降水变量,也是我对绘图感兴趣的变量。任何想法都会有所帮助,谢谢!
这是我的代码
from mpl_toolkits.basemap import Basemap, cm
from netCDF4 import Dataset as NetCDFFile
import matplotlib.pyplot as plt
nc = NetCDFFile('filename.nc','r')
p = nc.variables['prc']
data = p[:,:,0]
fig = plt.figure(figsize=(8,8))
ax = fig.add_axes([0.1,0.1,0.8,0.8])
m = Basemap(projection='cyl',lon_0=180,lat_0=0,resolution='l')
m.drawcoastlines()
m.drawstates()
m.drawcountries()
ny = data.shape[0]; nx = data.shape[1]
lons, lats = m.makegrid(nx,ny)
x,y = m(lons, lats) # compute map proj coordinates.
cs=plt.contourf(x,-y,data,range(0,1000,10),cmap=cm.s3pcpn,latlon=True)
cbar = m.colorbar(cs,location='bottom',pad="5%")
cbar.set_label('mm')
plt.show()