Fix the SIA diffusive flux computation.

Set diffusivity to 0 in ghosted cells outside the domain (i,j < 0
and i >= grid.Mx, j >= grid.My). (See the diff for more.)