Compute heads and velocities on unstructured grids

[dbrakenhoff]

Adapt head_grid and velocity_grid to use head_array() and velocity_array() that accept 1D x, y inputs.

• head_array is done in add head_array() methods for unstructured computation of heads #139
• velocity_array still to be added