The regions where boundary conditions are applied are set through the mesh nodes defined by the layers. In the example of Fig. 1 there are 3 nodes in the x and z directions and 4 nodes in the y direction. To apply a boundary condition in the superior face of the domain (perpendicular to the z axis) the nodes are defined as shown in Table II.