A time dependent coupled cluster approach to the calculation of Resonance Raman excitation profiles on general anharmonic surfaces is presented. The vibrational wave functions on the ground electronic surface are obtained by the coupled cluster method (CCM). It is shown that the propagation of the vibrational ground state on the upper surface is equivalent to propagation of the vacuum state by an effective hamiltonian generated by the similarity transformation of the vibrational hamiltonian of that surface by the CCM wave operator of the lower surface up to a normalization constant. This time propagation is carried out by the time-dependent coupled cluster method in a time dependent frame. Numerical studies are presented to asses the validity of the approach.