PRO rdfit
  ;*** This routine reads in data from a scan for "B" calls, and data
  ;*** from Michael's fitting of the "B" calls.
  str1=''; input buffer

  infile='../bfit01_fc15p8.fit'; with kernel moved to 15.8 Hz
  fname='bfit01_fc15p8.fit'

  infile='/usr/data/bland03/pioneer/2091/244/bfit01.fit.21113'; AGU poster
  fname='bfit01.fit.30109'

  infile='/usr/data/bland03/pioneer/2091/244/bfit_rcen.fit.30109'; revised fit
  fname='bfit_rcen.fit.30109'

  infile='/usr/data/bland03/pioneer/2091/244/bfit_rcen5gate.fit.30110'; 2nd ref
  fname='bfit_rcen5gate.fit.30110'

  ;*** Result of mods with thrice iterated centering and the "cheating" fit.
  infile='/usr/data/bland03/pioneer/2091/244/bfit_5gcheat.fit.30122'; 3nd ref
  fname='bfit_5gcheat.fit.30122'

  ;*** Result of mods with thrice iterated centering and the "cheating" fit, 
  ;*** cheating fit used only if nchi2 better by .02, status flags revised.
  infile='/usr/data/bland03/pioneer/2091/244/bfit_5gch_p02.fit.30123'; 3nd ref
  fname='bfit_5gch_po2.fit.30123'
  ;*** Roger's file.
  rinfile='/usr/data/bland03/pioneer/2091/244/r9124413.bwc.21113'; AGU2002
  rfname='r9124413.b2c.21113'

  ;*** Michael's latest program, run on results from 16Hz-only kernel.
  infile='/usr/data/bland03/pioneer/2091/244/bfit_30122bwc.fit'; 3nd ref
  fname='bfit_30122bwc.fit'
  ;*** Roger's file.
  rinfile='/usr/data/bland03/pioneer/2091/244/r9124413.bwc.30122'; AGU2002
  rfname='r9124413.b2c.30122'

  ;*** Baby "B" calls scanned by hand, and fit to the first harmonic.
  ;*** Michael's file.
  infile='/Users/guest/seamount/project/0502/baby/bbscan/bbscan_feb2305.fit'
  fname='bbscan_feb2305.fit'
  ;*** Roger's file.
  rinfile='/usr/data/bland03/pioneer/2091/284/scanbb.bwc'
  rfname='scanbb.bwc'


  print,'Opening file '+infile
  openr,lunin,infile,/get_lun
  npt=long(0)
  readu,lunin,npt
  print,'npt = ',npt

  exits=lonarr(npt)
  index=lonarr(npt)
  f0=dblarr(npt)
  alpha=dblarr(npt)
  phi0=dblarr(npt)
  res2=dblarr(npt)
  nchi2=dblarr(npt)
  readu,lunin,exits,index,f0,alpha,phi0,res2,nchi2
  close,lunin
  free_lun,lunin

  IF 1 eq 0 THEN BEGIN
    print,'exits:  ',exits[0:10]
    print,'index: ',index[0:10]
    print,'f0:     ',f0[0:10]
    print,'alpha:  ',alpha[0:10]
    print,'phi0:   ',phi0[0:10]
    print,'res2:   ',res2[0:10]
    print,'nchi2: ',nchi2[0:10]
  ENDIF

  ;*** Read in Roger's file of results from "B" call scan.

  print,'Opening file '+rinfile
  openr,lunin,rinfile,/get_lun
  npeaks=long(0)
  readu,lunin,npeaks
  print,'npeaks = ',npeaks

  cfile=replicate('                                                  ',npeaks)
  iyear=intarr(npeaks)
  cday=dblarr(npeaks)
  cdb=fltarr(npeaks)
  ccen=fltarr(npeaks)
  csig=fltarr(npeaks)
  cthnsig=fltarr(npeaks)
  cthfx=fltarr(npeaks)

  readu,lunin,cfile,iyear,cday,cdb,ccen,csig,cthnsig,cthfx
  close,lunin
  free_lun,lunin
  save
END
  save
END