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() |
---|