Varying relaxation times and combinations of bounceback and particle reflections

I am trying to implement a varying relaxation time in Palabos for both BGK and MRT methods. Actually, I want to make my relaxation time a function of say the fluid density (Knudsen number). In addition, I want to implement a solid boundary condition that will be a (linear) combination of particle bounce back and reflection.
Please how do I go about this?