pro make_correlation_jpegs If (not keyword_set(Date)) then Date='10' If (not keyword_set(stats_directory)) then stats_directory='stats_temp' stats_directory='/a/pippin01/Volumes/u08/lsst/CPanalysis/2005-05-'+Date+$ '/ShackHartman/'+stats_directory+'/' jpeg_dir='/a/pippin01/Volumes/u08/lsst/CPanalysis/2005-05-'+Date+$ '/ShackHartman/jpegs/' Cxl=fltarr(25,3) Cxt=fltarr(25,3) Cyt=fltarr(25,3) Cyl=fltarr(25,3) Date_array=['10','11','12'] For i=0,2 do begin stats_directory='/a/pippin01/Volumes/u08/lsst/CPanalysis/2005-05-'+$ Date_array(i)+'/ShackHartman/stats_temp/' readcol,stats_directory+'average_xcorrelations.txt',Cxls,Cxts,$ format='(f,f)' readcol,stats_directory+'average_ycorrelations.txt',Cyls,Cyts,$ format='(f,f)' readcol,stats_directory+'average_offsets.txt',x_coord_avg,x_off_avg,$ y_coord_avg,y_off_avg,format='(x,f,f,f,f)' cxl(*,i)=cxls cyl(*,i)=cyls cxt(*,i)=cxts cyt(*,i)=cyts endfor set_plot,'Z' device, set_font='Courier' ;device, set_resolution=[8000,6000] ;!P.CHARSIZE=10 ;!P.CHARTHICK=15. ;!X.THICK=20. ;!Y.THICK=20. ;!P.THICK=50. ;!P.color=0 device,set_resolution=[800,600] !p.charsize=1 !p.charthick=1.5 !x.thick=2 !y.thick=2 !p.thick=5 !p.color=0 white='FFFFFF'x black='000000'x red='FF0000'x loadct,39 items=['5/10/05','5/11/05','5/12/05'] syms=[2,2,2] color_arr=[0,100,200] ;*************FIRST PLOT !P.multi=[0,1,1] plot,Cxt(*,0),psym=1,title="Average Transverse correlation along X direction",$ xtitle="n",ytitle="Covariance (pixels_squared)",$ background=white,color=black,$ yrange=[min(Cxt)-.5,max(Cxt)+.5],$ POSITION=[0.10,0.10,0.8,0.9] oplot,Cxt(*,1),color=100,psym=1 oplot,Cxt(*,2),color=200,psym=1 legend,items,psym=syms,colors=color_arr,/normal,position=[0.82,0.9],$ charthick=2,charsize=1 jpgimg=tvrd() tvlct,reds,greens,blues,/get write_png,jpeg_dir+'Average_XTtransverse_Correlations_05_'+Date+'.png',$ jpgimg,reds,greens,blues erase !p.multi=[0,1,1] ;***************SECOND PLOT plot,Cxl(*,0),psym=1,title="Average Longitudinal correlation along X direction",$ xtitle="n",ytitle="Covariance (Pixels Squared)",$ background=white,color=black,$ yrange=[min(Cxl)-.5,max(Cxl)+.5],$ POSITION=[0.10,0.10,0.8,0.9] oplot,Cxl(*,1),color=100,psym=2 oplot,Cxl(*,2),color=200,psym=2 legend,items,psym=syms,colors=color_arr,/normal,position=[0.82,0.9],$ charthick=2,charsize=1 jpgimg=tvrd() tvlct,reds,greens,blues,/get write_png,jpeg_dir+'Average_XLongitudinal_correlations_05_'+Date+'.png',$ jpgimg,reds,greens,blues erase !P.multi=[0,1,1] ;**************THIRD PLOT plot,Cyt(*,0),psym=2,title="Average Transverse correlation along Y direction",$ xtitle="n",ytitle="Covariance (Pixels Squared)",$ background=white,color=black,$ yrange=[min(Cyt)-.5,max(Cyt)+.5],$ POSITION=[0.10,0.10,0.8,0.9] oplot,Cyt(*,1),color=100,psym=2 oplot,Cyt(*,2),color=200,psym=2 legend,items,psym=syms,colors=color_arr,/normal,position=[0.82,0.9],$ charthick=2,charsize=1 jpgimg=tvrd() tvlct,reds,greens,blues,/get write_png,jpeg_dir+'Average_YTransverse_correlations_05_'+Date+'.png',$ jpgimg,reds,greens,blues erase ;***************FOURTH PLOT !p.multi=[0,1,1] plot,Cyl(*,0),psym=2,title="Average Longitudinal correlation along Y direction",$ xtitle="n",ytitle="Covariance (Pixels Squared)",$ background=white,color=black,$ yrange=[min(Cyl)-.5,max(Cyl)+.5],$ POSITION=[0.10,0.10,0.8,0.9] oplot,Cyl(*,1),color=100,psym=2 oplot,Cyl(*,2),color=200,psym=2 legend,items,psym=syms,colors=color_arr,/normal,position=[0.82,0.9],$ charthick=2,charsize=1 jpgimg=tvrd() tvlct,reds,greens,blues,/get write_png,jpeg_dir+'Average_YLongitudinal_correlations_05_'+Date+'.png',$ jpgimg,reds,greens,blues erase jpgimg=tvrd() write_jpeg,jpeg_dir+'Average_YLongitudnal_correlations_05_'+Date+'.jpg',congrid(jpgimg, 1600, 1200, /interp, /center), quality=100 erase ;!p.charthick=3 !p.thick=10 !p.multi=[0,1,1] fitted=where(x_off_avg ne 0) partvelvec,x_off_avg(fitted),y_off_avg(fitted),$ x_coord_avg(fitted),y_coord_avg(fitted),$ yrange=[100,820],xrange=[180,880],$ background=black,color=white,$ POSITION=[0.10,0.10,0.9,0.9],$ xtitle='x coordinate',ytitle='y coordinate',$ title='Average Offsets from Reference Grid' jpgimg=tvrd() write_jpeg,jpeg_dir+'Average_Offsets_05_'+Date+'.jpg',congrid(jpgimg, 1600, 1200, /interp, /center), quality=100 erase device,/close set_plot,'X' end