Symmetric half-grid diffusion consistency (source code)

= Symmetric half-grid diffusion consistency

Evaluating a variable diffusion coefficient at both adjacent half-grid points makes odd spatial Taylor terms cancel. The conservative spatial bracket is $h^2(au_x)_x+O(h^4)$, yielding a fourth-order one-step residual when the forward time step is proportional to $h^2$.