Mass conservation

The mass conservation is guaranteed by

$$\frac{3}{4\pi}\frac{\partial M_{\rm r}}{\partial r^3}=\rho_{\star}+\rho_{\rm g}$$ (65)