code: def dodirot(D, I, Dbar, Ibar):
      
      d, irot = dogeo(D, I, Dbar, 90. - Ibar)
      drot = d - 180.
      if drot < 360.:
          drot = drot + 360.
      if drot > 360.:
          drot = drot - 360.
      return drot, irot
docs: Rotate a direction (declination, inclination) by the difference between
      dec=0 and inc = 90 and the provided desired mean direction
  
      Parameters
      ----------
      D : declination to be rotated
      I : inclination to be rotated
      Dbar : declination of desired mean
      Ibar : inclination of desired mean
  
      Returns
      ----------
      drot, irot : rotated declination and inclination