Just a caution - the formula applies to the true harmonic frequency ωe. Because of anharmonicity, the observed transition frequency ω0 may be somewhat different, and an apparent k derived from this frequency will vary slightly between isotopomers. Check whether your data are ω0 (observed) or ωe (derived by applying anharmonicity correction).