A variational approach for calculating Franck-Condon factors including mode-mode anharmonic coupling