Take the boundary free energy into account

[urbainvaes]

In addition to the bulk free energy, the code should calculate the free energy related to the boundary.

Suggested workflow:

• In the problem specification step, one should create two arrays: one with the labels of all the surfaces / lines of the problem, and another one with the corresponding spreading parameters.
• In the main solver file, one can them loop over them.