*return IDs of all triangles contained within and/or intersecting a circle centered on a given ra and dec*

        **Key Arguments:**
            - ``ra`` -- RA of central point in decimal degrees or sexagesimal
            - ``dec`` -- DEC of central point in decimal degrees or sexagesimal
            - ``radius`` -- radius of circle in degrees
            - ``inclusive`` -- include IDs of triangles that intersect the circle as well as those completely inclosed by the circle. Default *True*
            - ``convertCoordinates`` -- convert the corrdinates passed to intersect. Default *True*
            -

        **Return:**
            - ``trixelArray`` -- a numpy array of the match trixel IDs

        **Usage:**

            To return the trixels overlapping a circle with a 10 arcsec radius centred at 23:25:53.56, +26:54:23.9

            .. code-block:: python

                overlappingTrixels = mesh16.intersect(
                    ra="23:25:53.56",
                    dec="+26:54:23.9",
                    radius=10 / (60 * 60),
                    inclusive=True
                )

            Or to return the trixels completing enclosed by a circle with a 1 degree radius centred at 23:25:53.56, +26:54:23.9

            .. code-block:: python

                overlappingTrixels = mesh16.intersect(
                    ra="23:25:53.56",
                    dec="+26:54:23.9",
                    radius=1,
                    inclusive=False
                )

def intersect(self, ra, dec, radius, inclusive=True, convertCoordinates=True):
        """*return IDs of all triangles contained within and/or intersecting a circle centered on a given ra and dec*

        **Key Arguments:**
            - ``ra`` -- RA of central point in decimal degrees or sexagesimal
            - ``dec`` -- DEC of central point in decimal degrees or sexagesimal
            - ``radius`` -- radius of circle in degrees
            - ``inclusive`` -- include IDs of triangles that intersect the circle as well as those completely inclosed by the circle. Default *True*
            - ``convertCoordinates`` -- convert the corrdinates passed to intersect. Default *True*
            -

        **Return:**
            - ``trixelArray`` -- a numpy array of the match trixel IDs

        **Usage:**

            To return the trixels overlapping a circle with a 10 arcsec radius centred at 23:25:53.56, +26:54:23.9

            .. code-block:: python

                overlappingTrixels = mesh16.intersect(
                    ra="23:25:53.56",
                    dec="+26:54:23.9",
                    radius=10 / (60 * 60),
                    inclusive=True
                )

            Or to return the trixels completing enclosed by a circle with a 1 degree radius centred at 23:25:53.56, +26:54:23.9

            .. code-block:: python

                overlappingTrixels = mesh16.intersect(
                    ra="23:25:53.56",
                    dec="+26:54:23.9",
                    radius=1,
                    inclusive=False
                )
        """
        # CONVERT RA AND DEC DECIMAL DEGREES

        if convertCoordinates == True:
            converter = unit_conversion(
                log=self.log
            )
            ra = converter.ra_sexegesimal_to_decimal(
                ra=ra
            )
            dec = converter.dec_sexegesimal_to_decimal(
                dec=dec
            )

        if inclusive:
            inc = 1
        else:
            inc = 0
        return super(HTM, self).intersect(ra, dec, radius, inc)