The complex variable reproducing kernel particle method (CVRKPM) of solving two-dimensional variable coefficient advection—diffusion problems is presented in this paper. The advantage of the CVRKPM is that the shape function of a two-dimensional problem is formed with a one-dimensional basis function. The Galerkin weak form is employed to obtain the discretized system equation, and the penalty method is used to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional variable coefficient advection—diffusion problems are obtained. Two numerical examples are given to show that the method in this paper has greater accuracy and computational efficiency than the conventional meshless method such as reproducing the kernel particle method (RKPM) and the elementfree Galerkin (EFG) method.