;+
; NAME:
;    TERM_QUICK
; PURPOSE:
;    Show the SFC pointing locations graphically for a specific
;      time range. This program feeds information to QUICK.HTML.
;      To monitor the OP functions of the terminator acquisition,
;      use OP_TERM_SCORE
; CATEGORY:
; CALLING SEQUENCE:
; INPUTS:
; OPTIONAL (KEYWORD) INPUT PARAMETERS:
;    radius, makes a plot of apparent solar radius at each 
;      terminator time
;    list, makes a list of SFC filter positions
;    tstart, start time (default currently 11-Feb-00)
;    tstop, stop time (default !stime)
; ROUTINES CALLED
;    FINDFILE, GT_FILTA, GT_FILTB, RD_ROADMAP, !stime, TIM2DSET
; OUTPUTS:
; COMMON BLOCKS:
; SIDE EFFECTS:
;    LINECOLORS
; RESTRICTIONS:
; MODIFICATION HISTORY:
;    HSH, written 5-Jul-99 based on TERM_REVIEW_GIF
;    11-Jul-99, /list added (HSH)
;    16-Jul-99, plottime + charsize=1.25 (BNH)
;    19-Jul-99, fixed bug (HSH) and added /stats
;    14-Aug-99, adjusted range and title of plot (HSH)
;    14-Feb-00, revised according to new normal pointing (HSH)
;    19-Aug-99, adjusted plot (HSH)
;    12-jan-00, y2k fix
;    12-feb-00, update plot range (HSH)
;     6-mar-00, update plot range (HSH)
;    14-apr-00, update (HSH)
;-

pro term_quick, tstart=tstart, tstop=tstop, qstop = qstop, $
  radius = radius, list = list, stats = stats, set = set
sfcs_t=file_list('$DIR_SXT_SFC','sfc*',/cd)
sfcs=sfcs_t(sort(file2time(sfcs_t))) ;y2k
rd_roadmap,sfcs,rm_sfc
if n_elements(tstart) eq 0 then tstart = '2-jul-99'
if n_elements(tstop)  eq 0 then tstop  = !stime
t1_sfc  = tim2dset(rm_sfc, tstart)
t2_sfc  = tim2dset(rm_sfc, tstop)
sfcs = sfcs(t1_sfc:t2_sfc)
n_sfc = n_elements(sfcs)
rd_sda,sfcs,-1,ind_sfc,/nodata
if keyword_set(list) then begin
  filt = intarr(6,2)
  for i = 0, n_sfc-1 do begin
    ii = gt_filta(ind_sfc(i))
    jj = gt_filtb(ind_sfc(i))
    if ii eq 1 then filt(jj-1,0) = filt(jj-1,0)+1
    if ii eq 6 then filt(jj-1,1) = filt(jj-1,1)+1
  endfor
  print,'             ',gt_filtb(1+indgen(6),/str)
  print,'Open ',filt(*,0)
  print,'ND   ',filt(*,1)
endif
xy = ind_sfc.his.sun_center
xy = xy(0:1,*)
xyref = [500,590]
if keyword_set(radius) eq 0 then begin
  plot,[xyref(0),xyref(0)],[xyref(1),xyref(1)],psym=2,symsiz=3,$
    xr=[490,510],yr=[580,600],xst=1,yst=1,$
    charsi=1.5
  oplot,[490,510],[xyref(1),xyref(1)],lines=1
  oplot,[xyref(0),xyref(0)],[580,600],lines=1
  xyold = [537,585]
  oplot,[xyold(0),xyold(0)],[xyold(1),xyold(1)],psym=1,symsiz=3
  oplot,[xy(0,*),xy(0,*)],[xy(1,*),xy(1,*)],psym=5,symsiz=2
  thinal = thinal_index(ind_sfc)
  oplot,[xy(0,thinal),xy(0,thinal)],[xy(1,thinal),xy(1,thinal)],$
    psym=5,symsiz=2,color=7
  lastsfd,i_sfd,d_sfd,/notv,30
  xy_sfd = sxt_cen(i_sfd)
  oplot,[xy_sfd(0,*),xy_sfd(0,*)],[xy_sfd(1,*),xy_sfd(1,*)],psym=4
  plottime	; Added 6-jan-99 by DMcK
endif else begin
  rr = get_rb0p(ind_sfc,/r)
  utplot,ind_sfc,rr,psym=2,/ynoz,symsiz=2,$
    tit='Apparent solar radius',charsi=1.25,xst=0,$
    ytit='Apparent radius, arc sec'
  ss = thinal_index(ind_sfc)
  outplot,ind_sfc(ss),rr(ss),psym=2,symsiz=2,color=7
  outplot,['11-feb-00','11-feb-00'],[940,980],lines=1
  plottime
endelse
if keyword_set(stats) then begin
  sig_sfd_x = stdev(xy_sfd(0,*),sfd_x)
  sig_sfd_y = stdev(xy_sfd(1,*),sfd_y)
  print,sfd_x,' +- ',sig_sfd_x,' ?=> ',xyref(0)
  print,sfd_y,' +- ',sig_sfd_y,' ?=> ',xyref(1)
endif
if keyword_set(set) then begin
  i = 9
;  print,' '
;  print,"Showing the last ThinAl SFC's (. and *)  and the last 10 SFD's (+) for reference"
;  print,' '
  while i lt 10 do begin
    i = i+1
    wdef,1
    thal = thinal_index(rm_sfc)
    plot,xy(0,*),xy(1,*),psym=5,symsiz=2,xr=[490,525],yr=[585,605]
    oplot,xy_sfd(0,*),xy_sfd(1,*),psym=4,symsiz=2
    nn = n_sfc-1
    oplot,[xy(0,nn),xy(0,nn)],[xy(1,nn),xy(1,nn)],$
      psym=2
    print,'Select best guess at current pointing'
    cursor,x1,y1,/data & print,x1,y1
    pause
    print,'Select target'
    cursor,x2,y2,/data & print,x2,y2
    pause
    pix = gt_pix_size()
    print,'Increment Euler angles (X, Y) by: '
    print,   (y2-y1)*pix,  (x2-x1)*pix, ' arc sec'
    arrow,x1,y1,x2,y2,/data
  endwhile
endif
clearplot
clear_utplot
if keyword_set(qstop) then stop
end