LB Temperature Boundary Condition Implementation

Could anyone tell me how to implement Temperature for a boundary for say D2Q9?
If the pseudoparticles are tIn[0] to tIn[8] how do we interrelate them?