Fixed a bug in the linear solvers - residual was initialized with rhs in FT, which could lead to adding up rhs multiple times in the residual during communication, rhs is now added in the LT instead.
Fixed gravity to some degree, saturation can still surpass 1 on the boundary and in ugly spots for some reason