Abstract: The edge of fusion plasma confinement devices is under intense scrutiny. Workers wishing to understand the flow across boundary regions most often take the challenge of employing the Braginskii collisional fluid model. The validity of this model is lost in most regimes of modern day tokamaks, but the Braginskii framework offers useful qualitative insight, and as recent Alcator C-Mod comparisons suggest, even quantitatively relevant answers. A global numerical solution of these equations, however, is difficult. In this work we develop a stable discrete solution of this model for global studies with a fair portrait of the small guys living on the edge, the electrons. We describe some numerical challenges and gains from its implementation in the 3D electromagnetic GDB code. The importance of some global effects will be motivated with simulations of a simple magnetized torus. These global features also arise in tokamaks and play an important role in symmetry breaking of C-Mod simulations. The accuracy of the GDB model is evaluated by comparing specific C-Mod experiments, whose profiles are well reproduced at low fields.