Dust aerosols have a strong scattering peak at small scattering angles (near the sun). The usual single or double Henyey-Greenstein phase functions are convenient in general but do not peak enough to emulate that strong peak in the circumsolar region. Is there any other simple phase function that would provide better results?