source: trunk/src/python_script/draw_map_tup.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
from pylab import *
from mpl_toolkits.basemap import Basemap
from mpl_toolkits.basemap import shiftgrid, cm




#################################
#carte cont ant
#amsua ch2


icont= (mpgrid_t60>=0.5)
tuppgrid_ta2ant=tuppgrid_ta260*icont
#tuppgrid_ta2ant=tuppgrid_ta260*mpgrid_t60
zz1=tuppgrid_ta2ant
tt1='tup cont amsua ch2'
t1='tup_cont_amsua_ch2.png'

#ssmis ch16
tuppgrid_t16ant=tuppgrid_t1660*mpgrid_t60
zz1=tuppgrid_t16ant
tt1='tup cont ssmis ch16'
t1='tup_cont_ssmis_ch16.png'

##########################"
zz1=tuppgrid_t
tt1='tup ssmis ch15'
t1='tup_ssmis_ch15.png'

zz1=tuppgrid_t16
tt1='tup ssmis ch16'
t1='tup_ssmis_ch16.png'

################
tupdif_15116=tuppgrid_t-tuppgrid_t16#dif presque cero
zz1=tupdif_15116
tt1='tup ssmis ch16'
t1='tup_ssmis_ch16.png'

zz1=tuppgrid_ta2
tt1='tup amsua ch2'
t1='tup_amsua_ch2.png'


zz1=tssigma_grid_ta2
tt1='ecart type tup amsua ch2'
t1='ecart_tipe_tup_amsua_ch2.png'
a1=0
a2=7
a3=0.1


zz1=tssigma_grid_t16
tt1='ecart type tup ssmis ch16'
t1='ecart_tipe_tup_sssmis_ch16.png'
a1=0
a2=7
a3=0.1


zz1=tssigma_grid_t
tt1='ecart type tup ssmis ch15'
t1='ecart_tipe_tup_sssmis_ch15.png'
a1=0
a2=7
a3=0.1

###################################
#ssmis
import string
import numpy as np
import matplotlib.pyplot as plt
from pylab import *
from mpl_toolkits.basemap import Basemap
from mpl_toolkits.basemap import shiftgrid, cm
import drawmap

tupssmisgird_max=np.array([tuppgrid_tss1.max(),tuppgrid_tss2.max(),tuppgrid_tss3.max(),tuppgrid_tss4.max(),tuppgrid_tss5.max(),tuppgrid_tss8.max(),tuppgrid_tss9.max(),tuppgrid_tss10.max(),tuppgrid_tss11.max(),tuppgrid_tss12.max(),tuppgrid_tss13.max(),tuppgrid_tss14.max(),tuppgrid_tss15.max(),tuppgrid_tss16.max(),tuppgrid_tss17.max(),tuppgrid_tss18.max()])

coef=0.6
cax1=0
cax2=300
cax3=1
my_cmap=cm.s3pcpn_l_r

zz=tuppgrid_tss1
txt='tup ssmis ch1'
outfile='tup_ssmis_ch1.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)


zz=tuppgrid_tss2
txt='tup ssmis ch2'
outfile='tup_ssmis_ch2.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tss3
txt='tup ssmis ch3'
outfile='tup_ssmis_ch3.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tss4
txt='tup ssmis ch4'
outfile='tup_ssmis_ch4.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tss5
txt='tup ssmis ch5'
outfile='tup_ssmis_ch5.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tss8
txt='tup ssmis ch8'
outfile='tup_ssmis_ch8.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tss9
txt='tup ssmis ch9'
outfile='tup_ssmis_ch9.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tss10
txt='tup ssmis ch10'
outfile='tup_ssmis_ch10.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tss11
txt='tup ssmis ch11'
outfile='tup_ssmis_ch11.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)


