// Compute julian date (must be used with UTC+00:00 Time Zone)
y=openyear
m=openmonth
d=openday
// The most negative year in the current Julian period is -4713 Gregorian (4714 BC), so we want to add something to y so we can simplify reasoning about our calculations by considering only positive year values. Any value divisible by 4000 could be chosen, since the pattern of leap years in the reformed Gregorian calendar repeats every 4000 years.
y=y+8000
// Since leap days are added at the end of February, it makes sense to consider January and February as the last two months (numbered 13 and 14) of the previous year. That way, leap days are always added at the end of the notional year.
if(m<3) then
y=y-1
m=m+12
endif
// The Julian calendar has normal years of 365 days, plus extra leap days every fourth year.
jd1 = (y*365) + (y - (y mod 4)) /4
// The Gregorian calendar adds the correction that a century year is not a leap year, unless it is also divisible by 400.
jd2 = -(y - (y mod 100)) /100 + (y - (y mod 400)) /400
// This correction allows for the 3287-year difference between the arbitrary positive year offset we added earlier and the real Julian Day epoch of 1st. January, 4713 BC, Julian calendar (24th. November, 4714 BC, Gregorian calendar).
jd3 = -1200820
// With the months starting with March=3 and ending with February=14, there is a simple pattern of month lengths: 31, 30, 31, 30, 31, 31, 30, 31, 30, 31, 31 That is, there are always 30=150/5 days in each month with an additional day every 3/5 month. Since we started with March=3, the first value given would be (3*153+3)/5 = 92, but we want to start from zero so we make a correction.
jd4 = ((m*153+3) - ((m*153+3) mod 5)) /5 - 92
// The Julian Day number goes up once a day. Days of the month are conventionally numbered from 1, so we subtract 1.
jd5 = +d-1
// julian day
jd = jd1 + jd2 + jd3 + jd4 + jd5
frac = opentime/100
rounded = round(frac)
if rounded > frac then
floorBase = (rounded-1) * 100
else
floorBase = rounded * 100
endif
openseconds = opentime - floorBase
secs = jd*86400 + openhour*3600 + 60*openminute + openseconds
return secs - secs[1] as "duration in seconds"
.
