In addition, a shadowing function applicable to multiple scattering is taken into account, which extends the calculation to the surface with anisotropic slope distribution.