PRO sp2rd, file, x, y, header

; Read a 1d function from a spectrum format 1 file

  openr,lun,file,/get_lun
  header  = 'a'
  header2 = 'b'
  ndat = 1L
  ReadF, lun, format='(A)', header
  ReadF, lun, format='(A)', header2
  ReadF, lun, ndat

  x=DblArr(ndat)
  y=DblArr(ndat)
  z=DblArr(2L*ndat)

  i=0L & valid = 0
  ReadF, lun, z

  FOR i = 0L,ndat-1 DO BEGIN
    x(i)= z(2L*i) &  y(i)= z(2L*i+1)
  ENDFOR

  close,lun
  free_lun, lun

END