Fix the SIA diffusive flux computation (in 0.5).

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.)