  Topological superconductors, as one of the most important research directions
at present, have attracted much attention because of their potential to realize
topological quantum computation. However, a universal computational tool based
on first-principle calculations for topological superconductivity was not yet
fully developed, and eventually significant challenges in predicting
topological superconducting materials. It is difficult to calculate the
topological superconducting properties of the system in a self-consistent
manner. In this paper, we develop a numerical method to characterize the
superconducting band spectrum and superconducting topological invariants of
two-dimensional (2D) slab system from first-principles calculations and
implemented in an open-source software WannierTools. We hope that it would
accelerate the discovery of the topological superconductor candidates.
