In this paper, we consider a weakly coupled system of convection–diffusion equations subject to Robin boundary conditions, and having boundary and interior layers. The diffusion term of each equation is multiplied by a small singular perturbation parameter, but these parameters are assumed to be different in magnitude, and the source term is having a discontinuity at a point in the interior of the domain. An upwind scheme is used for the considered problem in conjunction with piecewise uniform Shishkin mesh. It is proved that the numerical approximations produced by this method are almost first order uniformly convergent with respect to both small parameters. Numerical results are presented to validate the theoretical results.