  We construct a contour function for the entanglement entropies in generic
harmonic lattices. In one spatial dimension, numerical analysis are performed
by considering harmonic chains with either periodic or Dirichlet boundary
conditions. In the massless regime and for some configurations where the
subsystem is a single interval, the numerical results for the contour function
are compared to the inverse of the local weight function which multiplies the
energy-momentum tensor in the corresponding entanglement hamiltonian, found
through conformal field theory methods, and a good agreement is observed. A
numerical analysis of the contour function for the entanglement entropy is
performed also in a massless harmonic chain for a subsystem made by two
disjoint intervals.
