On the calculation of infrared intensities in solution within the polarizable continuum model