Pro Hist_plot_3

inpath1='E:\csr_usr1\fq\3_CMP\csr_cod\'
inpath2='E:\csr_usr1\fq\3_CMP\csr_clr\'
inpath3='E:\csr_usr1\fq\3_CMP\usr\'

;---------------SW-sum--------------------------------
psname =  'H:\Olympic\Cloud\Fig\clear\CLregion2\version2\'+'SW'+'_sum'+'.eps'
!p.multi = [0,2,2]; for six figs
tek_color
psopen,psname
Openr,oun1,path1+'SW_sum.dat',/get_lun
readf,oun1,tempsum0
Free_lun,oun1
Openr,oun2,path2+'SW_sum.dat',/get_lun
readf,oun2,tempsum1
Free_lun,oun2

 n0=mean(tempsum0[*])
 m0=stddev(tempsum0[*])
 n1=mean(tempsum1[*])
 m1=stddev(tempsum1[*])

openr,oun,path1+'SW_sum_frequency.dat',/get_lun
readf,oun,temp0
SWsum1=temp0
Free_lun,oun
print,SWsum1

openr,oun,path2+'SW_sum_frequency.dat',/get_lun
readf,oun,temp1
SWsum2=temp1
Free_lun,oun
print,SWsum2

range=[100,2,3,4,5,6,7,8,9,10]
color0=[2,2,2,2,2,2,2,2,2,2]
color1=[4,4,4,4,4,4,4,4,4,4]
xname=['120','130','140','150','160','170','180','190','200','210','220']
x = indgen(11)*1
y = indgen(10)*10
plot,x,y,/nodata,yminor = 1,xminor=1,xticks=10,$
yrange=[0,60],xtickname=xname,ytitle='Frequency (%)',xtitle='  2002-2007 clear sky  SW Flux'
xyouts,5.5,56,'Reg1:     Mean=',charsize = 0.55,color=2
xyouts,7.75,56,n0,charsize=0.55,color=2
xyouts,6.8,52,'Stddev=',charsize=0.55,color=2
xyouts,7.75,52,m0,charsize=0.55,color=2

xyouts,5.5,48,'Reg2:     Mean=',charsize = 0.55,color=4
xyouts,7.75,48,n1,charsize=0.55,color=4
xyouts,6.8,44,'Stddev=',charsize=0.55,color=4
xyouts,7.75,44,m1,charsize=0.55,color=4
wide=0.263
space=0.263
range=1.9
BAR_PLOT,SWsum1[*],/outline,/overplot,background=0,$
        barwidth=wide,baserange=range,$
         baroffset=0.0001,barspace=space,$
         colors=color0
BAR_PLOT,SWsum2[*],/outline,/overplot,background=0,$
        barwidth=wide,baserange=range,$
         baroffset=1.0,barspace=space,$
         colors=color1
print,'End'
;----------------------------------------------------------------------
;---------------LW-sum--------------------------------
Openr,oun1,path1+'LW_sum.dat',/get_lun
readf,oun1,tempsum0
Free_lun,oun1
Openr,oun2,path2+'LW_sum.dat',/get_lun
readf,oun2,tempsum1
Free_lun,oun2


 n0=mean(tempsum0[*])
 m0=stddev(tempsum0[*])
 n1=mean(tempsum1[*])
 m1=stddev(tempsum1[*])

openr,oun,path1+'LW_sum_frequency.dat',/get_lun
readf,oun,temp0
LWsum1=temp0
Free_lun,oun
print,LWsum1

openr,oun,path2+'LW_sum_frequency.dat',/get_lun
readf,oun,temp1
LWsum2=temp1
Free_lun,oun
print,LWsum2

range=[100,2,3,4,5,6,7,8,9,10]
color0=[2,2,2,2,2,2,2,2,2,2]
color1=[4,4,4,4,4,4,4,4,4,4]
xname=['250','260','270','280','290','300','310','320','330','340','350']
x = indgen(11)*1
y = indgen(10)*10
plot,x,y,/nodata,yminor = 1,xminor=1,xticks=10,$
yrange=[0,60],xtickname=xname,ytitle='Frequency (%)',xtitle='  2002-2007 clear sky  LW Flux'
xyouts,5.5,56,'Reg1:     Mean=',charsize = 0.55,color=2
xyouts,7.75,56,n0,charsize=0.55,color=2
xyouts,6.8,52,'Stddev=',charsize=0.55,color=2
xyouts,7.75,52,m0,charsize=0.55,color=2

