perl - libnova 北纬 89.5 度的奇怪行为

Question

我写了这个程序:

#!/usr/bin/perl 

use Astro::Nova qw(get_solar_equ_coords get_lunar_equ_coords get_hrz_from_equ 
                   get_solar_rst_horizon get_timet_from_julian 
                   get_julian_from_timet get_lunar_rst get_lunar_phase); 

$observer = Astro::Nova::LnLatPosn->new("lng"=>0,"lat"=>89.5); 

for $i (2456623..2456624) { 
  print "DAY: $i\n"; 
  ($status,$rst) = get_lunar_rst($i, $observer); 
  print "STATUS: $status\n"; 

  $rst->get_transit(); 

  $rise = $rst->get_rise(); 
  print "RISE: $rise\n"; 
  $set = $rst->get_set(); 
  print "SET: $set\n\n"; 
} 

并得到这些结果:

DAY: 2456623 
STATUS: 1 
RISE: 5.5664193588601e-309 
SET: 1.55254159695923e-267 

DAY: 2456624 
STATUS: -1 
RISE: 1.9634470382202e-153 
SET: 2.26294632209635e+137 

换句话说，月球从极地(始终向上)移动到极地下方地平线实际上从未升起或落下。

我意识到 89.5 度是一种极端情况，但为什么这段代码不应该工作吗？

正如 http://aa.usno.navy.mil/data/docs/RS_OneYear.php 所示的注释是此间隔内的月落 (JD 2456623 = 2013-11-26 12:00:00):

另一个奇怪的地方:如果我注释掉“$rst->get_transit();”线，我得到这些结果:

DAY: 2456623 
STATUS: 1 
RISE: 5.5664193588601e-309 
SET: 1.14372958360957e-268 

DAY: 2456624 
STATUS: -1 
RISE: 6.80740365931403e+199 
SET: 4.81766816905579e+151 

上升/下降仍然很奇怪(这很好，因为状态为 -1)，但是他们是不同的。我总是假设 get_lunar_rst() 返回一个固定的结构，但显然不是？检查交通的行为本身时间会改变结构吗？

编辑:好的，我运行了该程序两次，根本没有更改它并得到:

DAY: 2456623
STATUS: 1
RISE: 5.5664193588601e-309
SET: 2.99352717623831e-264

DAY: 2456624
STATUS: -1
RISE: 1.9634470382202e-153
SET: 2.26294632209635e+137

DAY: 2456623
STATUS: 1
RISE: 5.5664193588601e-309
SET: 3.04770606791278e-262

DAY: 2456624
STATUS: -1
RISE: 1.9634470382202e-153
SET: 2.26294632209635e+137

换句话说，2456623集合无缘无故地发生了变化。

Answer 1

我已联系 libnova 的作者，他们已确认这是一个错误并正在努力纠正它。

libnova 错误地认为如果天体在正北时位于地平线以上，则它们是绕极地的。这是不真实的:https://astronomy.stackexchange.com/q/963

我已经编写了一个修复程序，它不假设这一点，但如果物体的偏角是非单峰的，则可能仍然不准确:https://astronomy.stackexchange.com/questions/962/is-lunar-elevation-at-a-given-location-for-a-given-day-unimodal