Reason

From Joe Kelley at the NRC:

I am working on a 7-cell model that uses the GapHeatTransfer BC for every pair of sidesets as was done for the "reference" calculation in the extended scope task.

In that calculation, there were only two axial regions in the meshing-tools script, the support plate and everything else.

add_element!(mesh, elid, eltype, connectivity)
Add an element into the mesh. ´elid´ is the element id, ´eltype´ is the type of the element and ´connectivity´ is the connectivity of the element.

At least eltype needs to be defined more precisely. I guess it's a symbol :Seg2, :Tri3 and so on.