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

#hist tb moyenne smmis ch15

tbp_ts16=np.reshape(tbnngrid_t,-1)

bins=arange(100,300,2)
h=plt.hist(tbp_ts16, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='blue')
plt.title('histogramme tb moyenne ssmis ch16')
#plt.show()
plt.savefig('hist_tb_moyenne_ssmis_ch16.png')
close()

##########################
#hist tb moyenné ant

#mm_t matrix mask


tbp_t60=np.reshape(tbpgrid_t60,-1)
mmp_t60=np.reshape(mzzpgrid_t60,-1)

iant = (mmp_t60 >= 0.5)
tbp_t60ant=tbp_t60[iant]

bins=arange(100,300,2)
h=plt.hist(tbp_t60ant, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='blue')
plt.title('histogramme tb moyenne cont ssmis ch16')
#plt.show()
plt.savefig('hist_tb_moyenne_cont_ssmis_ch16.png')
close()

##########################
#hist tb moyenné mer et glace de mer lat 50
tbp_t50=np.reshape(tbpgrid_t,-1)
mmp_t50=np.reshape(mzzpgrid_t50,-1)

imer = (mmp_t50 < 0.5)
tbp_tmer=tbp_t50[imer]

bins=arange(100,300,2)
h=plt.hist(tbp_tmer, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='blue')
plt.title('histogramme tb moyenne mer ssmis ch16')
#plt.show()
plt.savefig('hist_tb_moyenne_mer_ssmis_ch16.png')
close()

##########################
#hist tb

bins=arange(100,300,2)
h=plt.hist(tb_t, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='blue')
plt.title('histogramme tb ssmis ch16')
plt.show()
plt.savefig('hist_tb_ssmis_ch16.png')
close()

##########################
#hist tb ant

#mm_t matrix mask

iantv= (mm_t60 >= 0.5)
tb_t60cont= tb_t60[iantv]

bins=arange(100,300,2)
h=plt.hist(tb_t60cont, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='blue')
plt.title('histogramme tb cont ssmis ch16')
#plt.show()
plt.savefig('hist_tb_cont_ssmis_ch16.png')
close()

##########################
#hist tb mer et glace de mer lat 50

imerv= (mm_t50 < 0.5)
tb_t50merv= tb_t50[imerv]

bins=arange(100,300,2)
h=plt.hist(tb_t50merv, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='blue')
plt.title('histogramme tb mer ssmis ch16')
#plt.show()
plt.savefig('hist_tb_mer_ssmis_ch16.png')
close()