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

#mod gabi para graficar AMSUA CH2
#!/usr/bin/env python
# -*- coding: utf-8 -*-
import string
import numpy as np
import matplotlib.pyplot as plt
import ffgrid2
from pylab import *
from mpl_toolkits.basemap import Basemap
from mpl_toolkits.basemap import shiftgrid, cm
#import netCDF4

#x=monthly_lon
#y=monthly_lat
#z=monthly_outz

zz1=zzpgrid_3n
zz2=sigma_grid_3n
zz3=nngrid_3n

t1='amsua_ch2_tb.png'
tt1='amsua_ch2 tb'
t2='amsua_ch2_ecart_type.png'
tt2='amsua_ch2 ecart type'
t3='amsua_ch2_nom_donnes_tb.png'
tt3='amsua ch2 nom donnes tb'

t1='amsua_ch2_pos3_tb.png'
tt1='amsua_ch2 pos3 tb'
t2='amsua_ch2_pos3_ecart_type.png'
tt2='amsua_ch2 pos3 ecart type'
t3='amsua_ch2_pos3_nom_donnes_tb.png'
tt3='amsua ch2 pos3 nom donnes tb'


# on recupere les coordonnes xy de la projection orth standard
m1 = Basemap(projection='ortho', lat_0 = -90, lon_0 = 0,
resolution = 'l')

xii, yii = m1(*np.meshgrid(xvec,yvec))
width = m1.urcrnrx - m1.llcrnrx
height = m1.urcrnry - m1.llcrnry

coef = 0.6
width = width*coef
height = height*coef

# on fait une nouvelle projection en zoomant sur l'antarctique

map = Basemap(projection='ortho',lon_0=0,lat_0=-90,resolution='l',\

llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height)
xii, yii = map(*np.meshgrid(xvec,yvec))
cs=m.contourf(xii,yii,zz1)
#, clevs, cmap=my_cmap) 
cs=map.pcolormesh(xii,yii,zz1,cmap=cm.s3pcpn_l_r)
#cs=map.pcolormesh(xii,yii,zz2)
#cmap=cm.s3pcpn_l_r)
#sstanom)
#s3pcpn_l_r)
cbar =colorbar(cs)
plt.title(tt1)

# draw coastlines, country boundaries, fill continents.
map.drawcoastlines(linewidth=1)
# draw the edge of the map projection region (the projection limb)
map.drawmapboundary()
# draw lat/lon grid lines every 30 degrees.
map.drawmeridians(np.arange(0, 360, 10))
map.drawparallels(np.arange(-90, 90, 10))
#plt.show()

plt.savefig(t1)
close()

#carte ecart type
zz21 = sigma_grid_t
zz22 = sigma_grid_1n
zz23 = sigma_grid_2n
zz24 = sigma_grid_3n
tt21='amsua ch2 t ecart type'
tt22='amsua ch2 1n ecart type'
tt23='amsua ch2 2n ecart type'
tt24='amsua ch2 3n ecart type'
t='amsua_ch2_ecart_type.png'

figure(1)
subplot(221)
m1 = Basemap(projection='ortho', lat_0 = -90, lon_0 = 0,
resolution = 'l')

xii, yii = m1(*np.meshgrid(xvec,yvec))
width = m1.urcrnrx - m1.llcrnrx
height = m1.urcrnry - m1.llcrnry

coef = 0.6
width = width*coef
height = height*coef

# on fait une nouvelle projection en zoomant sur l'antarctique
map = Basemap(projection='ortho',lon_0=0,lat_0=-90,resolution='l',\
llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height)
xii, yii = map(*np.meshgrid(xvec,yvec))
clevs=arange(0,40,0.1)#star, stop, step
cs=map.contourf(xii,yii,zz21, clevs, cmap=cm.s3pcpn_l)
#cs=map.pcolormesh(xii,yii,zz2)
#cmap=cm.s3pcpn_l_r)
#sstanom)
#s3pcpn_l_r)
cbar =colorbar(cs)
plt.title(tt21)
# draw coastlines, country boundaries, fill continents.
map.drawcoastlines(linewidth=1)
# draw the edge of the map projection region (the projection limb)
map.drawmapboundary()
# draw lat/lon grid lines every 30 degrees.
map.drawmeridians(np.arange(0, 360, 10))
map.drawparallels(np.arange(-90, 90, 10))

