use wakes to improve `optimize_northing.py`

[aclerc]

At present optimize_northing.py makes the turbine yaw direction data agree with the chosen reanalysis data (normally ERA5 100m wind direction) however this is probably no better than +/- 5deg accurate since the real wind direction at wind farm hub height will not perfectly match ERA5.
The calculate north offset could be improved by finding wake centers in the data and comparing these to the expected wind direction of the wake centers.