Added index of refraction calculation and accessor.

Added the calculation of the index of refraction as a function of
density and the accessor to the _refraction member.  Also fixed a few
const doubles in the helper methods.