A numerical model, based on Parker’s transport equation, is utilized to investigate several aspects of the modulation of cosmic rays and the anomalous component beyond the solar wind termination shock (TS), in particular in the heliosheath, in the nose direction of the heliosphere. These aspects include: The latitude dependence of (1) radial perpendicular diffusion, the dominant diffusion process at the TS; (2) of the compression ratio, and (3) of the injection efficiency of pick-up ions. (4) The ‘break’ in the energy spectrum for anomalous protons. (5) The unfolding of the anomalous proton spectra in the heliosheath. (6) Effects on cosmic rays of adiabatic heating beyond the TS. Modelling results will be shown and discussed.