Abstract We consider numerical boundary conditions for high order finite difference schemes for solving convection-diffusion equations on arbitrary geometry. The two main difficulties for numerical boundary conditions in such situations are: (1) the wide stencil of the high order finite difference operator requires special treatment for a few ghost points near the boundary; (2) the physical bou...