A single software framework is introduced to evaluate numerical accuracy of the A-grid (NICAM) versus C-grid (MPAS) shallow-water model solvers on icosahedral grids. The C-grid staggering scheme excels in numerical noise control and total energy conservation, which results in exceptional stability for long time integration. Its weakness lies in the lack of model error reduction with increasing resolution in specific test cases (especially the root-mean-square error). The A-grid method conserves well potential enstrophy and shows a linear reduction of error with increasing resolution. The gridpoint noise manifests itself clearly on A-grid, but much less on C-grid. We show that the Coriolis force term on C-grid has a larger error than on A-grid. To treat the Coriolis term and kinetic energy gradient on an equal footing on C-grid, we propose combining these two quantities into a single tendency term and computing its value by a linear combination operation. This modification alone reduces numerical errors but still fails to converge the maximum error with resolution. The method of Peixoto can solve the maximum-error nonconvergence problem on C-grid but degrades the numerical stability. For the steady-state thin-layer test (0.01 m in depth), the A-grid method is less susceptible than C-grid methods, which are presumably disrupted by the Hollingsworth instability. The effect of horizontal diffusion on model accuracy and energy conservation is shown in detail. Programming experience shows that software implementation and optimization can strongly influence computational performance for models, although memory requirement and computational load of the two schemes are comparable.