pro plot_track,flist,xrange=xrange,yrange=yrange,labax=labax,$
title=title,paper=paper,oplot=oplot

;Little script to plot tracks from proja2034.
;INPUTS: list of file names. The RA and Dec for each integration is plotted
;Xrange /Yrange  - the xrange and yrange can be specified, the default
;                  is the range of the first file
;labax: select to apply "RA" and "Dec" as x and y axis labels.

;This program is very slow, since the files are big and takes a long time to open them, 
;find the info, and then close them, this is useful to use in conjunction with 
;driftmap_planner, when trying to plan out mulit-day drift observations.


nfiles=n_elements(flist)

if keyword_set(paper) then set_plot,'ps'

for thisfile=0,nfiles-1 do begin
fname=flist(thisfile)
    ftab_ext,fname,"CRVAL2",ra
    ftab_ext,fname,"CRVAL3",dec

ra=ra/15.


if thisfile eq 0 and not keyword_set(oplot) then begin
    if not keyword_set(xrange) then xrange=reverse(minmax(ra))
    if not keyword_set(yrange) then begin
        yrange=minmax(dec)
        yrange(0)=yrange(0)-nfiles*13/(60.*2)
    endif
    if not keyword_set(title) then title=""
    if keyword_set(labax) then begin
        xtitle="RA hrs"
        ytitle="Dec deg"
    endif else begin
        xtitle=""
        ytitle=""
    endelse
    plot,ra,dec,xrange=xrange,yrange=yrange,psym=3,$
      xtitle=xtitle,ytitle=ytitle,title=title,xstyle=1,ystyle=1
endif else begin
    
    oplot,ra,dec,psym=3
endelse

endfor


if keyword_set(paper) then begin
    device,/close
    set_plot,'x'
endif

end