<s> """<STR_LIT>""" <EOL> import numpy as np <EOL> MJD_0 = <NUM_LIT> <EOL> def ipart ( x ) : <EOL> """<STR_LIT>""" <EOL> return np . modf ( x ) [ <NUM_LIT:1> ] <EOL> def jcal2jd ( year , month , day ) : <EOL> """<STR_LIT>""" <EOL> year = int ( year ) <EOL> month = int ( month ) <EOL> day = int ( day ) <EOL> jd = <NUM_LIT> * year <EOL> x = ipart ( ( month - <NUM_LIT:9> ) / <NUM_LIT> ) <EOL> jd -= ipart ( ( <NUM_LIT:7> * ( year + <NUM_LIT> + x ) ) / <NUM_LIT> ) <EOL> jd += ipart ( ( <NUM_LIT> * month ) / <NUM_LIT> ) <EOL> jd += day <EOL> jd += <NUM_LIT> <EOL> return jd <EOL> def jd2jcal ( jd ) : <EOL> """<STR_LIT>""" <EOL> j = jd + <NUM_LIT> <EOL> k = ipart ( ( j - <NUM_LIT:1> ) / <NUM_LIT> ) <EOL> l = j - ( <NUM_LIT> * k ) <EOL> n = ipart ( ( l - <NUM_LIT:1> ) / <NUM_LIT> ) - ipart ( l / <NUM_LIT> ) <EOL> i = l - ( <NUM_LIT> * n ) + <NUM_LIT> <EOL> j = ipart ( ( <NUM_LIT> * i ) / <NUM_LIT> ) <EOL> day = i - ipart ( ( <NUM_LIT> * j ) / <NUM_LIT> ) <EOL> i = ipart ( j / <NUM_LIT> ) <EOL> month = j + <NUM_LIT:2> - ( <NUM_LIT> * i ) <EOL> year = ( <NUM_LIT:4> * k ) + n + i - <NUM_LIT> <EOL> return int ( year ) , int ( month ) , int ( day ) </s>
