added volume filling fraction correction to baryon conservation
cleaned up code a bit, added some comments
cleaned up the code a bit more
solver converges with large table