;save,lat_data,lon_data,cloud_fr,cloud_per, succ_rn,succ_rn60, ncl60land,nclTland,tot_per,per_60,$
;ri_avg,ri_std,df_avg,df_std,mean_coe,std_coe,mean_co,std_co,$
;ri_avg60,ri_std60,df_avg60,df_std60,mean_coe60,std_coe60,mean_co60,std_co60,$
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;----------------------------------------------------
;modified May 2020 for the paper
@cgLayout
;------------------------------------------------------------------------------------------
FILENAME1='cloudInfoAug23.sav'
restore, filename1 ;to get the 0.95 case which is 2 index

default_succ_rn=succ_rn(2)
default_succ_rn60=succ_rn60(2)

default_tot_per=tot_per(2)
default_per_60=per_60(2)

default_cloud_per=cloud_per(*,2)
default_cloud_fr=cloud_fr(*,2)

;save,cloud_fr,coud_per, succ_rn,succ_rn60, ncl60day,nclTday,tot_per,per_60,$
;ri_avg,ri_std,df_avg,df_std,mean_coe,std_coe,mean_co,std_co,FILENAME='cloudInfoAug23.sav'

;FILENAME='cloudInfoAug23MopRad_NIR.sav'
FILENAME='cloudInfoAug23MopRad.sav'
restore, filename
test_succ_rn=succ_rn
test_succ_rn60=succ_rn60

test_tot_per=tot_per
test_per_60=per_60
test_cloud_per=cloud_per
test_cloud_fr=cloud_fr
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

;we need to get rid off some tests to show the varations 

;!P.Multi = [0, 1, 2]
;Test=['0.01-modis-99.9','0.01','0.1','0.5','0.5-modis-99.9','0.7','0.8','0.8-1.2','0.9','0.95','1.0','1.2']

Test=['0.01','0.1','0.5','0.7','0.8','0.9','0.95','1.0']

nt=12
ntt=8
help,Test
succ_rn=Lonarr(nt) ;successfull reterieval number
succ_rn60=Lonarr(nt)


new_succ_rn=Lonarr(8)
new_succ_rn60=Lonarr(8)
new_tot_per=fltarr(8)
new_per_60=fltarr(8)
new_cloud_per=fltarr(6,8)
new_cloud_fr=fltarr(6,8)

tot_per=fltarr(nt)
per_60=fltarr(nt)
cloud_per=fltarr(6,nt) ;flag number and ntests
cloud_fr=fltarr(6,nt) ;flag number and ntests

succ_rn(0:8)=test_succ_rn(0:8)
succ_rn(9)=default_succ_rn
succ_rn(10:11)=test_succ_rn(9:10)

new_succ_rn(0:2)=succ_rn(1:3)
new_succ_rn(3:4)=succ_rn(5:6)
new_succ_rn(5:7)=succ_rn(8:10)

succ_rn60(0:8)=test_succ_rn60(0:8)
succ_rn60(9)=default_succ_rn60
succ_rn60(10:11)=test_succ_rn60(9:10)

new_succ_rn60(0:2)=succ_rn60(1:3)
new_succ_rn60(3:4)=succ_rn60(5:6)
new_succ_rn60(5:7)=succ_rn60(8:10)


tot_per(0:8)=test_tot_per(0:8)
tot_per(9)=default_tot_per
tot_per(10:11)=test_tot_per(9:10)

new_tot_per(0:2)=tot_per(1:3)
new_tot_per(3:4)=tot_per(5:6)
new_tot_per(5:7)=tot_per(8:10)

per_60(0:8)=test_per_60(0:8)
per_60(9)=default_per_60
per_60(10:11)=test_per_60(9:10)

new_per_60(0:2)=per_60(1:3)
new_per_60(3:4)=per_60(5:6)
new_per_60(5:7)=per_60(8:10)

cloud_per(*,0:8)=test_cloud_per(*,0:8)
cloud_per(*,9)=default_cloud_per
cloud_per(*,10:11)=test_cloud_per(*,9:10)

new_cloud_per(*,0:2)=cloud_per(*,1:3)
new_cloud_per(*,3:4)=cloud_per(*,5:6)
new_cloud_per(*,5:7)=cloud_per(*,8:10)

cloud_fr(*,0:8)=test_cloud_fr(*,0:8)
cloud_fr(*,9)=default_cloud_fr
cloud_fr(*,10:11)=test_cloud_fr(*,9:10)

new_cloud_fr(*,0:2)=cloud_fr(*,1:3)
new_cloud_fr(*,3:4)=cloud_fr(*,5:6)
new_cloud_fr(*,5:7)=cloud_fr(*,8:10)
;;;;;;;;;;;;;;;


