潮位调和分析工具学习（2）——使用T_tide时遇到的问题

遇到的问题

在处理长时间（大于18.61年）的潮位时间序列时，我们发现 T_tide 分析所得迟角的逐年变化规律十分奇怪。
以我国南海某潮位站点的长时间序列为例，利用 T_tide 对潮位序列做逐年的分析，得到主要分潮M₂的振幅与迟角的逐年变化（如下图所示）。

原因

主要分潮参数存在着周期18.61年的波动。所以对于时间大于18.61年的潮汐序列，逐年调和分析所得到的振幅和迟角需要做nodal校准。我们需要对调和分析的振幅和迟角需要做18.61年的拟合，以得到实际的nodal cycle。
与之对应，在使用T_tide做逐年潮位分析时，我们输入的时间总是0 ~ 365 days（或366 days）。实际上只有第一年的时间可以使用0 ~ 365 days，而随后第二年的时间得用 365 ~ 730 days（后续年份的输入时间以此类推）。

长期潮位序列调和分析中对迟角度18.61年循环的处理（参考潘海东博士的博客）

为解决上述问题，我们需对 t_tide.m 做一些修改；我们需要把数据的起始与终止时刻也作为输入参数。修改部分如下：

1. 函数头可修改为：

[nameu,fu,tidecon,xout,z0]=t_tide2(xin,aa,bb,varargin)
% aa 和 bb 分别表示序列的起始和终止时刻

2. 原 t_tide 程序的298行。将原来的代码 " t=dt*([1:nobs]'-ceil(nobsu/2)); " 修改为：

t=dt*([aa:bb]');

之后，使用 t_tide2 函数进行调和分析。得到结果如下：

最后附上潘海东博士的 b站主页 和 ResearchGate主页。
潘博士开发了一个功能强大且十分好用的新潮汐分析工具包S_tide。相关内容可参考潘博士的论文和博客！