fix boundary domains for quadrilateral meshes
fix output for quadrilateral meshes