Split up the tides function into a series of generic functions (which can hopefully be reused for different sources of harmonics). Also interpolate the harmonics to the positions rather than searching for the nearest value. The interpolation should yield smoother harmonics over small distances.
Split up the tides function into a series of generic functions (which can hopefully be reused for different sources of harmonics). Also interpolate the harmonics to the positions rather than searching for the nearest value. The interpolation should yield smoother harmonics over small distances.

pwcazenave/PyFVCOM