mk3DP1_bcond(bc_list);

Parameters:

bc_list - a list of the form [u_bc, v_bc, ..], with one entry for each dependent variable in the problem. Each of the u_bc have the form [[c1,r1], [c2,r2], [c3,r3], [c4,r4], [c5,r5], [c6,r6]], There is one pair [c, r] for each face of the block. The first member can be either D or N: c=N will give Neumann boundary conditions, c=D: Dirichlet. The second member r should be a numerical value which is made available to the solver as qBCval. The faces are numbered acccording to 1: y = ylo, 2: x=xhi, 3: y=yhi, 4: x=xlo, 5: z=zlo, 6: z=zhi.

• Supply boundary condition information for uniform 3D grids generated by mk3DP1_cube. The information is saved in the file bcond.scratch.
• This is included as a demonstration of the 3D element. A coupling to a serious 3D mesh generator is still missing.

> # make boundary conditions file (bcond.scratch);
> # supply a pair [c, r] for each face of the block:
> # c=N: Neumann BC. The numerical value r is available as qBCval.
> # c=D: Dirichlet BC. The numerical value r is available as qBCval.
> # The faces are numbered acccording to
# 1: y = ylo, 2: x=xhi, 3: y=yhi, 4: x=xlo, 5: z=zlo, 6: z=zhi;
> u_bc:=[[N,0.],[N,0.],[N,0.],[N,0.],[D,0.],[N,1.]];
> v_bc:=[[D,0.],[D,0.],[D,0.],[D,0.],[D,0.],[N,1.]];
> mk3DP1_bcond( [u_bc, v_bc] );