xyouts,5.5,48,'Reg2:     Mean=',charsize = 0.55,color=4
xyouts,7.75,48,n1,charsize=0.55,color=4
xyouts,6.8,44,'Stddev=',charsize=0.55,color=4
xyouts,7.75,44,m1,charsize=0.55,color=4
wide=0.263
space=0.263
range=1.9
BAR_PLOT,LWsum1[*],/outline,/overplot,background=0,$
        barwidth=wide,baserange=range,$
         baroffset=0.0001,barspace=space,$
         colors=color0
BAR_PLOT,LWsum2[*],/outline,/overplot,background=0,$
        barwidth=wide,baserange=range,$
         baroffset=1.0,barspace=space,$
         colors=color1
print,'End'
psclose
;----------------------------------------------------------------------

;---------------Tau1-sum--------------------------------
psname =  'H:\Olympic\Cloud\Fig\clear\CLregion2\version2\'+'TAU'+'_sum'+'.eps'
!p.multi = [0,2,2]; for six figs
tek_color
psopen,psname
Openr,oun1,path1+'Tau1_sum.dat',/get_lun
readf,oun1,tempsum2
Free_lun,oun1
Openr,oun2,path2+'Tau1_sum.dat',/get_lun
readf,oun2,tempsum3
Free_lun,oun2


 n0=mean(tempsum2[*])
 m0=stddev(tempsum2[*])
 n1=mean(tempsum3[*])
 m1=stddev(tempsum3[*])

openr,oun,path1+'Tau1_sum_frequency.dat',/get_lun
readf,oun,temp0
Tau1sum1=temp0
Free_lun,oun
print,Tau1sum1

openr,oun,path2+'Tau1_sum_frequency.dat',/get_lun
readf,oun,temp1
Tau1sum2=temp1
Free_lun,oun
print,Tau1sum2


range=[100,2,3,4,5,6,7,8,9,10]
color0=[2,2,2,2,2,2,2,2,2,2]
color1=[4,4,4,4,4,4,4,4,4,4]
xname=['0','0.25','0.5','0.75','1.0','1.25','1.5','1.75','2.0','2.25','2.5']
x = indgen(11)*1
y = indgen(10)*10
plot,x,y,/nodata,yminor = 1,xminor=1,xticks=10,$
yrange=[0,80],xtickname=xname,ytitle='Frequency (%)',xtitle='  2002-2007 clear sky  aerosol OPD (0.47mm)'
xyouts,5.5,75,'Reg1:      Mean=',charsize = 0.75,color=2
xyouts,7.75,75,n0,charsize=0.75,color=2
xyouts,6.8,70,'Stddev=',charsize=0.75,color=2
xyouts,7.75,70,m0,charsize=0.75,color=2

xyouts,5.5,65,'Reg2:      Mean=',charsize = 0.75,color=4
xyouts,7.75,65,n1,charsize=0.75,color=4
xyouts,6.8,60,'Stddev=',charsize=0.75,color=4
xyouts,7.75,60,m1,charsize=0.75,color=4
wide=0.263
space=0.263
range=1.9
BAR_PLOT,Tau1sum1[*],/outline,/overplot,background=0,$
        barwidth=wide,baserange=range,$
         baroffset=0.0001,barspace=space,$
         colors=color0
BAR_PLOT,Tau1sum2[*],/outline,/overplot,background=0,$
        barwidth=wide,baserange=range,$
         baroffset=1.0,barspace=space,$
         colors=color1

print,'End'
;----------------------------------------------------------------------

;---------------Tau2-sum--------------------------------
Openr,oun1,path1+'Tau2_sum.dat',/get_lun
readf,oun1,tempsum2
Free_lun,oun1
Openr,oun2,path2+'Tau2_sum.dat',/get_lun
readf,oun2,tempsum3
Free_lun,oun2


 n0=mean(tempsum2[*])
 m0=stddev(tempsum2[*])
 n1=mean(tempsum3[*])
 m1=stddev(tempsum3[*])

openr,oun,path1+'Tau2_sum_frequency.dat',/get_lun
readf,oun,temp0
Tau2sum1=temp0
Free_lun,oun
print,Tau1sum2

openr,oun,path2+'Tau2_sum_frequency.dat',/get_lun
readf,oun,temp1
Tau2sum2=temp1
Free_lun,oun
print,Tau1sum2


range=[100,2,3,4,5,6,7,8,9,10]
color0=[2,2,2,2,2,2,2,2,2,2]
color1=[4,4,4,4,4,4,4,4,4,4]
xname=['0','0.25','0.5','0.75','1.0','1.25','1.5','1.75','2.0','2.25','2.5']
x = indgen(11)*1
y = indgen(10)*10
plot,x,y,/nodata,yminor = 1,xminor=1,xticks=10,$
yrange=[0,80],xtickname=xname,ytitle='Frequency (%)',xtitle='  2002-2007 clear sky  aerosol OPD (0.55mm)'
xyouts,5.5,75,'Reg1:      Mean=',charsize = 0.75,color=2
xyouts,7.75,75,n0,charsize=0.75,color=2
xyouts,6.8,70,'Stddev=',charsize=0.75,color=2
xyouts,7.75,70,m0,charsize=0.75,color=2

