fix virial radius definition
fix return value in rg, reported by Wilma
add bug