Моделирование физико-химических процессов в полярной ионосфере / Акад. наук СССР, Кол. фил. им. С. М. Кирова, Поляр. геофиз. ин-т. – Апатиты : Кольский филиал АН СССР, 1986.

A/tX* <0 re го г* і,ч /V, cm 3 І0Ъ І0г io' i -i ІО Рис.I. Распределениеконцентра цииэлектроноввдольконвективной траекториивэкваториальнойплоскости приL_: ----3 ; ------4 ; ----- 4,9; L_ - параметрМак-Илвейнадляточки траекториивOhLT. S и Рис.2. Распределениеконцент рацииэлектроновввечернемсекто реподанным работы/ 6 /. ехр(- bh» иі уравнение (3) сводится квиду: Ж (АЖ ) - * 00 +-Q “°* (4) Потокнаверхнейграницеопределяетсякак Ф - -D sin I ( f + jq) - -D Bin (5) Задаваяразностныесеткиповысотеипотраекторииконвекции{hj} , і = 1 , 2, іо; {tjj} , k = 0, I, 2, ... ко, можнорешитьуравнение (4) методом прогонки. Значениефункции U наверхнейграницерассчитываетсяизвыражений ( 2 ), ( 5 ) ипрогоночногосоотношения, записанногодляверхнейграницыhr=h:-0 ! ^ 0-1 (6 ) Умножая ( 6 ) на<рк0 , получим яІо-і -*5o-i*£o+£o-i* гдеN ^o_1 - і£ 0 -Г* 1 о* к з- к іо-1 дк,(рк •^іо Мо* Аппроксимируяввыражениях (5) и (2) производныеповременииповысотераз ностнымисоотношениями a<*loL V ) 6 т к к-1 _ Z -Z __ дт • ѵ ж I іо N io"H io-1 к _к /укчѴ — a R ... . ' УeNio(L ' • иисключаяизсистемы ( 6 ) и ( 2 ) величины ин^()_1, получимдвухточечное разностноеуравнениепоиндексу квдольтраекторииконвекции: -с^ук + с*ук- 1 + С* О; к * О) 1f»•»кО) (7) 8