Added broken tests to BoundaryDiscreteModelsTests.jl

The current implementation of BoundaryDiscreteModelsTests.jl
does not allow for a different ordering of the cells than the one
in the model from which it is created

test/GeometryTests/BoundaryDiscreteModelsTests.jl

@@ -60,4 +60,16 @@ n_∂Γ = get_normal_vector(∂Γ)
 @test ∑(∫(n_Γ⋅VectorValue(0.0,0.0,1.0))d∂Γ) ≈ 3.0
 @test ∑(∫(n_∂Γ⋅VectorValue(0.0,0.0,1.0))d∂Γ) ≈ -1.0
 
+domain   = (0,1,0,1,0,1)
+cells    = (1,1,1)
+bgmodel  = CartesianDiscreteModel(domain,cells)
+model    = BoundaryDiscreteModel(Polytope{2},bgmodel,Int32[3,4,5,6,1,2])
+topo     = get_grid_topology(model)
+cell_vertices = get_faces(topo,2,0)
+@test_broken cell_vertices[3] == [1,2,3,4]
+@test_broken cell_vertices[4] == [5,6,7,8]
+@test_broken cell_vertices[5] == [1,2,5,6]
+@test_broken cell_vertices[6] == [3,4,7,8]
+@test_broken cell_vertices[1] == [1,3,5,7]
+@test_broken cell_vertices[2] == [2,4,6,8]
 end