Add support for Gmsh high-order output

This MR primarily adds support for Gmsh output at arbitrarily high-order for quads, hexes, triangles and tetrahedra, and up to order 4 output for prisms.