Hi,

I have a domain and I want to solve the pde to get the displacement and phase field. meanwhile I want to set the value of phase field in 3 surfaces near the supports and load line as 1 (s=1)

I wrote this in julia using gridap in Drichlet boundary conditions. Is it a correct way or I should define it in another way?

```
order = 1
reffe_PF = ReferenceFE(lagrangian ,Float64,order)
V0_PF = TestFESpace(model ,reffe_PF;conformity =:H1,
dirichlet_tags = ["Left","Right", "Load","LeftSupport","RightSupport", "LoadLine"],
dirichlet_masks =[true,true,true,true,true,true])
sh = zero(V0_PF)
```

```
PF1(x)=1.0
PF2(x)=1.0
PF3(x)=1.0
PF4(x)=1.0
PF5(x)=1.0
PF6(x)=1.0
U_PF = TrialFESpace(V0_PF,[PF1 ,PF2 ,PF3,PF4 ,PF5 ,PF6])
```

Thank you