Computational accuracy of radiation emerging from the ocean surface in the model atmosphere-ocean system

Computational accuracy of the Cox-Munk ocean surface model is discussed using the energy conservation of the incident radiation on it. The computational scheme in the atmosphere-ocean system is based on the doubling-adding method. Firstly, with the aid of the simple modification of the refraction matrix, computational accuracy is significantly improved even in the small discrete system of (15, 24), where the first and second numerals in the parentheses donote the discrete coordinates in the atmosphere and the ocean, respectively. This modification is necessary especially for the computations of the radiation including polarization effect which require a great deal of computer time and memory space. Secondly, the effects of shadows, multiple reflections, and refractions by the ocean surface are discussed to improve the accuracy. Finally, a computational example of albedo, the upwelling radiance and the degree of linear polarization at the top of the atmosphere-ocean system are discussed at wavelength 0.5μm.