Now using 18 energy groups.
Updated GPQ interval finder for integration limits. Now use an auxillary gaussian to find the high energy bound. With this method less roots (32 instead of 128). Have also updated the readrates routine to function with only the approximation (no shell model rates), and have added a line rejecting any calculations for nuclei that have number_densities(i)=0.0d0."
"
Applied a few cosmetic changes to the code structure.