add boundary conditions when storing advected velocities in target field
minor changes to MACGrid::computeVorticityConfinement()
minor changes to MACGrid::computeBouyancy()