cormullion / astro.jl Goto Github PK

using Astro, Dates

d = DateTime(2015,3,20, 22,39,0):Minute(1):DateTime(2015,3,22,19,49,0)

for day in d
    eq = equation_time(Dates.datetime2julian(day)) # in fractions of day
    println(day, " ", eq)
end

minutes -7.433569 eq return -0.032435 asc 0.000007 2015-03-20T22:40:00
minutes 1432.566635 eq return 6.250751 asc 0.000018 2015-03-20T22:41:00
minutes 1432.566839 eq return 6.250752 asc 0.000029 2015-03-20T22:42:00
...
minutes 1433.123109 eq return 6.253179 asc 0.029916 2015-03-22T19:47:00
minutes 1433.123316 eq return 6.253180 asc 0.029927 2015-03-22T19:48:00
minutes -6.876478 eq return -0.030004 asc 0.029938 2015-03-22T19:49:00
minutes -6.876271 eq return -0.030003

Some rounding error or something.