Calculating the frequencies and intensities of strongly anharmonic modes of adsorbates on surfaces: A low-cost but accurate computational approach