Treatment of incompressibility and boundary conditions in 3-D numerical spectral simulations of plane channel flows