output:
Description: As an alternative way to determine the local elemental to global numbering scheme than using the B rightarrowedge and B rightarrowface cummalative list the bmap array is defined in this function. The bma is an integer array of dimension nel by Nbmodes (number of boundary modes) which contain the global id of all the boundary modes within any element.