In this routine the unrefined tetrahedral mesh is catalogued. It consists of the folowing tasks:

Notice that the tetrahedral elements are stored as if all of them were linear, ie. C3D4-elements. The only quadratic information of the unrefined mesh is stored in fields iedgext and ifacext, to be discussed below.