code: def ecliptic_latlon(self, epoch=None):
          
          vector = self.ecliptic_xyz(epoch)
          d, lat, lon = to_polar(vector.au)
          return (Angle(radians=lat, signed=True),
                  Angle(radians=lon),
                  Distance(au=d))
docs: Compute J2000 ecliptic coordinates (lat, lon, distance)
  
          If you instead want the coordinates referenced to the dynamical
          system defined by the Earth's true equator and equinox, provide
          an epoch time.