We model the dynamics of magnetization in an artificial analogue of spin ice specializing to the case of a honeycomb network of connected magnetic nanowires. The inherently dissipative dynamics is mediated by the emission and absorption of domain walls in the sites of the lattice, and their propagation in its links. These domain walls carry two natural units of magnetic charge, whereas sites of the lattice contain a unit magnetic charge. Magnetostatic Coulomb forces between these charges play a major role in the physics of the system, as does quenched disorder caused by imperfections of the lattice. We identify and describe different regimes of magnetization reversal in an applied magnetic field determined by the orientation of the applied field with respect to the initial magnetization. One of the regimes is characterized by magnetic avalanches with a 1/n distribution of lengths.