A posteriori tests of large-eddy simulations for the temporal mixing latyer are performed using a variety of numerical methods in conjunction with the dynamic mixed subgrid model for the turbulent stress tensor. The results of the large-eddy simulations are compared with filtered direct numerical simulation (DNS) results. Five numerical methods are considered. The cell vertex scheme (A) is a weighted second-order central difference. The transverse weighting is shown to be necessary, since the standard second-otder central difference (A') gives rise to instabilities. By analogy, a new weighted fourth-order central difference (B) is constructed in order to overcome the instability in simulations with the standard fourth-order central method (B'). Furthermore, a spectral scheme (C) is tested. Simulations using these schemes have been performed for the case where the filter width equals the grid size (I) and the case where the filter width equals twice the grid size (11). The filtered DNS results are best approximated in case I1 for each of the numerical methods A, B and C. The deviations from the filtered DNS data are decomposed into modelling error effects and discretization error effects. In case I the absolute modelling error effects are smaller than in case I1 owing to the smaller filter width, whelleas the discretization error effects are larger, since the flow field contains more small-scale contributions. In case I scheme A is preferred over scheme B, whereas in case I1 the situation is the reverse. In both cases the spectral scheme C provides the most accurate results but at the expense of a considerably increased computational cost. For the prediction of some quantities the discretization errors are observed to eliminate the modelling errors to some extent and give rise to reduced total errors.