read_elvlc,'mg_2.elvlc_sigut',l1,term,conf,ss,ll,jj,ecm,eryd,ecmth,erydth,ref

read_wgfa_str,'mg_2.wgfa_ff',wgfastr,ref

openr,lin,'mg2_sigut_data_fs.txt',/get_lun
openw,lout,'mg_2.upsdat_sigut_fs',/get_lun

t=[0.1,0.3,0.5,1.0,2,3]*1e4

str1=''
WHILE eof(lin) NE 1 DO BEGIN
  readf,lin,str1
  bits=str_sep(str1,' ')
  l1=fix(bits[0])
  l2=fix(bits[1])
  de=abs(eryd[l2-1]-eryd[l1-1])
  ups=float(bits[2:7])
  printf,lout,format='(2i5,f10.6,6e10.3)',l1,l2,de,t
  k=where(wgfastr.lvl1 EQ l1 AND wgfastr.lvl2 EQ l2,nk)
  IF nk EQ 1 THEN gf=wgfastr[k].gf ELSE gf=0.
  
  printf,lout,format='(2i5,f10.6,6e10.3)',l1,l2,gf,ups
ENDWHILE

free_lun,lin,lout

END