Numerical analysis of the nonlinear rarefied gas flow caused by temperature gradient in a direction tangential to the wall through a planar channel of finite length is carried out based on the S- model kinetic equation under a range of Knudsen number from free-molecular to slip regime. Cercignani- Lampis solid boundary condition is chosen at channel walls. An implicit scheme is applied; the algorithm is optimized via massive parallelization in both physical and velocity spaces. Results showed that the gas-surface scattering model plays distinguished role for highly rarefied conditions.