pro Tracking_hr5132a

;This script will attempt to figure out how large of an offset there is between
;the position of hr5132a and the pointing of the telescope

CP_dir='/a/pippin01/Volumes/u08/lsst/CPanalysis/'

;Celestial Coordinates of hr5132a:
;				
;RA=       13 39 53.26		204.972
;DEC=     -53 27 59.00	       -53.466389

RA=204.972
DEC=-53.466389

;Latitude and Longitude of Cerro Pachon
;
;Latitude=	30 14 16.8	30.238000
;Longitude=	70 44 01.4	70.733722
;Altitude=	2738 m

Latitude=-30.238000
Longitude=-70.733722
height=2738

readcol,CP_dir+'SOAR_tracking.txt',dates,times,azimuth,altitude,$
	format='a,a,f,f'


Times_arr=fltarr(6,N_elements(dates))
jd_arr=dblarr(N_elements(dates))
ra_arr=fltarr(N_elements(dates))
dec_arr=fltarr(N_elements(dates))
ra_arr(*)=ra
dec_arr(*)=dec

for j=0, N_elements(times)-1 do begin
                                                                                
     time_arr=strsplit(times(j),':',/extract)
     date_arr=strsplit(dates(j),'-',/extract)
     Times_arr(0,j)=date_arr(0)
     Times_arr(1,j)=date_arr(1)
     Times_arr(2,j)=date_arr(2)
     Times_arr(3,j)=time_arr(0)
     Times_arr(4,j)=time_arr(1)
     Times_arr(5,j)=time_arr(2)
     yr_val=valid_num(date_arr(0),year)
     month_val=valid_num(date_arr(1),month)
     date_val=valid_num(date_arr(2),obsdate)
     sec_val=valid_num(time_arr(2),secs)
     min_val=valid_num(time_arr(1),mins)
     hr_val=valid_num(time_arr(0),hrs)

     if (hrs-4) lt 0 then begin
	obsdate=obsdate-1
        hrs=hrs+24
     endif else begin
	hrs=hrs
     endelse
     juldate,[year,month,obsdate,hrs,mins-3,secs],jd
     jd_arr(j)=jd+2400000
                                                     
   endfor

eq2hor,ra_arr,dec_arr,jd_arr,star_alt,star_az,$
      lat=latitude,lon=longitude,abberation_=0,refract_=0
;hadec2altaz,ra_arr,dec_arr,latitude,alt,az

set_plot,'ps'
device,filename='/a/pippin01/Volumes/u08/lsst/CPanalysis/Tracking.ps'
!p.multi=[0,1,2]

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

date=[11,12,13]

for i=0,2 do begin
  plot,star_alt(where(times_arr(2,*) eq date(i))),yrange=[40,80],$
	xtitle='time',ytitle='Altitude'
  oplot,altitude(where(times_arr(2,*) eq date(i))),psym=2
  plot,star_az(where(times_arr(2,*) eq date(i))),yrange=[140,230],$
	xtitle='time',ytitle='Azimuth'
  oplot,azimuth(where(times_arr(2,*) eq date(i))),psym=2
endfor

device,/close

set_plot,'z'
erase
device, set_font='Courier'
device,set_resolution=[700,600]
!p.charsize=.8
!p.charthick=1.2
!x.thick=2
!y.thick=2
!p.thick=1
!p.noerase=0

;For Presentation Purposes
device,set_resolution=[1400,1200]
!p.charsize=3
!p.charthick=5
!x.thick=4
!y.thick=4
!p.thick=3


red='FF0000'x
black='000000'x
white='FFFFFF'x

for i=0,2 do begin
 !p.charsize=3
 !p.charthick=5
 !x.thick=4
 !y.thick=4
 !p.thick=3
 
 plot,star_alt(where(times_arr(2,*) eq date(i))),yrange=[40,80],$
        xtitle='Time (arbitrary units)',ytitle='Altitude',color=black,$
	background=white,$
	title='Tracking Errors for 5/'+strtrim(Date(i)-1,2)+'/05'
 oplot,altitude(where(times_arr(2,*) eq date(i))),psym=2,color=black
  plot,star_az(where(times_arr(2,*) eq date(i))),yrange=[140,230],$
        xtitle='Time (arbitrary units)',ytitle='Azimuth',color=black,$
	background=white
  oplot,azimuth(where(times_arr(2,*) eq date(i))),psym=2,color=black


items=['Star Coordinate','Tracking Coordinate']
sym_arr=[0,2]

!p.charsize=2.5
!p.charthick=3
!x.thick=2
!y.thick=2
!p.thick=1
!p.noerase=0
legend,items,psym=sym_arr,colors=[0,0],textcolors=[0,0],pos=[.55,.9],$
	/norm,outline_color=0
jpgimg=tvrd()
tvlct,reds,greens,blues,/get
write_png,'/a/pippin01/Volumes/u08/lsst/CPanalysis/'+strtrim(date(i)-1,2)+$
	'Tracking.png',$
jpgimg,reds,greens,blues

endfor


erase

set_plot,'x'

stop

end