[6] | 1 | #mod gabi para graficar AMSUA CH2 |
---|
| 2 | #!/usr/bin/env python |
---|
| 3 | # -*- coding: utf-8 -*- |
---|
| 4 | import string |
---|
| 5 | import numpy as np |
---|
| 6 | import matplotlib.pyplot as plt |
---|
| 7 | from pylab import * |
---|
| 8 | from mpl_toolkits.basemap import Basemap |
---|
| 9 | from mpl_toolkits.basemap import shiftgrid, cm |
---|
| 10 | |
---|
| 11 | |
---|
| 12 | |
---|
| 13 | |
---|
| 14 | ################################# |
---|
| 15 | #carte cont ant |
---|
| 16 | #amsua ch2 |
---|
| 17 | |
---|
| 18 | |
---|
| 19 | icont= (mpgrid_t60>=0.5) |
---|
| 20 | tuppgrid_ta2ant=tuppgrid_ta260*icont |
---|
| 21 | #tuppgrid_ta2ant=tuppgrid_ta260*mpgrid_t60 |
---|
| 22 | zz1=tuppgrid_ta2ant |
---|
| 23 | tt1='tup cont amsua ch2' |
---|
| 24 | t1='tup_cont_amsua_ch2.png' |
---|
| 25 | |
---|
| 26 | #ssmis ch16 |
---|
| 27 | tuppgrid_t16ant=tuppgrid_t1660*mpgrid_t60 |
---|
| 28 | zz1=tuppgrid_t16ant |
---|
| 29 | tt1='tup cont ssmis ch16' |
---|
| 30 | t1='tup_cont_ssmis_ch16.png' |
---|
| 31 | |
---|
| 32 | ##########################" |
---|
| 33 | zz1=tuppgrid_t |
---|
| 34 | tt1='tup ssmis ch15' |
---|
| 35 | t1='tup_ssmis_ch15.png' |
---|
| 36 | |
---|
| 37 | zz1=tuppgrid_t16 |
---|
| 38 | tt1='tup ssmis ch16' |
---|
| 39 | t1='tup_ssmis_ch16.png' |
---|
| 40 | |
---|
| 41 | ################ |
---|
| 42 | tupdif_15116=tuppgrid_t-tuppgrid_t16#dif presque cero |
---|
| 43 | zz1=tupdif_15116 |
---|
| 44 | tt1='tup ssmis ch16' |
---|
| 45 | t1='tup_ssmis_ch16.png' |
---|
| 46 | |
---|
| 47 | zz1=tuppgrid_ta2 |
---|
| 48 | tt1='tup amsua ch2' |
---|
| 49 | t1='tup_amsua_ch2.png' |
---|
| 50 | |
---|
| 51 | |
---|
| 52 | zz1=tssigma_grid_ta2 |
---|
| 53 | tt1='ecart type tup amsua ch2' |
---|
| 54 | t1='ecart_tipe_tup_amsua_ch2.png' |
---|
| 55 | a1=0 |
---|
| 56 | a2=7 |
---|
| 57 | a3=0.1 |
---|
| 58 | |
---|
| 59 | |
---|
| 60 | zz1=tssigma_grid_t16 |
---|
| 61 | tt1='ecart type tup ssmis ch16' |
---|
| 62 | t1='ecart_tipe_tup_sssmis_ch16.png' |
---|
| 63 | a1=0 |
---|
| 64 | a2=7 |
---|
| 65 | a3=0.1 |
---|
| 66 | |
---|
| 67 | |
---|
| 68 | zz1=tssigma_grid_t |
---|
| 69 | tt1='ecart type tup ssmis ch15' |
---|
| 70 | t1='ecart_tipe_tup_sssmis_ch15.png' |
---|
| 71 | a1=0 |
---|
| 72 | a2=7 |
---|
| 73 | a3=0.1 |
---|
| 74 | |
---|
| 75 | ################################### |
---|
| 76 | #ssmis |
---|
| 77 | import string |
---|
| 78 | import numpy as np |
---|
| 79 | import matplotlib.pyplot as plt |
---|
| 80 | from pylab import * |
---|
| 81 | from mpl_toolkits.basemap import Basemap |
---|
| 82 | from mpl_toolkits.basemap import shiftgrid, cm |
---|
| 83 | import drawmap |
---|
| 84 | |
---|
| 85 | 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()]) |
---|
| 86 | |
---|
| 87 | coef=0.6 |
---|
| 88 | cax1=0 |
---|
| 89 | cax2=300 |
---|
| 90 | cax3=1 |
---|
| 91 | my_cmap=cm.s3pcpn_l_r |
---|
| 92 | |
---|
| 93 | zz=tuppgrid_tss1 |
---|
| 94 | txt='tup ssmis ch1' |
---|
| 95 | outfile='tup_ssmis_ch1.png' |
---|
| 96 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 97 | |
---|
| 98 | |
---|
| 99 | zz=tuppgrid_tss2 |
---|
| 100 | txt='tup ssmis ch2' |
---|
| 101 | outfile='tup_ssmis_ch2.png' |
---|
| 102 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 103 | |
---|
| 104 | zz=tuppgrid_tss3 |
---|
| 105 | txt='tup ssmis ch3' |
---|
| 106 | outfile='tup_ssmis_ch3.png' |
---|
| 107 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 108 | |
---|
| 109 | zz=tuppgrid_tss4 |
---|
| 110 | txt='tup ssmis ch4' |
---|
| 111 | outfile='tup_ssmis_ch4.png' |
---|
| 112 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 113 | |
---|
| 114 | zz=tuppgrid_tss5 |
---|
| 115 | txt='tup ssmis ch5' |
---|
| 116 | outfile='tup_ssmis_ch5.png' |
---|
| 117 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 118 | |
---|
| 119 | zz=tuppgrid_tss8 |
---|
| 120 | txt='tup ssmis ch8' |
---|
| 121 | outfile='tup_ssmis_ch8.png' |
---|
| 122 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 123 | |
---|
| 124 | zz=tuppgrid_tss9 |
---|
| 125 | txt='tup ssmis ch9' |
---|
| 126 | outfile='tup_ssmis_ch9.png' |
---|
| 127 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 128 | |
---|
| 129 | zz=tuppgrid_tss10 |
---|
| 130 | txt='tup ssmis ch10' |
---|
| 131 | outfile='tup_ssmis_ch10.png' |
---|
| 132 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 133 | |
---|
| 134 | zz=tuppgrid_tss11 |
---|
| 135 | txt='tup ssmis ch11' |
---|
| 136 | outfile='tup_ssmis_ch11.png' |
---|
| 137 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 138 | |
---|
| 139 | |
---|
| 140 | zz=tuppgrid_tss12 |
---|
| 141 | txt='tup ssmis ch12' |
---|
| 142 | outfile='tup_ssmis_ch12.png' |
---|
| 143 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 144 | |
---|
| 145 | zz=tuppgrid_tss13 |
---|
| 146 | txt='tup ssmis ch13' |
---|
| 147 | outfile='tup_ssmis_ch13.png' |
---|
| 148 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 149 | |
---|
| 150 | zz=tuppgrid_tss14 |
---|
| 151 | txt='tup ssmis ch14' |
---|
| 152 | outfile='tup_ssmis_ch14.png' |
---|
| 153 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 154 | |
---|
| 155 | zz=tuppgrid_tss15 |
---|
| 156 | txt='tup ssmis ch15' |
---|
| 157 | outfile='tup_ssmis_ch15.png' |
---|
| 158 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 159 | |
---|
| 160 | zz=tuppgrid_tss16 |
---|
| 161 | txt='tup ssmis ch16' |
---|
| 162 | outfile='tup_ssmis_ch16.png' |
---|
| 163 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 164 | |
---|
| 165 | zz=tuppgrid_tss17 |
---|
| 166 | txt='tup ssmis ch17' |
---|
| 167 | outfile='tup_ssmis_ch17.png' |
---|
| 168 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 169 | |
---|
| 170 | zz=tuppgrid_tss18 |
---|
| 171 | txt='tup ssmis ch18' |
---|
| 172 | outfile='tup_ssmis_ch18.png' |
---|
| 173 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 174 | |
---|
| 175 | ############################################# |
---|
| 176 | #amsub |
---|
| 177 | tupamsubgrid_max=np.array([tuppgrid_tab1.max(),tuppgrid_tab2.max(),tuppgrid_tab3.max(),tuppgrid_tab4.max(),tuppgrid_tab5.max()]) |
---|
| 178 | |
---|
| 179 | coef=0.6 |
---|
| 180 | cax1=0 |
---|
| 181 | cax2=300 |
---|
| 182 | cax3=1 |
---|
| 183 | my_cmap=cm.s3pcpn_l_r |
---|
| 184 | |
---|
| 185 | zz=tuppgrid_tab1 |
---|
| 186 | txt='tup amsub ch1' |
---|
| 187 | outfile='tup_amsub_ch1.png' |
---|
| 188 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 189 | |
---|
| 190 | |
---|
| 191 | zz=tuppgrid_tab2 |
---|
| 192 | txt='tup amsub ch2' |
---|
| 193 | outfile='tup_amsub_ch2.png' |
---|
| 194 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 195 | |
---|
| 196 | zz=tuppgrid_tab3 |
---|
| 197 | txt='tup amsub ch3' |
---|
| 198 | outfile='tup_amsub_ch3.png' |
---|
| 199 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 200 | |
---|
| 201 | zz=tuppgrid_tab4 |
---|
| 202 | txt='tup amsub ch4' |
---|
| 203 | outfile='tup_amsub_ch4.png' |
---|
| 204 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 205 | |
---|
| 206 | zz=tuppgrid_tab5 |
---|
| 207 | txt='tup amsub ch5' |
---|
| 208 | outfile='tup_amsub_ch5.png' |
---|
| 209 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 210 | |
---|
| 211 | |
---|
| 212 | |
---|
| 213 | ################################### |
---|
| 214 | #amsua |
---|
| 215 | |
---|
| 216 | 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()]) |
---|
| 217 | |
---|
| 218 | 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()]) |
---|
| 219 | |
---|
| 220 | zz1=tuppgrid_ta1 |
---|
| 221 | tt1='tup amsua ch1' |
---|
| 222 | t1='tup_amsua_ch1.png' |
---|
| 223 | |
---|
| 224 | zz1=tupnngrid_ta2 |
---|
| 225 | tt1='nom données tup ssmis ch15' |
---|
| 226 | t1='nom_données_tup_sssmis_ch15.png' |
---|
| 227 | |
---|
| 228 | zz1=tuppgrid_ta3 |
---|
| 229 | tt1='tup amsua ch3' |
---|
| 230 | t1='tup_amsua_ch3.png' |
---|
| 231 | |
---|
| 232 | |
---|
| 233 | zz1=tuppgrid_ta4 |
---|
| 234 | tt1='tup amsua ch4' |
---|
| 235 | t1='tup_amsua_ch4.png' |
---|
| 236 | |
---|
| 237 | import drawmap |
---|
| 238 | |
---|
| 239 | zz=tuppgrid_ta1 |
---|
| 240 | txt='tup amsua ch1' |
---|
| 241 | outfile='tup_amsua_ch11.png' |
---|
| 242 | coef=0.6 |
---|
| 243 | cax1=0 |
---|
| 244 | cax2=230 |
---|
| 245 | cax3=1 |
---|
| 246 | my_cmap=cm.s3pcpn_l_r |
---|
| 247 | |
---|
| 248 | zz=tuppgrid_ta2 |
---|
| 249 | txt='tup amsua ch2' |
---|
| 250 | outfile='tup_amsua_ch22.png' |
---|
| 251 | |
---|
| 252 | zz=tuppgrid_ta3 |
---|
| 253 | txt='tup amsua ch3' |
---|
| 254 | outfile='tup_amsua_ch33.png' |
---|
| 255 | |
---|
| 256 | zz=tuppgrid_ta4 |
---|
| 257 | txt='tup amsua ch4' |
---|
| 258 | outfile='tup_amsua_ch44.png' |
---|
| 259 | |
---|
| 260 | zz=tuppgrid_ta5 |
---|
| 261 | txt='tup amsua ch5' |
---|
| 262 | outfile='tup_amsua_ch55.png' |
---|
| 263 | |
---|
| 264 | zz=tuppgrid_ta6 |
---|
| 265 | txt='tup amsua ch6' |
---|
| 266 | outfile='tup_amsua_ch66.png' |
---|
| 267 | |
---|
| 268 | zz=tuppgrid_ta7 |
---|
| 269 | txt='tup amsua ch7' |
---|
| 270 | outfile='tup_amsua_ch77.png' |
---|
| 271 | |
---|
| 272 | zz=tuppgrid_ta15 |
---|
| 273 | txt='tup amsua ch15' |
---|
| 274 | outfile='tup_amsua_ch1515.png' |
---|
| 275 | |
---|
| 276 | drawmap.drawmap(zz, xvec, yvec, outfile, txt, cax1,cax2,cax3, my_cmap, coef) |
---|
| 277 | |
---|
| 278 | |
---|
| 279 | |
---|
| 280 | |
---|
| 281 | |
---|
| 282 | # on recupere les coordonnes xy de la projection orth standard |
---|
| 283 | m1 = Basemap(projection='ortho', lat_0 = -90, lon_0 = 0, |
---|
| 284 | resolution = 'l') |
---|
| 285 | |
---|
| 286 | xii, yii = m1(*np.meshgrid(xvec,yvec)) |
---|
| 287 | width = m1.urcrnrx - m1.llcrnrx |
---|
| 288 | height = m1.urcrnry - m1.llcrnry |
---|
| 289 | |
---|
| 290 | coef = 0.6 |
---|
| 291 | width = width*coef |
---|
| 292 | height = height*coef |
---|
| 293 | |
---|
| 294 | # on fait une nouvelle projection en zoomant sur l'antarctique |
---|
| 295 | |
---|
| 296 | map = Basemap(projection='ortho',lon_0=0,lat_0=-90,resolution='l',\ |
---|
| 297 | |
---|
| 298 | llcrnrx=-0.5*width,llcrnry=-0.5*height,urcrnrx=0.5*width,urcrnry=0.5*height) |
---|
| 299 | xii, yii = map(*np.meshgrid(xvec,yvec)) |
---|
[7] | 300 | #, clevs, cmap=my_cmap) |
---|
[6] | 301 | clevs=arange(a1,a2,a3)#star, stop, step |
---|
| 302 | cs=map.contourf(xii,yii,zz1, clevs, cmap=cm.s3pcpn_l_r) |
---|
| 303 | #cmap=cm.s3pcpn_l_r) |
---|
| 304 | #sstanom) |
---|
| 305 | #s3pcpn_l_r) |
---|
| 306 | cbar =colorbar(cs) |
---|
| 307 | plt.title(tt1) |
---|
| 308 | |
---|
| 309 | # draw coastlines, country boundaries, fill continents. |
---|
| 310 | map.drawcoastlines(linewidth=1) |
---|
| 311 | # draw the edge of the map projection region (the projection limb) |
---|
| 312 | map.drawmapboundary() |
---|
| 313 | # draw lat/lon grid lines every 30 degrees. |
---|
| 314 | map.drawmeridians(np.arange(0, 360, 30), labels=[0, 0, 0, 1]) |
---|
| 315 | map.drawparallels(np.arange(-90, 90, 10), labels=[1, 0, 0, 0]) |
---|
| 316 | #plt.show() |
---|
| 317 | plt.savefig(t1) |
---|
| 318 | close() |
---|