A comprehensive ring current model (CRCM) has been developed that couples the Rice Convection Model (RCM) and the kinetic model of Fok and coworkers. The coupled model is able to simulate, for the first time using a self-consistently calculated electric field, the evolution of an inner magnetosphere plasma distribution that conserves the first two adiabatic invariants. The traditional RCM calculates the ionospheric electric fields and currents consistent with a magnetospheric ion distribution that is assumed to be isotropic in pitch angle. The Fok model calculates the plasma distribution by solving the Boltzmann equation with specified electric and magnetic fields. To combine the RCM and the Fok model, the RCM Birkeland current algorithm has been generalized to arbitrary pitch angle distributions. Given a specification of height-integrated ionospheric conductance, the RCM component of the CRCM computes the ionospheric electric field and currents. The Fok model then advances the ring current plasma distribution using the electric field computed by the RCM and at the same time calculates losses along particle drift paths. We present the logic of CRCM and the first validation results following the H+ distribution during the previously studied magnetic storm of May 2, 1986. The H+ fluxes calculated by the coupled model agree very well with observations by AMPTE/CCE. In particular, the coupled model is able to reproduce the high H+ flux seen on the dayside at L~2.3 that the previous simulation, which employed a Stern-Volland convection model with shielding factor 2, failed to produce. Though the Stern-Volland and CRCM electric fields differ in several respects, the most notable difference is that the CRCM predicts strong electric fields near Earth in the storm main phase, particularly in the dusk-midnight quadrant. Thus the CRCM injects particles more deeply and more quickly. ¿ 2001 American Geophysical Union