Square elements in a mesh

I would not use dolfin-xml, as it is not suited for parallel computing. I suggest using: Converter from GMSH to XDMF (with physical groups)