output: Gmap,E is also altered
Description: Set up a global numbering scheme for the whole mesh and return the salient information in the Gmap structure and the global element structure passed to the routine
