Added preliminary get_density function
Defined constants used in functions