;+
; file: ft_chipx_plot.pro


; Read focus trends file directly:
;;ft = rdb_read(!DDIDL+"/marx/focus_trends.rdb")

; Or get it from focus_trends.pro before running this:
;
;  IDL> focus_trends, /NOM_AIM, ft=ft
;  IDL> @ft_chipx_plot
;
;-


; Make a year value
juldayval = LONG(indgen(n_elements(ft))+1.0)
for id=0, n_elements(ft)-1 do juldayval(id) = JULDAY(ft(id).month, ft(id).day, ft(id).year)
julrefdate = JULDAY(1,1,1999)
juldayval = DOUBLE(juldayval - julrefdate)
yrval = 1999.0 + juldayval/365.25

; Plot the CHIPX values in the focus trends file...

plot, yrval, ft.chipx, psym=4, yrange=[100.,350.],/ysty, $
	XRANGE=[1999.5, 2005.0], XSTYLE=1, $
	XTITLE='Observation Date ( ' + $
	STRCOMPRESS(FIX(ft(0).month),/REMOVE) + '/' + $
		STRCOMPRESS(FIX(ft(0).year),/REMOVE) + ' to ' + $
	STRCOMPRESS(FIX(ft(n_elements(ft)-1).month),/REMOVE) + '/' + $
		STRCOMPRESS(FIX(ft(n_elements(ft)-1).year),/REMOVE) + $
	' )', $
	YTITLE='CHIPX location of object on S3', $
	TITLE='CHIPX of Object versus Observation Date'

; Show where 256 (quad boundary) is as well as
; +/-20" from some point...
oplot, [0.,5000.0], 256.0*[1.,1.], linestyle=2
oplot, [0.,5000.0], 20.0*2.0+246.0*[1.,1.], linestyle=1
oplot, [0.,5000.0], -20.0*2.0+246.0*[1.,1.], linestyle=1

; 
xyouts, 2002.0, 260.0, 'Quadrant boundary (CHIPX=256)'

xyouts, 2000.2, 160.0, 'Dotted lines are +/- 20 arc sec from CHIPX=246'

; show July 1 2003
oplot, 2003.0+(7.-1)/12.0+[0.,0.], [0.,1024.], linestyle=1

; and the boundary where consistent pointing started
oplot, 2001.06+[0.,0.], [0.,1024.], linestyle=1

xyouts, 2003.5-0.15, 130.0, '7/1/03'