Calculation of temperature‐dependent multimode resonance Raman line shapes for harmonic potential surfacesa)