source: trunk/src/python_script/hist tup.py @ 6

############
#histogrammes moyennées
#tup ssmis ch15, ch16 et amsua ch2


tupp_ta260=np.reshape(tuppgrid_ta260,-1)
mmp_ta260=np.reshape(mpgrid_t60,-1)

iant = (mmp_ta260 >= 0.5)
tupp_ta260ant=tupp_ta260[iant]

bins=arange(2,14,0.2)
h=plt.hist(tupp_ta260ant, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='blue')
plt.title('histogramme tup moyenne cont amsua ch2')
#plt.show()
plt.savefig('hist_tup_moyenne_cont_amsua_ch2.png')
close()

#ssmis

tupp_t1660=np.reshape(tuppgrid_t1660,-1)
tupp_t1660ant=tupp_t1660[iant]

bins=arange(7,24,0.2)
h=plt.hist(tupp_t1660ant, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='blue')
plt.title('histogramme tup moyenne cont ssmis ch16')
#plt.show()
plt.savefig('hist_tup_moyenne_cont_ssmis_ch16.png')
close()

########################
#hist pas moyennées
#amsua ch2

iant60v = (mm_t60 >= 0.5)
tup_ta2ant=tup_ta260[iant60v]

bins=arange(0,24,0.2)
h=plt.hist(tup_ta2ant, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='blue')
plt.title('histogramme tup cont amsua ch2')
#plt.show()
plt.savefig('hist_tup_cont_amsua_ch2.png')
close()

#ssmis ch16

iant60v = (mm_t1660 >= 0.5)
tup_t16ant=tup_t1660[iant60v]

bins=arange(7,30,0.2)
h=plt.hist(tup_t16ant, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='blue')
plt.title('histogramme tup cont ssmis ch16')
#plt.show()
plt.savefig('hist_tup_cont_ssmis_ch16.png')
close()