This paper outlines a complete and self-consistent cell model theory of the electrokinetics of dense spherical colloidal suspensions for general electrolyte composition, frequency of applied field, zeta potential, and particle size. The standard electrokinetic equations, first introduced for any given particle configuration, are made tractable to computation by averaging over particle configurations. The focus of this paper is on the systematic development of suitable boundary conditions at the outer cell boundary obtained from global constraints on the suspension. The approach is discussed in relation to previously published boundary conditions that have often been introduced in an ad hoc manner. Results of a robust numerical calculation of high-frequency colloidal transport properties, such as dynamic mobility, using the present model are presented and compared with some existing dense suspension models.