Lat.: 12.4848° S Long.: 142.7866° E	Partial Solar Eclipse Magnitude: 0

 

 

Event	Date	Time (UT)	Alt	Azi
Start of partial eclipse (C1) :	2009/07/22	03:26:18.5	054.9°	338.9°
Maximum eclipse :	2009/07/22	03:27:32.1	054.8°	338.4°
End of partial eclipse (C4) :	2009/07/22	03:28:52.8	054.7°	337.9°

http://eclipse.gsfc.nasa.gov/5MCSE/5MCSEplate.html

 

第477个文件，包含2009年7月22日的日食食界图。

http://eclipse.gsfc.nasa.gov/5MCSE/plate477.pdf

 

看来NASA的日月食精度也是一些问题。

 

NASA的五千年日食，那个txt的5000年日食表，是把地球看成正球形计算的。我计算过，当中心线进入地球时，使用正球计算，中心点的结果与之相同，使用椭球，结果相差较多。

 

 

而那个尖角问题，也许是算法不同造成的。不过，我觉得，在那个尖角区域是看不到日食的。 

 

应当说，使用椭球比正球的精度要高一些。

毕竟赤道半径比极半径要大20多公里。

 

但如果许兄计算的结果与NASA的结果相差几十甚至上百公里，那么就不仅仅是椭球与正球的问题了。

可能还有月亮直径大小的问题。地月距离的问题。

 

对于半影的边缘，我觉得不必过于认真。也不见得一定要与NASA的结果保持绝对的一致。

 

*******************

我用《寿星》4。75版计算了一下：

经：-142.7866 纬：-12.4848

时间表 (日偏食)
初亏: 2009-07-22 11:26:09 UTC
食甚: 2009-07-22 11:27:32 UTC
复圆: 2009-07-22 11:29:01 UTC
时长: 0分0.0秒
食分: 0.00014
月日视径比: 1.07732(全或环食分)
是否NASA径比(1是,0否): 0
食分指日面直径被遮比例

 

上面的计算结果（NASA为否，不打勾）表明是可以看到（应当说是可以计算到，否则这么小的食分，肉眼根本看不出来）日偏食的。

若NASA打勾，则无偏食。

 

查了一下地图（纸质印刷的），经：-142.7866 纬：-12.4848。这个位置，至少在那个尖角以南约2度（即约220公里）。所以半影南界确实是进入澳洲那个尖角很多的。不应当是仅仅与尖角相擦而过。

已发布v4.76版：

 

1、改正“日月食”中，“莆田、北京”等几个标记点地标移位的问题，早期版没有这个问题，后来因为图形居中对齐出问题。

2、改正南北界线计算的错误

 

一直认为自已没算错，结果还是算错了。当时为了编程方便，没有考虑地球地转的影响，后来又加入了自转的影响，但影速与自转速度叠加时又出现正负号的错误。

 

NASA日月食的精度不高，中国天文年历的可能更高一些。

 

高精度不一定有很大的实际价值，但计算时提高精度，有助于发现算法上的不足。我当时加入NASA径比选项的原因，就是为了与NASA比对数据，以便发现计算过程中的错误，而不是为了以NASA的标准严格输出数据报表。

 

我最近，设计程序比较马虎，老出错。我越来越觉得，我不适合编写程序。

 

日月食计算的原理并不复杂，但实际计算却很麻烦。

就象解二次方程那样，原理原自配方法，而实际求解时，却要分几种情况：“单根、无根、两根”等情况。而日食可能出现的情况，以及相应的方程更多。

许兄不必过谦。

 

《寿星》程序，是许兄一个人单兵作战完成的。编程中遇到一些问题，是可以理解的。

网络上大多数收费的天文软件，都是一个团队共同开发的。

参与编程的人比较多，互相探讨的人也多。遇到问题也容易被及时修正。

 

虽然《寿星》是许兄一个人编写的，在我们很多朋友的使用过程中，若发现有小缺陷，我们都会及时提醒许兄。

从这点讲，《寿星》也是经过长期测试的，程序的主体不存在大问题的。《寿星》也是经受得起考验的。

 

《寿星》是功能很全的程序，使用者对其中的功能各有所取。只要这套程序能给使用者带来便利，我们就可以说《寿星》是成功的。

这个问题，我一开始以为是“地图没画好”，所以才造成几个地标偏移，一直也没提出疑问。

因为我想等着把地图重新画一遍。

 

后来许兄说可以借用GOOGLE的功能提取出地图轮廓。我就没去专门画地图了。毕竟用手工做地图，实在太费时间。

南极洲的地图已经非常简单了，我也花了五个小时以上。

对于亚欧大陆，估计没二十个小时，是无法完成的。

再把非洲，南北美洲，澳洲都画出来，估计要连续做一个星期，每天做十个小时，才可能完成。

但就算这样，也只是得到轮廓而已，比起GOOGLE的地图，差得太远。

GOOGLE的地图经过放大后，可以查到某条路上的一个小卖部，这么详细的地图，不是个人能完成的。

等我有空把《中国天文年历》2009版的日食图扫描下来上传，看看紫金山天文台的计算结果与NASA之间有没差别。

已发布V4.77版

 

1、“日食概略”中增加0.5食分界线

2、重新整理一下“日食概略”部分的内核代码

谢谢楼主了. 辛苦了.