;+
;   BCS_LC_PLOT
;
;   To plot a better BCS light curve.  Uses output of GS_LC and can fit splines
;   as an overlay.
;
;  CDP August 92
;
;-

n1 = dset_str(gs_ds1).dset
n2 = dset_str(gs_ds2).dset

if !d.name ne 'PS' then begin
   anst = ' '
   read,'Give plot title> ',anst
endif

!y.title = 'counts/sec'
ymax = max(line)*1.2
ymin = -0.1*ymax

repeat begin
   utplot,roadmap(n1:n2),line,title=anst,yrange=[ymin,ymax]
   oplot,[0,100000],[0,0]
   if !d.name ne 'PS' then begin
      read,'Number of splines to fit? (0 to skip) > ',sn
   endif
   if sn gt 0 then begin
      x = findgen(n_elements(line))
      y = bells(x,line,x,sn)
      outplot,roadmap(n1:n2),y,thick=3
      if !d.name ne 'PS' then begin
         ans = ' '
         read, 'Fit OK? (Y/N/Quit) > ',ans
         if strupcase(ans) eq 'Q' then ans = 'y'
         yesno,ans
      endif else begin
         ans = 1
      endelse
   endif  else begin
      ans = 1
   endelse
endrep until ans eq 1


end