  Two sets of mutually commuting $q-$difference operators $x_i$ and $y_j$,
$i,j=1, ...,N$ such that $x_i$ and $y_i$ generate a homomorphic image of the
$q-$Onsager algebra for each $i$ are introduced. The common polynomial
eigenfunctions of each set are found to be entangled product of elementary
Pochhammer functions in $N$ variables and $N+3$ parameters. Under certain
conditions on the parameters, they form two `dual' bases of polynomials in $N$
variables. The action of each operator with respect to its dual basis is block
tridiagonal. The overlap coefficients between the two dual bases are expressed
as entangled products of $q-$Racah polynomials and satisfy an orthogonality
relation. The overlap coefficients between either one of these bases and the
multivariable monomial basis are also considered. One obtains in this case
entangled products of dual $q-$Krawtchouk polynomials. Finally, the `split'
basis in which the two families of operators act as block bidiagonal matrices
is also provided.
