Труды КНЦ вып.5. ГЕЛИОГЕОФИЗИКА. 8/2019(10)

излучения атмосферы в рамках приближения плоской и горизонтально однородной атмосферы. При этом численно решается одномерное по пространству уравнение переноса собственного излучения с помощью разработанной авторами модификации метода дискретных ординат, которая подробно описана в [24]. Для записи этого уравнения обозначим через и - косинус угла между направлением импульса фотона и вертикальным направлением. Этот угол будем называть зенитным. Иногда его отсчитывают от направления вниз. Через z обозначим высоту над поверхностью Земли в вертикальном столбе атмосферы, в котором производится расчёт поля излучения. Также для этого столба обозначим через T( z ) - температуру атмосферного газа на высоте z , через I ( v , z , u ) - интенсивность излучения с частотой V и зенитным углом, косинус которого равен и на высоте z , через er(v,z) и со(v ,z ) - коэффициенты объемного ослабления (экстинции), и альбедо 1-кратного рассеяния атмосферного газа на высоте z для излучения с частотой ѵ, через В ( у , Т ) = 2 Иѵъс 2 (ехр(/? ѵ /(к І: Г )) - 1 j В(ѵ,Т) - функцию Планка. Уравнение переноса собственного излучения атмосферы в рассматриваемом столбе можно записать в виде и dl(vz,u ) = _ I { y z u ^+ ( ^ _ C0 z))b(v, T{z)) + <j{y,z) d z N CD d w , / ч 1 /”2 ж 63 V,Z J ' I ( v , z ,w ) \ f x (v ,z ,m(w ,u ,(p ) )d (p 4 n j где w и и - косинусы зенитных углов фотонов до и после рассеяния, ср - разность между азимутальными углами фотонов до рассеяния и после рассеяния, m(w,u,<p) = w u + cos<p^j(l - w2) ( l - w2) - косинус угла рассеяния, a %( v , z ,m} - индикатрисса рассеяния для излучения с частотой ѵ на высоте z на угол, косинус которого равен м . Для уравнения (1) используются следующие граничные условия. На верхней границе для собственного излучения атмосферы является равенство нулю направленного вниз излучения: / ( ^ Zmax,W< 0) = 0 С2) На нижней границе направленное вверх излучение складывается из рассеянного поверхностью падающего излучения и из теплового излучения поверхности с температурой Тр . При изотропном рассеянии поверхностью условие на нижней границе имеет вид 226