fixed baryon conservation equation
analytical jacobian nearly converges
baryon conservation still non-convergent