使用S_TIDE计算分潮的交点因子和订正角

使用S_TIDE工具包里的s_nodal_cal.m函数来计算分潮的交点因子和订正角，调用程序如下：

[ff,uu]=s_nodal_cal({'K1';'O1';'M2'},[1920,01,01],[2020,01,01],365,10);

subplot(2,1,1);

plot(ff(1,:),'k');hold on;plot(ff(2,:),'r');plot(ff(3,:),'b')

legend('K1','O1','M2');

set(gca,'XTick',0:25:2020)

set(gca,'XTickLabel',{'1920','1945','1970','1995','2020'})

title('(a) nodal factor')

subplot(2,1,2);

plot(uu(1,:),'k');hold on;plot(uu(2,:),'r');plot(uu(3,:),'b')

set(gca,'XTick',0:25:2020)

set(gca,'XTickLabel',{'1920','1945','1970','1995','2020'})

title('(b) nodal angle')

结果如下：

需要注意的是输入的分潮名必须用cell数组，365代表计算的时间间隔为1年（365天）。10代表纬度为北纬10度，这个参数不重要，只要不为0即可。