docs: def julian_day(t: date) -> int:
      """Convert a Python datetime to a Julian day"""
      # Compute the number of days from January 1, 2000 to date t
      dt = t - julian_base_date
      # Add the julian base number to the number of days from the julian base date to date t
      return julian_base_number + dt.days