pro plot_cns_together,stats_directory=stats_directory 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