code: def get_interpolated_value(self, energy, integrated=False):
          
          inter = {}
          for spin in self.cohp:
              if not integrated:
                  inter[spin] = get_linear_interpolated_value(self.energies,
                                                              self.cohp[spin],
                                                              energy)
              elif self.icohp is not None:
                  inter[spin] = get_linear_interpolated_value(self.energies,
                                                              self.icohp[spin],
                                                              energy)
              else:
                  raise ValueError("ICOHP is empty.")
          return inter
docs: Returns the COHP for a particular energy.
  
          Args:
              energy: Energy to return the COHP value for.