The UCLA & IRD Brest users (e.g., Penven, Marchesiello and me) have a similar case. Up to now we use FFT on the output fields to get the tidal constituents and tidal mean. Since you then need a series that is long enough and of sufficiently high temporal resolution to separate the components, it...