辐射过程是大气中最重要的物理过程之一,近年来计算机的发展,为大气辐射传输数值计算的迅速发展提供了可能。球谐函数法在20世纪初就被提出了,最早应用于星际辐射领域,在中子传输领域应用的最为广泛。但将其和其他方法相结合的历史却很短。球谐函数-离散坐标法是Evens在1998年提出的关于三维大气辐射传输计算的数值方法,这种方法结合了球谐函数法和离散坐标法的优点,相比传统的球谐函数法其物理意义更加清晰,数学模型更加简单,且有精度和计算效率上的优势。本文研究了球谐函数-离散坐标法(SHDOM)的基本原理,将SHDOM程序的计算结果与BMC程序的计算结果进行比较,分别计算了以下几种情况的辐射强度的变化:(1)各向同性散射下,固定波长和衰减系数,变化反照率;(2)各项同性散射下,固定波长和反照率,变化衰减系数;(3)各向异性散射下,相函数为(1 + cosθ)和(1 ? cosθ)时,分别变化反照率和衰减系数。应用云的典型伽马分布参数计算了层云的散射性质,并利用该散射性质,结合SHDOM程序,对二维大气辐射模型进行计算。根据一维算例的计算结果,随着计算模型中反照率和衰减系数的增大,程序的计算误差有变大的趋势。在云滴谱采用伽马分布典型参数的情况下,计算得到的衰减系数和反照率随有效半径的变化曲线会出现一个峰值,该峰值在衰减系数等于3附近出现。二维大气辐射传输算例中,顶面出射辐射强度的分布存在多个峰值,其成因有两个,一个是由于无云层衰减的底面直射产生,另一个是由对顶面向下辐射的强烈散射产生。