# Reading total energies energies
e1 = nlread('Pd100_CO.hdf5', TotalEnergy)[0]
e2 = nlread('Pd100.hdf5', TotalEnergy)[0]
e3 = nlread('CO.hdf5', TotalEnergy)[0]

# computing adsorption energy
E = e1.evaluate().inUnitsOf(eV) - e2.evaluate().inUnitsOf(eV) - e3.evaluate().inUnitsOf(eV)
print("Adsorption energy = %.2f eV" % E)