A numerical scheme for computation of the convective diffusion of solute in saturated-unsaturated soil layers with arbitrarily. shaped boundaries has been developed applying the boundary, fitted coordinate method. Examples of computation on convective diffusion of solute, with surface and point or localized sources in an S-shaped hillslope are given. For a plain source like acid rain, the solute concentration can heardly observed in the throughflow discharge caused by the first rainfall, but it increases largely as soil moisture increases by the convection of rain water in unsaturated soils. For point sources, even in the initial stage of infiltration near the ground surface, the solute is diffused laterally, although the flow velocity is nearly vertically and the diffusing cloud of solute is convected downslope.