source: trunk/src/python_script/draw_map_mod_gc1.py @ 6

Last change on this file since 6 was 6, checked in by gaclod, 12 years ago

add GC python scripts

  • Property svn:executable set to *
File size: 1.6 KB
Line 
1#!/usr/bin/env python
2# -*- coding: utf-8 -*-
3import string
4import numpy as np
5import matplotlib.pyplot as plt
6import ffgrid2
7from pylab import *
8from mpl_toolkits.basemap import Basemap
9from mpl_toolkits.basemap import shiftgrid, cm
10#import netCDF4
11
12x=monthly_lon
13y=monthly_lat
14z=monthly_outz
15
16def 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()
Note: See TracBrowser for help on using the repository browser.