pro pc_read_saffman,tt,iir
;
;  Simple routine for reading Saffman integrals.
;  Here with hardwired file name and file size.
;
file='data/Iv_bcs_saffman_mag.dat'
n=1024L/2
;
ir=fltarr(n)
openr,1,file
;
;  Initialize counter, read, and append, until end of file (eof):
;
it=0L
while not eof(1) do begin
  readf,1,t
  readf,1,ir
  print,t
  if it eq 0 then begin
    tt=t
    iir=ir
  endif else begin
    tt=[tt,t]
    iir=[iir,ir]
  endelse
  it=it+1
endwhile
nt=it
;
;  Reform
;
iir=reform(iir,n,nt)
close,1
END