Streamlined main script, minor changes to solver. Need to check equilibrium calculation. Also need to reformat solution loop to include possible rate interpolation.