Evaluation of the internal radiation field in a composite planetary atmosphere with the lambert surface calculated by matrix method

If the atmosphere is simulated by a number of homogeneous sublayers, it was shown (Takashima, 1973a) that the intensity and polarization parameters emerging from any boundary of internal sublayer's field can be determined, provided that the diffuse reflection and transmission matrices of radiation of sublayers are known. Furthermore, if the surface (ground) is assumed to reflect light in accordance with the Lambert law, it is shown in this paper that the internal radiation field at boundary of any sublayer can be also computed in terms of the diffuse radiation matrices of sublayers rather than in terms of that of the entire atmosphere (Sekera, unpublished). The effect of polarization is included.