% file: bp_read_fits.i
%
% Read in the bpl structure values from a FITS bin table.
% A little tedious since I have bpl as an array
% of structures and not a structure of arrays...
%

%
% There are no string columns in the saved FITS file
% so makeup defaults...
%

define bp_read_fits (filename)
{
  % setup the column arrays...
  variable columns = struct { 
  id,
  ignore,
  lam,
  grp,
%  desc,
  lfit,
  lclo,
  lchi,
  afit,
  aclo,
  achi,
  wfit,
  wclo,
  wchi,
  cfit,
  cclo,
  cchi,
  chifit  };
%  rgfit   };

  % Read in the table
  columns = fits_read_table(filename);

  % Loop over the elements filling bpl structure...
  variable iloop=0;
  while (iloop < length(columns.id)) {

  bpl[iloop].id = columns.id[iloop];
  bpl[iloop].ignore = columns.ignore[iloop];
  bpl[iloop].lam = columns.lam[iloop];
  bpl[iloop].grp = columns.grp[iloop];
%  bpl[iloop].desc = bpl[iloop].desc;
  bpl[iloop].lfit = columns.lfit[iloop];
  bpl[iloop].lclo = columns.lclo[iloop];
  bpl[iloop].lchi = columns.lchi[iloop];
  bpl[iloop].afit = columns.afit[iloop];
  bpl[iloop].aclo = columns.aclo[iloop];
  bpl[iloop].achi = columns.achi[iloop];
  bpl[iloop].wfit = columns.wfit[iloop];
  bpl[iloop].wclo = columns.wclo[iloop];
  bpl[iloop].wchi = columns.wchi[iloop];
  bpl[iloop].cfit = columns.cfit[iloop];
  bpl[iloop].cclo = columns.cclo[iloop];
  bpl[iloop].cchi = columns.cchi[iloop];
  bpl[iloop].chifit = columns.chifit[iloop];
  bpl[iloop].rgfit = filename;

  iloop = iloop +1;
  }

}