Add tags to faces of 3D models using CartesianDiscreteModel()

[sy-nguyen-van]

Hello everyone,
I am new to Gridap.jl. I am creating a 3D model of a cantilever beam to apply a distributed load on the right face, and the left face is fixed.
Do you know how to add tags to LEFT and RIGHT faces ?
Here is my code:

using Gridap

L = 1 # Length of 3D cantilever beam
W = 0.2 # Width of 3D cantilever beam

domain = (0,L,0,W,0,W)
partition = (20,6,6)
model = CartesianDiscreteModel(domain,partition) # Mesh
writevtk(model,"model_3D")

add_tag_from_tags!(labels,"left",....) # How to add tag to LEFT face ?
add_tag_from_tags!(labels,"right",....) # How to add tag to RIGHT face ?

Ω = Triangulation(model)
Γ = Boundary(Ω,tags="right") # I want to a load on this face
reffe = ReferenceFE(lagrangian,Float64,1)
Vₕ = TestFESpace(Ω,reffe,dirichlet_tags="left") # This face is fixed

Thank you so much!

[JordiManyer]

just use vtk to visualize them, then select the ones you want. You should include edges and corners too, so 9 tags per side.

[sy-nguyen-van]

Thank you for your suggestions! It works now.
Here is my final code:

using Gridap
L = 1.0 # Length
W = 0.2 # Width
domain = (0,L, 0,W, 0,W)
partition = (40,12,12)
model = CartesianDiscreteModel(domain, partition)
labels = get_face_labeling(model)
# Tag the right face (x = L)
add_tag_from_tags!(labels,"left",[1,3,5,7,13,15,17,19,25]) # left face: corners (1,3,5,7); edges (13,15,17,19); others (25) 
add_tag_from_tags!(labels, "right", [2,4,6,8,14,16,18,20,26]) # Right face: corners (2,4,6,8); edges (14,16,18,20); others (26) 
writevtk(model,"3D_Problem")