code: def datetime_to_djd(time):
      
  
      if time.tzinfo is None:
          time_utc = pytz.utc.localize(time)
      else:
          time_utc = time.astimezone(pytz.utc)
  
      djd_start = pytz.utc.localize(dt.datetime(1899, 12, 31, 12))
      djd = (time_utc - djd_start).total_seconds() * 1.0/(60 * 60 * 24)
  
      return djd
docs: Converts a datetime to the Dublin Julian Day
  
      Parameters
      ----------
      time : datetime.datetime
          time to convert
  
      Returns
      -------
      float
          fractional days since 12/31/1899+0000