“…This model is used due to its unique ability to calculate upwelling and downwelling irradiances, radiances, reflectance, mean cosines, shape factors, and directional attenuation coefficients in realistic natural waters and can be applied to water columns containing fluid mud lutocline layers. This model has been vigorously tested [1,3,4,5] and compared to selected results of previous Monte Carlo models reported by Gordon [9], Kattawar [15], Morel and Gentilli [19], Stavn [13], Kirk [14], etc. Comparisons of the models repeated in the literature are shown in Table 1.…”