  This article devises a new primal-dual weak Galerkin finite element method
for the convection-diffusion equation. Optimal order error estimates are
established for the primal-dual weak Galerkin approximations in various
discrete norms and the standard $L^2$ norms. A series of numerical experiments
are conducted and reported to verify the theoretical findings.
