In this paper, an improved numerical solver to evaluate the time-dependent radiative transfer equation (RTE) for underwater optical wireless communications (UOWC) is investigated. The RTE evaluates the optical path-loss of light wave in an underwater channel in terms of the inherent optical properties related to the environments, namely the absorption and scattering coefficients as well as the phase scattering function (PSF). The proposed numerical algorithm was improved based on the ones proposed in [1]-[4], by modifying the finite difference scheme proposed in [1] as well as an enhancement of the quadrature method proposed in [2] by involving a more accurate 7-points quadrature scheme in order to calculate the quadrature weight coefficients corresponding to the integral term of the RTE. Furthermore, the scattering angular discretization algorithm used in [3] and [4] was modified, based on which the receiver's field of view discretization was adapted correspondingly. Interestingly, the RTE solver has been applied to three volume scattering functions, namely: the single-term HG phase function, the two-term HG phase function [5], and the Fournier-Forand phase function [6], over Harbor-I and Harbor-II water types. Based on the normalized received power evaluated through the proposed algorithm, the bit error rate performance of the UOWC system is investigated in terms of system and channel parameters. The enhanced algorithm gives a tightly close performance to its Monte Carlo counterpart improved based on the simulations provided in [7], by adjusting the numerical cumulative distribution function computation method as well as optimizing the number of scattering angles. Matlab codes for the proposed RTE solver are presented in [8].