In the present work, the total magnetization in superconducting state is separated into critical state and paramagnetic components in terms of an H(x)-dependent magnetic flux density. Utilizing this model, we reproduce successfully M-H curves measured by Sandu et al. [Phys. Rev. B 74 (2006) 184511] and Sandu et al. [J. Supercond. Incorp. Novel Magn. 17 (2004) 701] for different forms of J(c).