The article presents an implementation of the Frenkel exciton coupling model into the OpenMolcas program package. The implemented protocol enables calculations of the monomers' electronic structure on a multiconfigurational level and direct computations of the Coulomb intermonomer couplings, thus, avoiding costly supermolecular calculations followed by a diabatization scheme. Additional efficiency of the computational scheme is provided by the Cholesky decomposition of the two-electron integrals of the aggregate system.