Model wave functions constructed from (1+1)D conformal field theory (CFT) have played a vital role in studying chiral topologically ordered systems. There usually exist multiple degenerate ground states when such states are placed on the torus. The common practice for dealing with this degeneracy within the CFT framework is to take a full correlator on the torus, which includes both holomorphic and antiholomorphic sectors, and decompose it into several conformal blocks. In this paper, we propose a pure chiral approach for the torus wave function construction. By utilizing the operator formalism, the wave functions are written as chiral correlators of holormorphic fields restricted to each individual topological sector. This method is not only conceptually much simpler, but also automatically provides us the anyon eigenbasis of the degenerate ground states (also known as the “minimally entangled states”). As concrete examples, we construct the full set of degenerate ground states for SO(n) ₁ and SU(n) ₁ chiral spin liquids on the torus, the former of which provide a complete wave function realization of Kitaev's sixteenfold way of anyon theories. We further characterize their topological orders by analytically computing the associated modular S and T matrices.