Cleanup september

Starting to do the utility calculations in the model class in numpy arrays for hopefully faster computation. Solved several bugs and improved readability and reduced duplicate code.