Fix typo in solve documentation
Use surface_velocity if velocity is NaN