基本的pyephem代码

Question

我有一个非常简单的代码来获取太阳的经度但是当我将输出与Astrolog和Astrodienst的比较不正确时，有13分钟的差异。我没有添加Observer，因为我认为默认是午夜GMT（这是我想要的）。我做错了什么？

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

Astrodienst - 7°56'39 Cap

由于

Answer 1

13分钟差异的原因是由于纪元设置，当我添加

sun = ephem.Sun（start，epoch = start）

结果与瑞士星历相同。