# not sure if this is correct

import datetime

def position(now=None):
   if now is None:
      now = datetime.datetime.now()

   days = (now - datetime.datetime(2001, 1, 1)) / datetime.timedelta(days=1)
   lunations = float.fromhex("0x1.a29b0e8eb915ep-3") + days * float.fromhex("0x1.1568432488c6ap-5")

   return round((lunations % 1) * 360)