gfoclt can be used to find eclipses (occultations that involve our Moon) and other occultations.

gfoclt
Abstract
Determine time intervals when an observer sees one target occulted by, or in transit across, another.
The surfaces of the target bodies may be represented by triaxial ellipsoids or by topographic data provided by DSK files.