fix boundary domains for quadrilateral meshes
fix output for quadrilateral meshes
_cell_closure hack