[6] | 1 | #!/usr/bin/env python |
---|
| 2 | # -*- coding: utf-8 -*- |
---|
| 3 | import string |
---|
| 4 | import numpy as np |
---|
| 5 | import matplotlib.pyplot as plt |
---|
| 6 | import ffgrid2 |
---|
| 7 | from pylab import * |
---|
| 8 | from mpl_toolkits.basemap import Basemap |
---|
| 9 | from mpl_toolkits.basemap import shiftgrid, cm |
---|
| 10 | #import netCDF4 |
---|
| 11 | |
---|
| 12 | x=monthly_lon |
---|
| 13 | y=monthly_lat |
---|
| 14 | z=monthly_outz |
---|
| 15 | |
---|
| 16 | def draw(z, x, y, outfile, txt, cax1,cax2,cax3, my_cmap): |
---|
| 17 | |
---|
| 18 | m = Basemap(llcrnrlon=-180,urcrnrlon=180,llcrnrlat=-90,urcrnrlat=90,projection='cyl',resolution='c',fix_aspect=True) |
---|
| 19 | #m = Basemap(projection='ortho',lon_0=0,lat_0=-90,resolution='c') |
---|
| 20 | #m = Basemap(resolution='c',projection='sinu',lon_0=60) |
---|
| 21 | #m = Basemap(llcrnrlon=-180,urcrnrlon=180,llcrnrlat=-90,urcrnrlat=90,projection='cyl',resolution='i',fix_aspect=True) |
---|
| 22 | #m = Basemap(llcrnrlon=-180,urcrnrlon=180,llcrnrlat=-60,urcrnrlat=60,projection='cyl',resolution='c',fix_aspect=True) |
---|
| 23 | |
---|
| 24 | #ax = plt.gca() |
---|
| 25 | xii, yii = m(*np.meshgrid(x,y)) |
---|
| 26 | #ax = plt.gca() |
---|
| 27 | #clevs=arange(0,0.8,0.005) |
---|
| 28 | clevs=arange(cax1,cax2,cax3) |
---|
| 29 | #boxColors= ['darkkhaki','royalblue'] |
---|
| 30 | |
---|
| 31 | m.drawcoastlines(linewidth=1) |
---|
| 32 | #m.drawlsmask(land_color='white',ocean_color='#F4EAF1', lakes=False) |
---|
| 33 | cs=m.contourf(xii,yii,z,20) # emissivity |
---|
| 34 | cbar =colorbar(cs) |
---|
| 35 | cbar.set_label(txt) |
---|
| 36 | #cbar.set_clim(0.65,1) |
---|
| 37 | #m.drawmapboundary(fill_color='aqua') |
---|
| 38 | #m.drawparallels(np.arange(-90.,90.,10.), labels=[1,0,0,0]) |
---|
| 39 | #m.drawmeridians(np.arange(-180.,180.,30.), labels=[0,0,0,1]) |
---|
| 40 | m.drawparallels(np.arange(-90.,90.,45.), color='#909090', linewidth=1.0, zorder=None, dashes=[1, 1]) |
---|
| 41 | #m.drawmeridians(np.arange(-180.,180.,50.)) |
---|
| 42 | # stations: Thulé, Resolute, Longyearbyen, lqaluit, murmonsk, arkangelek |
---|
| 43 | |
---|
| 44 | plt.savefig(outfile) |
---|
| 45 | #close() |
---|