xyouts,5.5,65,'Reg2:      Mean=',charsize = 0.75,color=4
xyouts,7.75,65,n1,charsize=0.75,color=4
xyouts,6.8,60,'Stddev=',charsize=0.75,color=4
xyouts,7.75,60,m1,charsize=0.75,color=4
wide=0.263
space=0.263
range=1.9
BAR_PLOT,Tau2sum1[*],/outline,/overplot,background=0,$
        barwidth=wide,baserange=range,$
         baroffset=0.0001,barspace=space,$
         colors=color0
BAR_PLOT,Tau2sum2[*],/outline,/overplot,background=0,$
        barwidth=wide,baserange=range,$
         baroffset=1.0,barspace=space,$
         colors=color1

print,'End'
;----------------------------------------------------------------------

;---------------Tau2-sum--------------------------------
Openr,oun1,path1+'Tau3_sum.dat',/get_lun
readf,oun1,tempsum2
Free_lun,oun1
Openr,oun2,path2+'Tau3_sum.dat',/get_lun
readf,oun2,tempsum3
Free_lun,oun2


 n0=mean(tempsum2[*])
 m0=stddev(tempsum2[*])
 n1=mean(tempsum3[*])
 m1=stddev(tempsum3[*])

openr,oun,path1+'Tau3_sum_frequency.dat',/get_lun
readf,oun,temp0
Tau3sum1=temp0
Free_lun,oun
print,Tau1sum2

openr,oun,path2+'Tau3_sum_frequency.dat',/get_lun
readf,oun,temp1
Tau3sum2=temp1
Free_lun,oun
print,Tau1sum2


range=[100,2,3,4,5,6,7,8,9,10]
color0=[2,2,2,2,2,2,2,2,2,2]
color1=[4,4,4,4,4,4,4,4,4,4]
xname=['0','0.25','0.5','0.75','1.0','1.25','1.5','1.75','2.0','2.25','2.5']
x = indgen(11)*1
y = indgen(10)*10
plot,x,y,/nodata,yminor = 1,xminor=1,xticks=10,$
yrange=[0,80],xtickname=xname,ytitle='Frequency (%)',xtitle='  2002-2007 clear sky  aerosol OPD (0.659mm)'
xyouts,5.5,75,'Reg1:      Mean=',charsize = 0.75,color=2
xyouts,7.75,75,n0,charsize=0.75,color=2
xyouts,6.8,70,'Stddev=',charsize=0.75,color=2
xyouts,7.75,70,m0,charsize=0.75,color=2

xyouts,5.5,65,'Reg2:      Mean=',charsize = 0.75,color=4
xyouts,7.75,65,n1,charsize=0.75,color=4
xyouts,6.8,60,'Stddev=',charsize=0.75,color=4
xyouts,7.75,60,m1,charsize=0.75,color=4
wide=0.263
space=0.263
range=1.9
BAR_PLOT,Tau3sum1[*],/outline,/overplot,background=0,$
        barwidth=wide,baserange=range,$
         baroffset=0.0001,barspace=space,$
         colors=color0
BAR_PLOT,Tau3sum2[*],/outline,/overplot,background=0,$
        barwidth=wide,baserange=range,$
         baroffset=1.0,barspace=space,$
         colors=color1

print,'End'
;----------------------------------------------------------------------

psclose
END
;---------------for 4 Figs-------------------
PRO psopen, psfname
set_plot,'ps'
device,filename=psfname
device,/landscape
device, xoffset=2., yoffset=11.5,/inch
device, xsize=11.0, ysize=8.0,/inch
device,/Helvetica,/bold
device,/encapsulated
 device,/color
!p.font = 0
 !p.thick = 6
 !x.thick = 6
 !y.thick = 6
 !p.charsize=0.8
End
;-------------------------------------
;---------------------------------------------
;---------------for 6 Figs-------------------
;Pro psopen,psname
; set_plot,'ps'
; device, /encapsulated,preview=2
; device,filename = psname
; device,/portrait
; device,/helvetica,font_size=20,/bold
; device, /color, bits_per_pixel=8  ; 2**8=256 colors
; device, xsize=8, ysize=11.2,/inch
; device,xoff=3,yoff=3
; !p.font = 0
; !p.thick = 6
; !x.thick = 6
; !y.thick = 6
; !p.charsize=1.0
;End
;-------------------------------------

Pro psclose
 device,/close
 !p.thick = 1
 !x.thick = 1
 !y.thick = 1

End