  In this paper we discuss the inclusive $J/\psi$ production in proton-proton
collisions from fusion of three pomerons. We demonstrate that this mechanism
gets dominant contribution from the region which can be theoretically described
by CGC/Saturation approach. Numerically, it gives a substantial contribution to
the $J/\psi$ production, and is able to describe the experimentally observable
shapes of the rapidity, momenta and multiplicity distributions. The latter fact
provides a natural explanation of the experimentally observed enhancement of
multiplicity distribution in $J/\psi$ production.