loadct,39
SET_PLOT, 'PS'

pos=cgLayout([2, 2], OXMargin=[10, 10], OYMargin=[8, 4], XGap=10, YGap=4)
;ck=nt-1 ;number of ticks
ck=ntt-1 ;number of ticks
z=1

epsfile='plot_Aug23TestMopRad_v2.eps'
  print,epsfile
 device, file =epsfile,/bold,/color,/encapsulated,xsize=18,xoffset=0.2,/PORTRAIT,ysize=12

ytitle='MOPITT pixel number'
xtitle='MOPITT rad. ratio'
position=pos[*,0]
cgplot,new_succ_rn,xtickname=Test,xticks=ck,charsize=z-0.2,/nodata,xtitle=xtitle,ytitle=ytitle,font=0,xstyle=1,/noerase,$
yrange=[80000,700000],ystyle=1,xticklen=1,Position=position,XTickFormat='(A1)';[0.1, 0.55,0.9,0.95]

 cgoplot,new_succ_rn60,psym=-16,symsize=z-0.4,SymColor='blu7',font=0
 cgoplot,new_succ_rn,psym=-16,symsize=z-0.4,SymColor='olive',font=0
    
AL_Legend, ['Global number of successful pixels','Pixels from 60N-60S'] ,$
PSym=[16,16], Position=[0.2, 0.9], /Normal, Color=['olive','blu7'],charsize=z-0.5,font=0,symsize=[z-0.5,z-0.5]
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

ytitle='MOPITT pixel percent'
xtitle=' '
position=pos[*,2]
cgplot,new_tot_per,xtickname=Test,xticks=ck,charsize=z-0.2,/nodata,xtitle=xtitle,ytitle=ytitle,font=0,xstyle=1,/noerase,$
yrange=[15,100],ystyle=1,xticklen=1,Position=position,XTickFormat='(A1)';[0.1, 0.1,0.9,0.5]

 cgoplot,new_per_60,psym=-16,symsize=z-0.5,SymColor='blu7',font=0
 cgoplot,new_tot_per,psym=-16,symsize=z-0.5,SymColor='olive',font=0
ypos = Replicate(!Y.Window[0] - 0.04, ck+1)
   xpos = !X.Window[0] + (!X.Window[1] - !X.Window[0]) * $
      Findgen(ck+ 1) / ck

FOR j=0, ck DO XYOutS, xpos[j], ypos[j], $
      Test[j], Alignment=0.0, Orientation=-45,charsize=z-0.5,font=0, /Normal

AL_Legend, ['Global percent of successful pixels','Pixel percent from 60N-60S'] ,$
PSym=[16,16], Position=[0.2, 0.48], /Normal, Color=['olive','blu7'],charsize=z-0.5,font=0,symsize=[z-0.5,z-0.5]

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
ytitle='Pixel number for each flag'
flag=findgen(6)+1
position=pos[*,1]
cgplot,flag,new_cloud_fr(*,0),xticks=5,charsize=z-0.2,/nodata,xtitle='Flag Num.',ytitle=ytitle,font=0,xstyle=1,/noerase,$
yrange=[0,400000],ystyle=1,xticklen=1,Position=position;[0.1, 0.1,0.9,0.5]
 ;colors=['gold','blue','dark green',

colors = ['cornflower blue','blue','sea green','green','gold','goldenrod','maroon','hot pink','plum','dark orchid','charcoal','black']
  for f=0, ntt-1 do begin
 cgoplot,flag,new_cloud_fr(*,f),psym=-16,symsize=z-0.5,SymColor=colors(f),font=0
  endfor
help,Test
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

ytitle='Pixel percent for each flag'
flag=findgen(6)+1
position=pos[*,3]
cgplot,flag,new_cloud_per(*,0),xticks=5,charsize=z-0.2,/nodata,xtitle='Flag Num.',ytitle=ytitle,font=0,xstyle=1,/noerase,$
yrange=[0,80],ystyle=1,xticklen=1,Position=position;[0.1, 0.1,0.9,0.5]
 ;colors=['gold','blue','dark green',

  for f=0, ntt-1 do begin
 cgoplot,flag,new_cloud_per(*,f),psym=-16,symsize=z-0.5,SymColor=colors(f),font=0
  endfor
help,Test
pp=intarr(ntt)+16
ss=fltarr(ntt)+z-0.5
AL_Legend, Test ,$
PSym=pp, Position=[0.87, 0.65], /Normal,$
 Color=colors,charsize=z-0.5,font=0,symsize=ss
!P.Multi = 0
device,/close
print,epsfile
end