% file: bp_write_fits.i
%
% Output the bpl structure to a FITS bin table.
% A little tedious since I have bpl as an array
% of structures and not a structure of arrays...
%

%
% Doesn't work for string columns (array_struct_field? the FITS write?)
%

define bp_write_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   };

  columns.id = array_struct_field(bpl,"id");
  columns.ignore = array_struct_field(bpl,"ignore");
  columns.lam = array_struct_field(bpl,"lam");
  columns.grp = array_struct_field(bpl,"grp");
%  columns.desc = array_struct_field(bpl,"desc");
  columns.lfit = array_struct_field(bpl,"lfit");
  columns.lclo = array_struct_field(bpl,"lclo");
  columns.lchi = array_struct_field(bpl,"lchi");
  columns.afit = array_struct_field(bpl,"afit");
  columns.aclo = array_struct_field(bpl,"aclo");
  columns.achi = array_struct_field(bpl,"achi");
  columns.wfit = array_struct_field(bpl,"wfit");
  columns.wclo = array_struct_field(bpl,"wclo");
  columns.wchi = array_struct_field(bpl,"wchi");
  columns.cfit = array_struct_field(bpl,"cfit");
  columns.cclo = array_struct_field(bpl,"cclo");
  columns.cchi = array_struct_field(bpl,"cchi");
  columns.chifit = array_struct_field(bpl,"chifit");
%  columns.rgfit = array_struct_field(bpl,"rgfit");

  variable keywords = struct {creator};
  keywords.creator = "bp_write_fits.i";
  
  % and write it out...
  fits_write_binary_table(filename, "BUMPS", columns, keywords);

}