Add advection mask, and correct x,y,z order for advection.
Set the populations for non-advected lattice sites too.
Mask density calculations.
Add mask to collision.
Add mask to momentum.
Add velocity mask.