Most of the programs use rectilinear normal coordinates in frequency calculations. The molecular vibrations defined by rectilinear normal coordinates are based on the infinitesimal approximation of the nuclear displacements. Therefore, such coordinates are not the useful for describing the molecular vibrations involving large displacements from the equilibrium. For such vibrations curvilinear coordinates system is best choice.
Let me know program that do vibrational calculations in curvilinear coordinates system?