pro plot_cns_together,stats_directory=stats_directory
;what
If (not keyword_set(stats_directory)) then stats_directory='stats'
Date=['10','11','12']

cxl=fltarr(3,25)
cxt=fltarr(3,25)
cyt=fltarr(3,25)
cyl=fltarr(3,25)

for i=0,2 do begin

stats_directory='/a/pippin01/Volumes/u08/lsst/CPanalysis/2005-05-'+Date(i)+$
          '/ShackHartman/stats/'
postscript_Dir='/a/pippin01/Volumes/u08/lsst/CPanalysis/2005-05-'+Date(0)+$
          '/ShackHartman/distorted_postscripts/'

readcol,stats_directory+'Depth_0Correlations.txt',cxl_temp,cxt_temp,$
	cyl_temp,cyt_temp,format='(f,f,f,f)'

cxl(i,*)=cxl_temp
cxt(i,*)=cxt_temp
cyl(i,*)=cyl_temp
cyt(i,*)=cyt_temp

endfor

set_plot,'ps'
device,filename=Postscript_dir+'All_Correlations.ps'

loadct,39                       ;load color table 39
device,/color                   ;allow color on the postscript
device,ysize=8.5,/inches        ;Height of plot in y
device,xsize=10.0,/inches        ;Width of plot in x
device,yoffset=1.0,/inches      ;Y position of lower left corner

white='FFFFFF'x
black='000000'x
!P.CHARSIZE=1.2
!P.THICK=4.

items=['10','11','12']
colors=[0,100,200]
sym_arr=[2,2,2]

plot,[0,24],[0,1.5],$
     /xstyle,/ystyle,/nodata,background='FFFFFF'x,color='000000'x,$
     title='X Longitudinal',ytitle='Covariance',xtitle='n'

for i=0,2 do begin
    oplot,cxl(i,*),color=100*i,psym=2
endfor

legend,items,psym=sym_arr,color=colors

plot,[0,24],[0,1.5],$
     /xstyle,/ystyle,/nodata,background='FFFFFF'x,color='000000'x,$
     title='Y Transverse',ytitle='Covariance',xtitle='n'

for i=0,2 do begin
    oplot,cyt(i,*),color=100*i,psym=2
endfor

legend,items,psym=sym_arr,color=colors

plot,[0,24],[0,1.5],$
     /xstyle,/ystyle,/nodata,background='FFFFFF'x,color='000000'x,$
     title='Y Longitudinal',ytitle='Covariance',xtitle='n'

for i=0,2 do begin
    oplot,cyl(i,*),color=100*i,psym=2
endfor

legend,items,psym=sym_arr,color=colors

plot,[0,24],[0,1.5],$
     /xstyle,/ystyle,/nodata,background='FFFFFF'x,color='000000'x,$
     title='X Transverse',ytitle='Covariance',xtitle='n'

for i=0,2 do begin
    oplot,cxt(i,*),color=100*i,psym=2
endfor

legend,items,psym=sym_arr,color=colors

plot,[0,24],[0,1.5],$
     /xstyle,/ystyle,/nodata,background='FFFFFF'x,color='000000'x,$
     title='Night of the 10th',ytitle='Covariance',xtitle='n'

oplot,cxl(0,*),color=0,psym=2
oplot,cyl(0,*),color=70,psym=2
oplot,cyt(0,*),color=140,psym=2
oplot,cxt(0,*),color=210,psym=2

colors=[0,70,140,210]
sym_arr=[2,2,2,2]
items=['X Longitudinal','Y Longitudinal','Y Transverse','X Transverse']

legend,items,color=colors,psym=sym_arr

plot,[0,24],[0,1.5],$
     /xstyle,/ystyle,/nodata,background='FFFFFF'x,color='000000'x,$
     title='Night of the 11th',ytitle='Covariance',xtitle='n'

oplot,cxl(1,*),color=0,psym=2
oplot,cyl(1,*),color=70,psym=2
oplot,cyt(1,*),color=140,psym=2
oplot,cxt(1,*),color=210,psym=2

colors=[0,70,140,210]
sym_arr=[2,2,2,2]
items=['X Longitudinal','Y Longitudinal','Y Transverse','X Transverse']

legend,items,color=colors,psym=sym_arr

plot,[0,24],[0,1.5],$
     /xstyle,/ystyle,/nodata,background='FFFFFF'x,color='000000'x,$
     title='Night of the 12th',ytitle='Covariance',xtitle='n'

oplot,cxl(2,*),color=0,psym=2
oplot,cyl(2,*),color=70,psym=2
oplot,cyt(2,*),color=140,psym=2
oplot,cxt(2,*),color=210,psym=2

colors=[0,70,140,210]
sym_arr=[2,2,2,2]
items=['X Longitudinal','Y Longitudinal','Y Transverse','X Transverse']

legend,items,color=colors,psym=sym_arr

device,/close
set_plot,'X'

end