wip
keep the same element order as given mesh if np==1
Merge branch 'fix-unidom'