A numerical model of the cosmic ray storm in the two‐dimensional heliosphere is constructed incorporating the drift effect. We estimate the effect of a flare‐associated interplanetary shock and the disturbed region behind it (characterized by enhancement in velocity and magnetic field, and decrease in mean free path) on the density and anisotropy of cosmic rays in the heliosphere. As the disturbance propagates outward, a density enhancement appears on the front side, and a density depression region is produced on the rear side. The effect of drift on the cosmic ray storm appears most clearly in the higher‐latitude region. For the parallel (antiparallel) state of the solar magnetic field which corresponds to the pre‐(post‐) 1980 period, the density in the higher‐latitude region decreases (increases) before the shock arrival. The maximum density depression near the earth for the parallel state is greater than for the antiparallel state, and the energy spectrum of the density depression in percentage is softer for the parallel state than for the antiparallel state. Prior to the arrival of the shock, the phase of solar diurnal anisotropy begins to shift to the earlier hours, and its amplitude becomes greater for both polarity states. North‐south anisotropy also becomes greater because of the enhanced drift for both po