Within the framework of further development of unified Monte Carlo code for the needs of biomedical optics and biophotonics, we present an approach for modeling of coherent polarized light propagation in highly scattering turbid media, such as biological tissues. The temporal coherence of light, linear and circular polarization, interference, and the helicity flip of circularly polarized light due to reflection at the medium boundary and/or backscattering events are taken into account. To achieve higher accuracy in the results and to speed up the modeling, the implementation of the code utilizes parallel computing on NVIDIA graphics processing units using Compute Unified Device Architecture. The results of the simulation of coherent linearly and circularly polarized light are presented in comparison with the results of known theoretical studies and the results of alternative modelings.