This paper develops a Fortran code which is capable to construct the simplest LS eigcnfunctions for desired symmetry and determine all permitted atomic LS spectral terms under a given orbital occupancy by implementing and extending the Schaefer and Harris method. Examples (in some cases the most complete set to date) of multiple spectroscopic terms of LS coupling of atomic states for both non-equlvalent and equivalent electronic configurations arc given. It also corrects a few observed errors from the recent literature.