zz=tuppgrid_tss12
txt='tup ssmis ch12'
outfile='tup_ssmis_ch12.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tss13
txt='tup ssmis ch13'
outfile='tup_ssmis_ch13.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tss14
txt='tup ssmis ch14'
outfile='tup_ssmis_ch14.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tss15
txt='tup ssmis ch15'
outfile='tup_ssmis_ch15.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tss16
txt='tup ssmis ch16'
outfile='tup_ssmis_ch16.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tss17
txt='tup ssmis ch17'
outfile='tup_ssmis_ch17.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tss18
txt='tup ssmis ch18'
outfile='tup_ssmis_ch18.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

#############################################
#amsub
tupamsubgrid_max=np.array([tuppgrid_tab1.max(),tuppgrid_tab2.max(),tuppgrid_tab3.max(),tuppgrid_tab4.max(),tuppgrid_tab5.max()])

coef=0.6
cax1=0
cax2=300
cax3=1
my_cmap=cm.s3pcpn_l_r

zz=tuppgrid_tab1
txt='tup amsub ch1'
outfile='tup_amsub_ch1.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)


zz=tuppgrid_tab2
txt='tup amsub ch2'
outfile='tup_amsub_ch2.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tab3
txt='tup amsub ch3'
outfile='tup_amsub_ch3.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tab4
txt='tup amsub ch4'
outfile='tup_amsub_ch4.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)

zz=tuppgrid_tab5
txt='tup amsub ch5'
outfile='tup_amsub_ch5.png'
drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)



###################################
#amsua

tupmax=np.array([tuppgrid_ta1.max(), tuppgrid_ta2.max(), tuppgrid_ta3.max(),tuppgrid_ta4.max(),tuppgrid_ta5.max(),tuppgrid_ta6.max(),tuppgrid_ta7.max(),tuppgrid_ta15.max()])

tupmin=np.array([tuppgrid_ta1.min(), tuppgrid_ta2.min(), tuppgrid_ta3.min(),tuppgrid_ta4.min(),tuppgrid_ta5.min(),tuppgrid_ta6.min(),tuppgrid_ta7.min(),tuppgrid_ta15.min()])

zz1=tuppgrid_ta1
tt1='tup amsua ch1'
t1='tup_amsua_ch1.png'

zz1=tupnngrid_ta2
tt1='nom données tup ssmis ch15'
t1='nom_données_tup_sssmis_ch15.png'

zz1=tuppgrid_ta3
tt1='tup amsua ch3'
t1='tup_amsua_ch3.png'


zz1=tuppgrid_ta4
tt1='tup amsua ch4'
t1='tup_amsua_ch4.png'

import drawmap

zz=tuppgrid_ta1
txt='tup amsua ch1'
outfile='tup_amsua_ch11.png'
coef=0.6
cax1=0
cax2=230
cax3=1
my_cmap=cm.s3pcpn_l_r

zz=tuppgrid_ta2
txt='tup amsua ch2'
outfile='tup_amsua_ch22.png'

zz=tuppgrid_ta3
txt='tup amsua ch3'
outfile='tup_amsua_ch33.png'

zz=tuppgrid_ta4
txt='tup amsua ch4'
outfile='tup_amsua_ch44.png'

zz=tuppgrid_ta5
txt='tup amsua ch5'
outfile='tup_amsua_ch55.png'

zz=tuppgrid_ta6
txt='tup amsua ch6'
outfile='tup_amsua_ch66.png'

zz=tuppgrid_ta7
txt='tup amsua ch7'
outfile='tup_amsua_ch77.png'

zz=tuppgrid_ta15
txt='tup amsua ch15'
outfile='tup_amsua_ch1515.png'

drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef)





# 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))
#, clevs, cmap=my_cmap) 
clevs=arange(a1,a2,a3)#star, stop, step
cs=map.contourf(xii,yii,zz1, clevs, cmap=cm.s3pcpn_l_r)
#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, 30), labels=[0, 0, 0, 1])
map.drawparallels(np.arange(-90, 90, 10), labels=[1, 0, 0, 0])
#plt.show()
plt.savefig(t1)
close()