Cosmetic changes to better match the format of the other functions
Add new output which stores the number of elements in each model boundary
Change the interpolation to occur at the boundary nodes to the centroid of the boundary elements (since FVCOM work on velocities at the centre of each element)
Transfer output from open boundary nodes to open boundary elements