In the energy range of RHIC and LHC the mechanisms of nuclear suppression ofcharmonia are expected to be strikingly different from what is known for theenergy of the SPS. One cannot think any more of charmonium produced on a boundnucleon which then attenuates as it passes through the rest of the nucleus. Thecoherence length of charmonium production substantially exceeds the nuclearradius in the new energy range. Therefore the production amplitudes ondifferent nucleons, rather than the cross sections, add up and interfere, i.e.shadowing is at work. So far no theoretical tool has been available tocalculate nuclear effects for charmonium production in this energy regime. Wedevelop a light-cone Green function formalism which incorporates the effects ofthe coherence of the production amplitudes and of charmonium wave functionformation, and is the central result of this paper. We found a substantialdeviation from QCD factorization, namely, gluon shadowing is much stronger forcharmonium production than it is in DIS. We predict for nuclear effects $x_2$scaling which is violated at lower energies by initial state energy loss whichmust be also included in order to compare with available data. In this paperonly the indirect J/Psi originating from decay of P-wave charmonia areconsidered. The calculated x_F-dependence of J/Psi nuclear suppression is in agood accord with data. We predict a dramatic variation of nuclear suppressionwith x_F in pA and a peculiar peak at x_F=0 in AA collisions at RHIC.