Non adiabatic coupling for internal conversion rate

I know I can calculate non adiabatic coupling for electronic part from Q-Chem . I am confused about to calculate for vibrational part of wave function.

Since under Franck condon approximation we can separate those two.

Can ezFCF help me in the calculation and if yes , how ?