Gluon shadowing and unitarity effects

New data from HERA experiment on deep inelastic scattering have been used to parametrize nucleon and Pomeron structure functions. Within the Gribov theory, the parameterizations were employed to calculate gluon shadowing for various heavy ions. The latter was compared with predictions from other models. Calculations of multiplicity reduction due to gluon shadowing for d+Au collisions at forward rapidities at $\sqrt{s}$=200 GeV are in good agreement with BRAHMS data on the nuclear modification factor.