The linear stability of a non-divergent barotropic parallel shear flow in a zonal and a non-zonal channel on the β plane was examined numerically. When the channel is non-zonal, the governing equation is slightly modified from the Orr-Sommerfeld equation. Numerical solutions were obtained by solving the discretized linear perturbation equation as an eigenvalue problem of a matrix. When the channel is zonal and lateral viscosity is neglected the problem is reduced to the ordinary barotropic instability problem described by Kuo's (1949) equation. The discrepancy between the stability properties of westward and eastward flows, which have been indicated by earlier studies, was reconfirmed. It has also been suggested that the unstable modes are closely related to the continuous modes discretized by a finite differential approximation. When the channel is non-zonal, the properties of unstable modes were quite different from those of the zonal problem in that: (1) The phase speed of the unstable modes can exceed the maximum value of the basic flow speed; (2) The unstable modes are not accompanied by their conjugate mode; and (3) The basic flow without an inflection point can be unstable. The dispersion relation and the spatial structure of the unstable modes suggested that, irrespective of the orientation of the channel, they have close relation to the neutral modes (Rossby channel modes) which are the solutions in the absence of a basic shear flow. The features mentioned above are not dependent on whether or not the flow velocity at the boundary is zero.