Boundaries are now update directly in shift_int. Potential resulting from selfgravity is now also shifted by integer number of cells to (possibly) optimize history recycling. NOT TESTED