docs: def mass_tot(self, rho0, Rs):
          """
          total mass within the profile
          :param rho0:
          :param a:
          :param s:
          :return:
          """
          m_tot = 2*np.pi*rho0*Rs**3
          return m_tot