In [1]:
from pyciss import _utils as u

In [31]:
u.which_epi_janus_resonance('Janus', '2005-01-01')


Out[31]:
'Janus2'

In [7]:
name='Janus'

In [9]:
from astropy.time import Time

In [30]:
time = '2004-02-01'

t1 = Time('2002-01-21').to_datetime()
delta = Time(time).to_datetime() - t1
print(delta.days)
yearfraction = delta.days / 365
print(yearfraction)
if int(yearfraction / 4) % 2 == 0:
    print( name + '2')
else:
    print( name + '1')


741
2.03013698630137
Janus2

In [22]:
int(-1/4) % 2


Out[22]:
0

In [ ]:
def test_which_epi_janus_resonance():
    assert u.which_epi_janus_resonance('Janus', '2005-01-01') == "Janus2"
    assert u.which_epi_janus_resonance('Epimetheus', '2007-01-01') == "Epimetheus1"
    assert u.which_epi_janus_resonance('Janus', '2011-01-01') == "Janus2"