Basic pyephem code

Question

I have a very simple code to get the longitude of the sun but when I compare the output to Astrolog and Astrodienst its incorrect, there is a 13 minute difference. I have not added Observer as I think default is midnight GMT (which is what I want). What am I doing wrong?

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()

Example of results for 2015/12/30:

code - 2015/12/30 00:00:00 277:43:36.6

Astrodienst - 7°56'39 Cap

Thanks

Albert · Accepted Answer · 2016-02-11 03:41:40Z

1

The reason why the 13 minute difference is because of the epoch setting, when I added

sun = ephem.Sun(start, epoch = start)

the results were the same as swiss ephemeris.

answered Feb 11, 2016 at 3:41

Albert

701 gold badge1 silver badge10 bronze badges

Sign up to request clarification or add additional context in comments.

1 Comment

Albert Over a year ago

I have only ticked this to close the question, as it drew no responses. Sorry for wasting peoples time.

Collectives™ on Stack Overflow

Basic pyephem code

1 Answer 1

1 Comment

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

1 Comment

Your Answer

Sign up or log in

Post as a guest

Related