A detailed algorithm is described that enables an implementation of a general valence bond (VB) method using the Clifford algebra unitary group approach (CAUGA). In particular, a convenient scheme for the generation and labeling of classical Rumer-Weyl basis (up to a phase) is formulated, and simple rules are given for the evaluation of matrix elements of unitary group generators, and thus of any spin-independent operator, in this basis. The case of both orthogonal and nonorthogonal atomic orbital bases is considered, so that the proposed algorithm can also be exploited in molecular orbital configuration interaction calculations, if desired, enabling a greater flexibility for N-electron basis-set truncation than is possible with the standard Gel'fandTsetlin basis. Finally, an exploitation of this formalism for the VB method, based on semiempirical Pariser-Parr-Pople (ppp)-type Hamiltonian and nonorthogonal overlap-enhanced atomic orbital basis, and its computer implementation, enabling us to carry out arbitrarily truncated or full VB calculations, is described in detail.