subplot(222)
# on fait une nouvelle projection en zoomant sur l'antarctique
map = Basemap(projection='ortho',lon_0=0,lat_0=-90,resolution='l',\
llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height)
xii, yii = map(*np.meshgrid(xvec,yvec))
clevs=arange(0,40,0.1)#star, stop, step
cs=map.contourf(xii,yii,zz22, clevs, cmap=cm.s3pcpn_l)
#cs=map.pcolormesh(xii,yii,zz2)
#cmap=cm.s3pcpn_l_r)
#sstanom)
#s3pcpn_l_r)
cbar =colorbar(cs)
plt.title(tt22)
# draw coastlines, country boundaries, fill continents.
map.drawcoastlines(linewidth=1)
# draw the edge of the map projection region (the projection limb)
map.drawmapboundary()
# draw lat/lon grid lines every 30 degrees.
map.drawmeridians(np.arange(0, 360, 10))
map.drawparallels(np.arange(-90, 90, 10))


subplot(223)
# on fait une nouvelle projection en zoomant sur l'antarctique
map = Basemap(projection='ortho',lon_0=0,lat_0=-90,resolution='l',\
llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height)
xii, yii = map(*np.meshgrid(xvec,yvec))
clevs=arange(0,40,0.1)#star, stop, step
cs=map.contourf(xii,yii,zz23, clevs, cmap=cm.s3pcpn_l)
#cs=map.pcolormesh(xii,yii,zz2)
#cmap=cm.s3pcpn_l_r)
#sstanom)
#s3pcpn_l_r)
cbar =colorbar(cs)
plt.title(tt23)
# draw coastlines, country boundaries, fill continents.
map.drawcoastlines(linewidth=1)
# draw the edge of the map projection region (the projection limb)
map.drawmapboundary()
# draw lat/lon grid lines every 30 degrees.
map.drawmeridians(np.arange(0, 360, 10))
map.drawparallels(np.arange(-90, 90, 10))


subplot(224)
# on fait une nouvelle projection en zoomant sur l'antarctique
map = Basemap(projection='ortho',lon_0=0,lat_0=-90,resolution='l',\
llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height)
xii, yii = map(*np.meshgrid(xvec,yvec))
clevs=arange(0,40,0.1)#star, stop, step
cs=map.contourf(xii,yii,zz24, clevs, cmap=cm.s3pcpn_l)
#cs=map.pcolormesh(xii,yii,zz2)
#cmap=cm.s3pcpn_l_r)
#sstanom)
#s3pcpn_l_r)
cbar =colorbar(cs)
plt.title(tt24)
# draw coastlines, country boundaries, fill continents.
map.drawcoastlines(linewidth=1)
# draw the edge of the map projection region (the projection limb)
map.drawmapboundary()
# draw lat/lon grid lines every 30 degrees.
map.drawmeridians(np.arange(0, 360, 10))
map.drawparallels(np.arange(-90, 90, 10))

plt.show()
plt.savefig(t)
close()






##########################"
#carte nom de donnees
m1 = Basemap(projection='ortho', lat_0 = -90, lon_0 = 0,
resolution = 'l')

xii, yii = m1(*np.meshgrid(xvec,yvec))
width = m1.urcrnrx - m1.llcrnrx
height = m1.urcrnry - m1.llcrnry

coef = 0.6
width = width*coef
height = height*coef

# on fait une nouvelle projection en zoomant sur l'antarctique

map = Basemap(projection='ortho',lon_0=0,lat_0=-90,resolution='l',\

llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height)
xii, yii = map(*np.meshgrid(xvec,yvec))

evels=[300, 290, 280, 270, 260, 250, 240, 230, 220, 210, 200, 190, 180, 170, 160, 150, 140, 130, 120, 110, 100]
cs=map.pcolormesh(xii,yii,zz3)
#cmap=cm.s3pcpn_l_r)
#sstanom)
#s3pcpn_l_r)
cbar =colorbar(cs)
plt.title(tt3)

# draw coastlines, country boundaries, fill continents.
map.drawcoastlines(linewidth=1)
# draw the edge of the map projection region (the projection limb)
map.drawmapboundary()
# draw lat/lon grid lines every 30 degrees.
map.drawmeridians(np.arange(0, 360, 10))
map.drawparallels(np.arange(-90, 90, 10))
#plt.show()

plt.savefig(t3)
close()