Identifying tidal constituents with spectral analysis

R code for finding the dominant tidal constituent from a time series of water depth data. Minor constituents can be found by converting the respective frequency components of local amplitude maxima into a period, which I leave up to you.