[Wien] total charge density

[Yundi Quan]

Hi,
According to userguide, the case.clmsum file contains the charge density
in harmonics representation. For l=0, m=0, charge density is given as
4*pi*r^2*\rho(r). For others, charge density is given as r^2*\rho(r). My
understanding is that the total charge density should be the sum of all
the harmonics. Therefore, I summed up all the lms ( I multiplied 4pi with
higher harmonics before summation). Is this a correct way of finding out
the radial charge density within the Muffin-tin sphere?
Another question is how is total charge (CTO) within the Muffin-tin sphere
calculated in scf file.

Thanks a lot.