Implement Jayne and St Laurent Internal wave drag parameterization (with tides)

Context

With explicit tides on, unless high spatial resolution is used, the internal wave spectrum is not fully resolved. As a matter of fact, only a fraction of the conversion from barotropic to baroclinic tides is simulated which substantially degrades the accuracy of external (and eventually internal) tides. The Jayne and St Laurent (2001) parameterization represents this missing energy sink for tides in the form of an additional bottom drag.

https://agupubs.onlinelibrary.wiley.com/doi/10.1029/2000GL012044

Proposal

• Read a new internal wave linear drag coefficient (units of m/s) in the following form: Ciwd = Hr^2 / L * Nb, where Hr is the unresolved topography roughness, Nb, the bottom stratification and L a characteristic length scale
• Implement additionnal drag in the barotropic equations on high frequency (tidal) motions only. This requires the use of an online filtering of barotropic velocities (25h time averaging) to consider the principal semidiurnal tides only.

---