Added adjustable BC to jacobi iter. function, cons. of mass not working when GAMMA>0