add routine to compute scale radius from r_Delta
abandon delta_s, rho_s. Fixes #1
consistent function naming, pep8 clean-up