added lattice pressure function
made n_b input to lattice pressure

this makes the simple solver converge but Ye is still off