fix units error in if(R<RBH) for enclosed mass