[プログラミング] TAIとUTCの差の近似計算らしきもの

id:yellow_73:20071030 の続き。
http://hpiers.obspm.fr/eop-pc/models/models.html
いるんですね、いろいろやってる人が。
この中で、TAI-UTC を計算する FORTRAN のコードがあります。

1960年1月1日から2003年12月31日までの、公表されている値を使うというものでした。
1973年1月1日のうるう秒制導入後は、ステップ状になっているので、年月日からTAI-UTCが出るテーブルを用意していました。
周波数・ステップ調整のころについては、刻みが不定でサンプリング的にデータをとって、線形補完しているようです。