phn / jdcal Goto Github PK

>>> jd2gcal(0, -68964)
(-4901, 2, 1, 0.5)
>>> jd2gcal(0, -68965)
(-4901, 2, 0, 0.5)
>>> jd2gcal(0, -68966)
(-4901, 2, -1, 0.5)

>>> jd2gcal(0, -100010)
(-4986, 2, 1, 0.5)
>>> jd2gcal(0, -100011)
(-4986, 2, 0, 0.5)
>>> jd2gcal(0, -100012)
(-4986, 2, -1, 0.5)
>>> jd2gcal(0, -100041)
(-4986, 2, -30, 0.5)
>>> jd2gcal(0, -100042)
(-4986, 1, -1, 0.5)

$ pip install jdcal==1.0
Collecting jdcal==1.0
  Could not find a version that satisfies the requirement jdcal==1.0 (from versions: 1.0.1)
No matching distribution found for jdcal==1.0

$ pip install jdcal
$ ./test_jdcal.py
 passed

$ cat test_jdcal.py
#!/usr/bin/env python
# -*- coding: utf-8 -*-

from __future__ import print_function
from jdcal import jd2gcal

def test():
    for mjd in range(146097):
        w = {jd2gcal: jd2gcal(2400000.5, mjd)[:3], mjd2date: mjd2date(mjd)}
        if w['jd2gcal']!=w['mjd2date']:
            print(w['jd2gcal'], "!=", w['mjd2date'])
            return
    print(' passed')

# conversion from modified Julian date to Gregorian calendar date 
def mjd2date(mjd):
    n = mjd + 678881  # Gregorian
    a = 4*n + 3
    a += 4*((3*((4*(n+1))//146097+1))//4)  # Gregorian
    b = 5*((a%1461)//4) + 2
    year, month, day = a//1461, b//153 + 3, (b%153)//5 + 1
    if month > 12:
        year, month = year+1, month%12
    return year, month, day

test()