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

#hist tb moyenne

tbp_t=np.reshape(tbpgrid_t,-1)

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

##########################
#mer
mmp_t=np.reshape(mzzpgrid_t,-1)
imer = (mmp_t < 0.5)
tbp_tmer=tbp_t[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 amsua ch2 new')
#plt.show()
plt.savefig('hist_tb_moyenne_mer_amsua_ch2_new.png')
close()
#################################################
#contienente pas moyenne
iant60 = (mm_t60 >= 0.5)
tb_tant=tb_t60[iant60]

bins=arange(100,300,2)
h=plt.hist(tb_tant, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='blue')
plt.title('histogramme tb cont amsua ch2 new')
plt.show()
plt.savefig('hist_tb_cont_amsua_ch2_new.png')
close()

#solo alt>

alt_tant=alt_t60[iant60]
ialt1= (alt_tant>=2000)
tb_talt1=tb_tant[ialt1]

ialt2= (alt_tant<2000)
tb_talt2=tb_tant[ialt2]


fig=figure(1)
bins=arange(100,300,2)
h1=plt.hist(tb_talt1, bins=bins,histtype='bar',rwidth=0.4,
label='>2000m',normed='True',color='red')
h2=plt.hist(tb_talt2, bins=bins,histtype='bar',rwidth=0.1,
label='<2000m',normed='True',color='blue')
plt.legend()
fig.suptitle('Histogrammes tb alt 2000m AMSUA CH2', fontsize=14)
plt.savefig('hist_tb_alt2000m_amsua_ch2.png')
close()


alt_tant=alt_t60[iant60]
ialt11= (alt_tant>=4000)
tb_talt11=tb_tant[ialt11]

ialt22= (alt_tant>=3000) & (alt_tant<4000)
tb_talt22=tb_tant[ialt22]

ialt33= (alt_tant>=2000) & (alt_tant<3000)
tb_talt33=tb_tant[ialt33]

ialt44= (alt_tant>=1000) & (alt_tant<2000)
tb_talt44=tb_tant[ialt44]

ialt55= (alt_tant<1000)
tb_talt55=tb_tant[ialt55]


fig=figure(1)
bins=arange(100,300,2)
h1=plt.hist(tb_talt11, bins=bins,histtype='step',rwidth=0.4,
label='>400m',normed='True',color='red')
h2=plt.hist(tb_talt22, bins=bins,histtype='step',rwidth=0.1,
label='400m - 3000m',normed='True',color='blue')
h3=plt.hist(tb_talt33, bins=bins,histtype='step',rwidth=0.1,
label='3000m - 2000m',normed='True',color='green')
h4=plt.hist(tb_talt44, bins=bins,histtype='step',rwidth=0.1,
label='2000m - 1000m',normed='True',color='orange')
h5=plt.hist(tb_talt55, bins=bins,histtype='step',rwidth=0.1,
label='>1000m',normed='True',color='black')
fig.suptitle('Histogrammes tb AMSUA CH2', fontsize=14)
plt.legend()
#plt.show()
plt.savefig('hist_tb_dif_alt_amsua_ch2.png')
close()


##############

#hist tb continent
tbp_tc=np.reshape(tbpgrid_tc,-1)

bins=arange(100,300,2)
h=plt.hist(tbp_tc, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb moyenne cont')
plt.show()
plt.savefig('hist_tb_moyenne_cont_amsua_ch2.png')
close()

#hist tb galce mer

tbp_mer=np.reshape(mertb,-1)
bins=arange(100,300,2)
h=plt.hist(tbp_mer, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb moyenne mer galce de mer')
#plt.show()
plt.savefig('hist_tb_moyenne_mer_glace_de_mer_amsua_ch2.png')
close()

#ឧist tb cont pos1 

tbp_1nc=np.reshape(tbpgrid_1nc,-1)
bins=arange(100,300,2)
h=plt.hist(tbp_1nc, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb moyenne cont pos1')
#plt.show()
plt.savefig('hist_tb_moyenne_cont_pos1_amsua_ch2.png')
close()

#ឧist tb cont pos2 

tbp_2nc=np.reshape(tbpgrid_2nc,-1)
bins=arange(100,300,2)
h=plt.hist(tbp_2nc, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb moyenne cont pos2')
#plt.show()
plt.savefig('hist_tb_moyenne_cont_pos2_amsua_ch2.png')
close()

#ឧist tb cont pos3 

tbp_3nc=np.reshape(tbpgrid_3nc,-1)
bins=arange(100,300,2)
h=plt.hist(tbp_3nc, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb moyenne cont pos3')
#plt.show()
plt.savefig('hist_tb_moyenne_cont_pos3_amsua_ch2.png')
close()


################################################################

mertb_1n

#ឧist tb mer pos1 

tbmer_1n=np.reshape(mertb_1n,-1)
bins=arange(100,300,2)
h=plt.hist(tbmer_1n, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb moyenne mer pos1')
#plt.show()
plt.savefig('hist_tb_moyenne_mer_pos1_amsua_ch2.png')
close()

#ឧist tb mer pos2

tbmer_2n=np.reshape(mertb_2n,-1)
bins=arange(100,300,2)
h=plt.hist(tbmer_2n, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb moyenne mer pos2')
#plt.show()
plt.savefig('hist_tb_moyenne_mer_pos2_amsua_ch2.png')
close()

#ឧist tb mer pos3

tbmer_3n=np.reshape(mertb_3n,-1)
bins=arange(100,300,2)
h=plt.hist(tbmer_3n, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb moyenne mer pos3')
#plt.show()
plt.savefig('hist_tb_moyenne_mer_pos3_amsua_ch2.png')
close()

#######################################################""
#hist tb pas moyenné cont

iterre = (mm_t >= 0.5)
tb_tcnm=tb_t[iterre]



bins=arange(100,300,2)
h=plt.hist(tb_tcnm, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb ant')
#plt.show()
plt.savefig('hist_tb_ant_amsua_ch2.png')
close()



#hist tb mer pas moyenné 
imer =  (mm_t50 < 0.5)
tb_tmer=tb_t50[imer]

bins=arange(100,300,2)
h=plt.hist(tb_tmer, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb mer')
plt.show()
plt.savefig('hist_tb_mer_amsua_ch2.png')
close()



#glace de mer et plataformes
iglace = tb_tmer>=190
tb_tglace= tb_tmer[iglace]

bins=arange(100,300,2)
h=plt.hist(tb_tglace, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb glace de mer (tb>190)')
#plt.show()
plt.savefig('hist_tb_glace_de_mer190_amsua_ch2.png')
close()




#hist tb ant pas moyenné  pos1
it_1n = (mm_tl1n >= 0.5)
tb_1naa=tb_1na[it_1n]

bins=arange(100,300,2)
h=plt.hist(tb_1naa, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb ant pos1')
#plt.show()
plt.savefig('hist_tb_ant_pos1_amsua_ch2.png')
close()

#hist tb ant pas moyenné  pos2
it_2n = (mm_tl2n >= 0.5)
tb_2naa=tb_2na[it_2n]

bins=arange(100,300,2)
h=plt.hist(tb_2naa, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb ant pos2')
#plt.show()
plt.savefig('hist_tb_ant_pos2_amsua_ch2.png')
close()

#hist tb ant pas moyenné  pos3
it_3n = (mm_tl3n >= 0.5)
tb_3naa=tb_3na[it_3n]

bins=arange(100,300,2)
h=plt.hist(tb_3naa, bins=bins,histtype='bar',rwidth=0.2,
label='hist',normed='True',color='green')
plt.title('histogramme tb ant pos3')
#plt.show()
plt.savefig('hist_tb_ant_pos3_amsua_ch2.png')
close()