code: def get_element_spd_dos(self, el):
          
          el = get_el_sp(el)
          el_dos = {}
          for site, atom_dos in self.pdos.items():
              if site.specie == el:
                  for orb, pdos in atom_dos.items():
                      orbital_type = _get_orb_type_lobster(orb)
                      if orbital_type not in el_dos:
                          el_dos[orbital_type] = pdos
                      else:
                          el_dos[orbital_type] = \
                              add_densities(el_dos[orbital_type], pdos)
  
          return {orb: Dos(self.efermi, self.energies, densities)
                  for orb, densities in el_dos.items()}
docs: Get element and spd projected Dos
  
  
          Args:
              el: Element in Structure.composition associated with LobsterCompleteDos
  
          Returns:
              dict of {Element: {"S": densities, "P": densities, "D": densities}}