source: trunk/src/python_script/hist_ssmis_ch15.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: 2.1 KB
Line 
1#hist tb moyenne smmis ch15
2
3tbp_ts15=np.reshape(tbnngrid_t,-1)
4
5bins=arange(100,300,2)
6h=plt.hist(tbp_ts15, bins=bins,histtype='bar',rwidth=0.2,
7label='hist',normed='True',color='blue')
8plt.title('histogramme tb moyenne ssmis ch15')
9#plt.show()
10plt.savefig('hist_tb_moyenne_ssmis_ch15.png')
11close()
12
13##########################
14#hist tb moyenné ant
15
16#mm_t matrix mask
17
18mzzpgrid_lat60=mzzpgrid_t
19mmvec_lat60=mm_t
20
21tbp_t60=np.reshape(tbpgrid_t60,-1)
22mmp_t60=np.reshape(mzzpgrid_lat60,-1)
23
24iant = (mmp_t60 >= 0.5)
25tbp_t60ant=tbp_t60[iant]
26
27bins=arange(100,300,2)
28h=plt.hist(tbp_t60ant, bins=bins,histtype='bar',rwidth=0.2,
29label='hist',normed='True',color='blue')
30plt.title('histogramme tb moyenne cont ssmis ch15')
31#plt.show()
32plt.savefig('hist_tb_moyenne_cont_ssmis_ch15.png')
33close()
34
35##########################
36#hist tb moyenné mer et glace de mer lat 50
37tbp_t50=np.reshape(tbpgrid_t,-1)
38mmp_t50=np.reshape(mzzpgrid_t50,-1)
39
40imer = (mmp_t50 < 0.5)
41tbp_tmer=tbp_t50[imer]
42
43bins=arange(100,300,2)
44h=plt.hist(tbp_tmer, bins=bins,histtype='bar',rwidth=0.2,
45label='hist',normed='True',color='blue')
46plt.title('histogramme tb moyenne mer ssmis ch15')
47#plt.show()
48plt.savefig('hist_tb_moyenne_mer_ssmis_ch15.png')
49close()
50
51##########################""
52#hist tb moyenne smmis ch15
53
54bins=arange(100,300,2)
55h=plt.hist(tb_t, bins=bins,histtype='bar',rwidth=0.2,
56label='hist',normed='True',color='blue')
57plt.title('histogramme tb ssmis ch15')
58#plt.show()
59plt.savefig('hist_tb_ssmis_ch15.png')
60close()
61
62##########################
63#hist tb ant
64
65#mm_t matrix mask
66
67iantv= (mmvec_lat60 >= 0.5)
68tb_t60cont= tb_t60[iantv]
69
70bins=arange(100,300,2)
71h=plt.hist(tb_t60cont, bins=bins,histtype='bar',rwidth=0.2,
72label='hist',normed='True',color='blue')
73plt.title('histogramme tb cont ssmis ch15')
74#plt.show()
75plt.savefig('hist_tb_cont_ssmis_ch15.png')
76close()
77
78##########################
79#hist tb mer et glace de mer lat 50
80
81imerv= (mm_t < 0.5)
82tb_t50merv= tb_t[imerv]
83
84bins=arange(100,300,2)
85h=plt.hist(tb_t50merv, bins=bins,histtype='bar',rwidth=0.2,
86label='hist',normed='True',color='blue')
87plt.title('histogramme tb mer ssmis ch15')
88#plt.show()
89plt.savefig('hist_tb_mer_ssmis_ch15.png')
90close()
Note: See TracBrowser for help on using the repository browser.