A nonlinear fitting procedure is presented that employs all measured line positions and iteratively compares their values with those calculated from numerically diagonalized model Hamiltonians with adjustable molecular constants. Case (a) wavefunctions with definite parity are introduced as a convenient basis set, and the effects of spin-orbit, spin-spin, spin-rotation, and centrifugal distortion interactions neglected in the Born-Oppenheimer separation are included simultaneously using the Van Vleck transformation. The spectroscopic constants found by this procedure represent the minimum-variance, unbiased set and maintain, to a high degree of accuracy, the separate mechanical and magnetic meanings of the molecular constants. Arguments are presented that such spectroscopic constants with mechanical meaning allow the most accurate construction of potential energy functions using the Rydberg-Klein-Rees (RKR) procedure. Residual mechanical and magnetic ambiguities, such as engendered by Λ-type doubling, are discussed.