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(a{u}_x{)}_x+O(h^4)$, yielding a fourth-order one-step residual when the forward time step is proportional to $h^2$.