Extract solution values on surface nodes

See for instance: Fenicsx: use an array of data for boundary conditions - #5 by dokken