Once the sites of anomalous scatterers are known, phase probability for the total scattering in the unit cell can be obtained. The output contains phase probability distribution that can be used for density modification.