PRO plt1arr
; Read in text file with RAFOS arrivals, sort the data and make some
; plots.
  !p.background='ffffff'x; white
  !p.color=0; black
  str1=''; input buffer
; Propagation times:	old V1:	11m 56s
;			new V1:	12m 30s
;			SS1:	4m 58s
;			SS2:	8m 43s
;			SS3:	4m 13s
;			Hoke:	7m 31s
  srcname=['Old V1', 'New V1', 'SS1', 'SS2', 'SS3', 'Hoke']
  ihrsrc=[4,4,4,5,5,6]
  tcensrc=[716.,750.,2098.,523.,2053.,451.]
  dtporm=20.; halfwidth of selection band to identify source
  snmin=5.; minimum signal to noise
  dtplot=2.; plot time width about tcensrc (plus and minus)
  nptmax=1000; allow for this many values
  year=intarr(nptmax)
  iyear=0
  day=intarr(nptmax)
  iday=0
  hour=intarr(nptmax)
  ihour=0
  min=intarr(nptmax)
  imin=0
  sec=fltarr(nptmax)
  fsec=0.
  width=fltarr(nptmax)
  fwidth=0.
  signal=fltarr(nptmax)
  fsignal=0.
  noise=fltarr(nptmax)
  fnoise=0.
  dayarr=dblarr(nptmax)
  ddayarr=double(0)

  infile='/usr/data/bland/pioneer/2001/r0124531.80a'; arrivals
  read,'infile = '+infile+':  ',str1; Change file name?
  IF strlen(str1) gt 0 THEN infile=str1; Use new file name.
  isrc=0; Default source.
  print,'Sources 0-'+strtrim(string(strlen(srcname)),2)+ $
    ' - choose one: ',srcname
  read,isrc
  read,'plot width (plus and minus) '+string(dtplot)+':  ', str1
  IF strlen(str1) gt 0 THEN dtplot=float(str1)
  read,'minimum signal to noise '+string(snmin)+':  ', str1
  IF strlen(str1) gt 0 THEN snmin=float(str1)

  openr,lun1,infile,/get_lun
  readf,lun1,str1; Titles
  print,'  ipt',str1
  ipt=0
  WHILE not(eof(lun1)) DO BEGIN
    readf,lun1,iyear,iday,ihour,imin,fsec, $
      fwidth,fsignal,fnoise,ddayarr
;    print,ipt,iyear,iday,ihour,imin,fsec, $
;      fwidth,fsignal,fnoise,ddayarr, $
;      format='(i5,i5,i4,2i3,2f8.4,2f8.2,f14.8)'

; **** Select some nice events.
    allsec=imin*60.+fsec; number of seconds after the hour
;    print,'ihour,allsec, s1SS1, s2SS1',ihour,allsec,s1SS1,s2SS1
;    IF fnoise eq 0 THEN BEGIN
;      print,'ERROR!!  zero noise'
;      print,ipt,iyear,iday,ihour,imin,fsec, $
;        fwidth,fsignal,fnoise,ddayarr, $
;        format='(i5,i5,i4,2i3,2f8.4,2f8.2,f14.8)'
;      read,'Press any key to continue',str1
;    ENDIF
    IF ihour mod 12 eq ihrsrc(isrc) and $
      allsec gt tcensrc(isrc)-dtporm and $
      allsec lt tcensrc(isrc)+dtporm and $
      fsignal/(fnoise+.01) gt snmin $
      THEN BEGIN
      sec(ipt)=allsec; 
      dayarr(ipt)=ddayarr; arrival time in days
      print,ipt,iyear,iday,ihour,imin,fsec, $
        fwidth,fsignal,fnoise,ddayarr, $
        format='(i5,i5,i4,2i3,2f8.4,2f8.2,f14.8)'
      ipt=ipt+1
    ENDIF
  ENDWHILE
  npt=max([ipt,1]); number of accepted points
  sec=sec(0:npt-1); truncate
  dayarr=dayarr(0:npt-1)
  plot,dayarr,sec,/ynozero,psym=2, $
    title=strmid(infile,strlen(infile)-12,12)+',  source '+srcname(isrc), $
    xtitle='day of 2001', $
    ytitle='travel time (sec)', $
    yrange=[tcensrc(isrc)-dtplot,tcensrc(isrc)+dtplot], $
    charsize=1.5

  save
END