In most chemistry methods developed for the direct simulation Monte Carlo (DSMC) technique, chemical reactions are computed as an integral part of the collision simulation routine. In the macroscopic chemistry method developed here, the simulation of collisions and the simulation of reactions are decoupled; reactions are computed independently, after the collision routine. The number of reaction events to perform in each cell is calculated using the macroscopic reaction rates k+, k- and equilibrium constant K*, calculated from the local macroscopic flow conditions. The macroscopic method is developed here for the symmetrical diatomic dissociating gas. For each dissociation event, a single diatomic simulator particle is selected with a probability based on its internal energy, and is replaced by two atomic particles. For each recombination event, two atomic particles are selected at random, and are replaced by a single diatomic particle. The dissociation energy is accounted for by adjusting the translational thermal energies of all particles in the cell. The macroscopic method gives density profiles in agreement with experimental data for the chemical relaxation region downstream of a strong shock in nitrogen. In the non-equilibrium regions within the shock, and along the stagnation streamline of a blunt cylinder in rarefied flow, the macroscopic method gives results in excellent agreement with those obtained using the most common conventional DSMC chemistry method in which reactions are calculated during the collision routine. The number of particles per computational cell has a minimal effect on the results provided by the macroscopic method. Unlike most DSMC chemistry methods, the macroscopic method is not limited to simple forms of k+, k- and K*. Any forms may be used, and these may be any function of the macroscopic conditions. This is demonstrated by using a two-temperature rate model, and a form of K* with a number density dependence. With the two-temperature model, the macroscopic method gives densities in the post-shock chemical relaxation region that also agree with the experimental data. For a form of K^* with a number density dependence, the macroscopic method can accurately reproduce chemical recombination behavior. In a primarily dissociative flow, the number density dependence of K* has very little effect on the flow. The macroscopic method requires slightly less computing time than the most common DSMC chemistry method.