Added spherical peak density.  Added colors on important quantities on the output.  Added possibility of calculating only column density and exiting before doing any fits.