code: def fhi_header(filename, ppdesc):
          
          lines = _read_nlines(filename, 4)
  
          try:
              header = _dict_from_lines(lines[:4], [0, 3, 6, 3])
          except ValueError:
              # The last record with rchrg ... seems to be optional.
              header = _dict_from_lines(lines[:3], [0, 3, 6])
  
          summary = lines[0]
  
          return NcAbinitHeader(summary, **header)
docs: Parse the FHI abinit header. Example:
  
          Troullier-Martins psp for element  Sc        Thu Oct 27 17:33:22 EDT 1994
           21.00000   3.00000    940714                zatom, zion, pspdat
             1    1    2    0      2001    .00000      pspcod,pspxc,lmax,lloc,mmax,r2well
          1.80626423934776     .22824404341771    1.17378968127746   rchrg,fchrg,qchrg