Hi,
n2to3 should preserve the boundaryID of the original 2D mesh and extrude it to all elements on the face 1 to face 4.
Face 5 and 6 will use user's input and store BC in CBC.
In userdat2, you can do something like this
do iel=1,nelt
do ifc=1,4 ! convert bc id to cbc based on original 2D mesh
if (boundaryID(ifc,iel).eq.1) cbc(ifc,iel,1) = 'W '
enddo
do ifc=5,6 ! convert face 5-6 (z direction) from n2to3
if (cbc(ifc,iel,1).eq.'v ') cbc(ifc,iel,1) = 'W '
enddo
enddo
Personally, I like putting use CBC to carry my boundary conditions. And, I use 'W ', 'W01', 'W02', ... to act as a the boundaryID.
Therefore, if you want to distinguish the newly generated BC on face 5-6 and some other 'W ', you can put 'W01' in n2to3 and convert it in userdat2.
Hope this helps,
Yu-Hsiang
--