有许多模型估计轨道参数耗时过长。他们的依赖,而复杂的公式,占其他行星产生的扰动,月亮,和许多其他因素。甚至一些模型考虑相对论效应中描述这个答案。
一个轨道的解决方案,我经常看到在论文中引用和使用处理气候迫使LA93,计算拉斯卡尔et al .(1993)和描述岁差的轨道,地球日晒量从-20年最高产量研究+ 10最高产量研究。然而,更准确的解决方案是可用的现在,更喜欢LA2004和特别是LA2010。
LA93的代码(在fortran)和输出的模型是可用的维齐尔目录,你可以查询模型输出,或进入这个FTP文件夹和下载的文件代码和长表与模型输出每一千年。
如果“三米兰柯维奇因素”你的意思是倾斜,离心率和旋进指数(偏心$ \ * $罪(近日点经度))。数据已经准备好这两个文件下载和使用:
CLIVAR0N。ASC(-20最高产量研究0最高产量研究)
CLIVAR0P。ASC(0 + 10最高产量研究最高产量研究)
的byte-per-byte描述这些文件中描述自述是
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -单位标签字节格式解释- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 4 - 11 F8.0 1000年t时间从1000年来J2000 13 - 31 D19.12——e偏心33 - 51 D19.12 rad eps倾斜(弧度)53 - 71 D19.12 - - - CP e *罪(经度的近日点移动equinox) - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -与代码还可以计算常见的气候代理如隔离在65°N或任何其他纬度。
LA2004(-50 + 20最高产量研究最高产量研究)研究团队实现了一个相同在线计算器,允许计算日晒量来源于轨道和岁差的数量对任何感兴趣的纬度。
LA2010解决方案描述La2010:新地球的轨道运动的长期解决方案250年,最高产量研究。你可以找到类似的数据文件和代码为这个解决方案天文地球古气候解决方案页面。我的理解,LA2010轨道参数代表了最先进的解决方案,它包括相对论修正。