Calculation of two-dimensional vibrational potential energy surfaces utilizing prediagonalized basis sets and Van Vleck perturbation methods