Triplett line energies

require("isisscripts");

%read in data from the drake 1998 table
variable dk = ascii_read_table("drake_table3_cm.dat",
                               [{"%d"},{"%F"},{"%F"},{"%F"},{"%F"},{"%F"}]); 

%conversion factor from cm^-1 to keV
variable cf  = 8065730;

%list of elements
variable elem = ["He","Li","Be","B","C","N","O","F","Ne","Na","Mg","Al",
                 "Si","P","S","Cl","Ar","K","Ca","Sc","Ti","V","Cr","Mn",
                 "Fe","Co","Ni","Cu","Zn","Ga","Ge","As","Se","Br","Kr"];

%define structure for the lines; line conversion from the values in
%Drake is according to what Natalie told me what the individual columns
%in the table represent
variable tri = struct{element = elem,
  w = dk.col3/cf,
  y = dk.col2/cf,
  x = (dk.col6-dk.col5+dk.col2)/cf,
  z = (dk.col2 - dk.col5)/cf};

%write fits fite and ascii file
fits_write_binary_table("tripletts.fits","DATA",tri);
print_struct("tripletts.dat",tri);

print("STACK:");
_print_stack;
exit;