;History
;2017/02/24   only data from April to October
;2016/05/02  SSA timeseries 2005- 2014 at 440nm and 368nm 
;			 To keep the consistency with monthly whisker plot, only choose data from March to October
;2015/10/20  temporily update for time series of SSA (v2)
;2014/07/08  modified because get_matchup_no2.pro was modified
;2014/01/23  quartile plot

;pro
;interquartile range(IQR)
;IQR = Q3 - Q1 : trimmed estimator (estimator derived from aother estimator by excluding some of extreme values)

inst='527'
number_of_years=12 ;2005 to 2016
time = findgen(number_of_years)+2005 ;from 2005 to 2014
max_numberdaysinyear=200
number_year = fltarr(number_of_years)

;		Filters
sza_min=30 & sza_max=75
;aot440_min = 0.40
aot440_min = 0.20
sphericity_min=95     ; minimal accepted sphericity [%]
sky_error_max = 5.0	; maximal sky error [%]
;delta_aot_max=0.05  ; use this condition for 440nm channel only
delta_aot_max=0.01  ; use this condition for 440nm channel only
mw_std_max = 0.02       ; max StDev of MFRSR retrievals of SSA  allowed
mw_error_max = 0.02
;mw_error_max = 0.05

dir='./' & name='gsfc' & title=' NASA/GSFC'
y_min=0.85 & y_max=1.0 & y_lebel=y_min-0.008

extension = '.txt'
;extension = '.it0'

!x.thick = 5.0
!y.thick = 5.0
!p.thick = 5.0
!p.charthick = 5.0

y0=0.85
mode='print' &  name_output='fig02.ps' ;name_output='yearly_boxplot_SSA_aeronet.ps'
;mode='screen'

;     Set plotting mode
if mode eq 'screen' then begin
  ;set_plot,'X' & device,pseudo_color=8
  device, decomposed=0
  window, 1, xs=1200, ys=900
  loadct, 39
endif else begin
  set_plot,'PS' & device,/color,/landscape,filename=name_output
endelse

!p.font=0
DEVICE, SET_FONT='Helvetica Bold Italic', /TT_FONT

x1 = number_of_years*2.+1.
index = indgen(number_of_years)*2

cgPlot, [0,1], /nodata, yrange=[y0,1], thick=2, xrange=[0,x1], $
  xtickformat='(A1)', ytitle='SSA', YStyle=1, XStyle=1, xticks=1 ;xtickformat='(A1)' means 'turnnig labels off'
width = ((!X.CRange[1] - !X.Crange[0]) / (number_of_years*2)) * 0.55


for h=0, number_of_years-1 do begin
iyear = time[h]
gsfc_read_matchup_allchan_adderr_v3_year, iyear, extension, inst, dir,name, delta_ssa440_max, sza_min, sza_max, aot440_min, sphericity_min, sky_error_max, delta_aot_max, mw_std_max, mw_error_max, aot_mfrsr_spec, w_mfrsr_spec, k_mfrsr_spec,  aot_aero_spec, w_aero_spec, k_aero_spec, n_aero_spec, cvf, cvc, rvf, rvc,  e_tot_aer,sf,sc,alb440,alb670,alb870,alb1020
data = w_mfrsr_spec  ;mfrsr
data2= w_aero_spec   ;AERONET
sample_size = n_elements(data[6,*])
number_year[h] = sample_size

if sample_size ge 10 then begin

c_data='black' ;AERONET
cgboxplot_nowhisker, data2[0,*],  /overplot, xlocation=index[h]+1, width=width, $
   Boxcolor=c_data, color=c_data
c_data='Red' ;MFRSR
cgboxplot_nowhisker, data[6,*], /overplot, xlocation=index[h]+1.75, width=width, $
   Boxcolor=c_data, color=c_data

print, 'AERONET'
print, index[h]+1.-width/2., index[h]+1, index[h]+1.+width/2., format='(3(F0.3, 2x))'
print, 'MFRSR'
print, index[h]+1.75-width/2., index[h]+1.75, index[h]+1.75+width/2., format='(3(F0.3, 2x))'


;if chan eq '440' then begin
;  xyouts,1.75,y0+0.03, 'AOD!d440!n>'+string(aot440_min,format='(f4.2)')+'  SZA>'+string(sza_min,format='(i2)')+'!uo!n'+'      Number = '+string(total(number_year),format='(i3)'),$
;    charthick=2,charsize=2,color=cgcolor('Blue')
;endif 
center_text = index[h]+(1.+1.75)/2.
cgText, center_text, y0+0.01, string(sample_size,format='(i3)'), Alignment=0.5, color='black', charthick=2, charsize=2 ;add labels of xaxis

endif
;stop
endfor ;h
print, number_year,total(number_year)

cgText, 1.0,  y0+0.05,'AERONET 440 nm',color=cgcolor('black'), charthick=3, charsize=2.5
cgText, 14.0, y0+0.05,'MFRSR 440 nm',color=cgcolor('red'), charthick=3, charsize=2.5
;cgText, 22.1, y0+0.06,'MFRSR '+chan+' nm',color=cgcolor('red'), charthick=2, charsize=1.5

labels = string(time, format='(I4)')
center_text = index+(1.+1.75)/2.
for j=0, number_of_years-1 do begin
  cgText, center_text[j], y0-0.01, labels[j], Alignment=0.5, charthick=2, charsize=1.3;, color=cgcolor('black')
endfor

close,/all
if mode eq 'print' then device,/close

stop
end