We describe a rapid, accurate method for calculating rovibrational distributions in diatomic products from elementary chemical reactions. The basis of the model is momentum interconversion at a critical configuration defined in terms of molecular dimensions of the species involved. This approach shares common elements with recent models of inelastic processes and the kinematic reactive model of Elsum and Gordon. We point out that these and related approaches represent a development of Newtonian mechanics equivalent to that followed in the conventional formulation of classical mechanics, but one in which motive force for change at the molecular level is attributed to dp/dt rather than to dV/dq. This leads to a particularly transparent form of mechanics that uses only familiar data such as bond length, mass, spectroscopic constants, and velocity, yet may be applied to the highly resolved single collision experiments of molecular reaction dynamics. We describe key aspects of the computational method, e.g., the definition of the critical configuration, the disposal of reaction enthalpy, the manner of assigning product vibrational states, and the way in which conservation of energy is ensured. Examples are chosen to illustrate the range of reactions to which the method may be applied. Each would represent a challenge to conventional theory. We show that velocity-angular momentum diagrams may be used to interpret data and to give physical insight into the origins of observed rotational distributions. Good agreement is obtained between experimental and calculated (v,j) distributions for a wide range of elementary reactions suggesting that our model, despite its simplicity, captures the principal physics of chemical change at the molecular level.