Plt.title('CAMS AOD forecast') # Set the name of the variable to plot Plt.colorbar(cs,orientation='vertical', shrink=0.5) Urcrnrlon=180.,llcrnrlat=lats.min(),urcrnrlat=lats.max(), \Ĭs = m.pcolormesh(x,y,data,shading='flat',cmap=plt.cm.gist_stern_r) M = Basemap(projection='cyl', llcrnrlon=-180, \ Grid_lon, grid_lat = np.meshgrid(lons, lats) #regularly spaced 2D grid Lats = np.linspace(float(grb), \ĭata, lons = shiftgrid(180., data, lons, start=False) Lons = np.linspace(float(grb), \įloat(grb), int(grb) ) # need to shift data grid longitudes from (0.360) to (-180.180) Grib = 'cams_aod.grib' # Set the file name of your input GRIB file From mpl_toolkits.basemap import shiftgrid
0 Comments
Leave a Reply. |