Fortran code-lattice boltzmann-boundary condition

Hi everyone
Anybody can help me to write boundary condition fortran code.
I have boundary between cylinder and neighboring lattice point.
rw is the point at the cylinder surface , rp is the neighboring point inside the cylinder , rl is the neighboring site in the fluid area.
Yu-Mei-Lue-shyy model:
this method uses following equation: