;xy_offsets ; ;PURPOSE: ; This script will return the difference in coordinates between a ; ShackHartman image centroid and the corresponding point on the ; reference grid ; ; Since the center of the grid seems to be drifting pro cn_full_log_jpegs,Date=Date,stats_directory=stats_directory,$ x_center_file=x_center_file,y_center_file=y_center_file,$ depth=depth,x_off=x_off,y_off=y_off If (not keyword_set(Date)) then Date='10' If (not keyword_set(stats_directory)) then stats_directory='stats' If (not keyword_set(x_center_file)) then x_center_file='x_centers.txt' If (not keyword_set(y_center_file)) then y_center_file='y_centers.txt' If (not keyword_set(depth)) then depth=1 If (not keyword_set(postscript_dir)) then postscript_dir='postscripts' If (not keyword_set(x_off)) then x_off=0 If (not keyword_set(y_off)) then y_off=0 stats_directory='/a/pippin01/Volumes/u08/lsst/CPanalysis/2005-05-'+Date+$ '/ShackHartman/'+stats_directory+'/' postscript_dir='/a/pippin01/Volumes/u08/lsst/CPanalysis/2005-05-'+Date+$ '/ShackHartman/'+postscript_dir+'/' Cor_directory='/a/pippin01/Volumes/u08/lsst/CPanalysis/2005-05-'+Date+$ '/ShackHartman/correlations/' jpeg_dir='/a/pippin01/Volumes/u08/lsst/CPanalysis/2005-05-'+Date+$ '/ShackHartman/jpegs/' num=25 Cxt_means=fltarr(num,depth) Cxl_means=fltarr(num,depth) Cyl_means=fltarr(num,depth) Cyt_means=fltarr(num,depth) for m=0, depth-1 do begin readcol,stats_directory+'Depth_'+string(strtrim(m,2))+'Correlations.ps',$ Cxl,Cxt,Cyl,Cyt,format='f,f,f,f' Cxt_means(*,m)=cxt Cxl_means(*,m)=cxl Cyl_means(*,m)=cyl Cyt_means(*,m)=cyt endfor set_plot,'Z' loadct,39 ;load color table 39 device, set_font='Courier' device,set_resolution=[800,600] !p.charsize=1 !p.charthick=1.5 !x.thick=2 !y.thick=2 !p.thick=3 !p.color=0 white='FFFFFF'x black='000000'x items_time=['t=0'] sym_time=[1] items_k_tran=['Data Points','n!e-7/32!n','n!e-13/64!n','n!e-27/128!n'] sym_k_tran=[1,0,0,0] line_k_tran=[0,1,2,3] items_k_lon=['Data Points','n!e-1/3!n'] sym_k_long=[0,1] colors_time=[0] ;*************FIRST PLOT !P.multi=[0,2,2] plot,Cxt_means(*,0),psym=1,title="Transverse X",$ xtitle="n",ytitle="Covariance (Pixels Squared)",$ background=white,color=black,$ yrange=[1e-3,1e1],$ position=[0.5,0.55,0.8,0.9], /ylog ;Over Plot the various separations in time and also generate legends for i=1,depth-1 do begin oplot,Cxt_means(*,i),psym=1,color=20*i sym_time=[sym_time,1] items_time=[items_time,string(strtrim(i*7,2))+' seconds'] colors_time=[colors_time,20*i] endfor ;***************FIFTH PLOT plot,Cxl_means(*,0),psym=1,title="Longitudinal X",$ xtitle="n",ytitle="Covariance (Pixels Squared)",$ background=white,color=black,$ yrange=[1e-3,1e1],$ position=[0.1,0.55,0.4,0.9],/ylog for i=1,depth-1 do begin oplot,Cxl_means(*,i),psym=1,color=20*i endfor ;**************NINTH PLOT plot,Cyt_means(*,0),psym=1,title="Transverse Y",$ xtitle="n",ytitle="Covariance (Pixels Squared)",$ background=white,color=black,$ yrange=[1e-3,1e1],$ position=[0.5,0.1,0.8,0.45],/ylog for i=1,depth-1 do begin oplot,Cyt_means(*,i),psym=1,color=20*i endfor ;*************THIRTEENTH plot,Cyl_means(*,0),psym=1,title="Longitudinal Y",$ xtitle="n",ytitle="Covariance (Pixels Squared)",$ background=white,color=black,$ yrange=[1e-3,1e1],$ position=[0.1,0.1,0.40,0.45],/ylog for i=1,depth-1 do begin oplot,Cyl_means(*,i),psym=1,color=20*i endfor legend,items_time,psym=sym_time,position=[.82,.9],/normal,colors=colors_time jpgimg=tvrd() tvlct,reds,greens,blues,/get write_png,jpeg_dir+'Log_time_correlations'+Date+'.png',$ jpgimg,reds,greens,blues erase device,/close set_plot,'X' end