基本pyephem码

Question

我有一个非常简单的代码来获取太阳的经度，但是当我将输出与Astrolog和Astrodienst不正确的输出进行比较时，有13分钟的差异。我没有添加观察者，因为我认为默认是午夜格林尼治标准时间(这是我想要的)。我做错了什么？

import ephem

start = ephem.date('2015/01/01')
end = ephem.date('2015/12/31')
f2 = open("Sun", 'w')
while start <= end:
    sun = ephem.Sun(start)
    ecl = ephem.Ecliptic(sun)
    f2.write(str(ephem.date(start))+' '+ str(ecl.lon) +'\n')
    start+=1
f2.close()

2015/12/30年度成果实例：

编号- 2015/12/30 00:00:00 277:43:36.6

星体- 7°56'39帽

谢谢

Answer 1

13分钟差异的原因是由于时代的原因，当我添加

sun =ephem.Sun(开始，时代=开始)

结果与瑞士星历相同。