Труды КНЦ вып.9 (ГЕЛИОГЕОФИЗИКА вып. 5/2018(9))
отдельных нитей использовались атомарные операции. Это мало влияет на скорость вычислений, поскольку на одну операцию записи в массив приходится несколько тысяч алгебраических операций. Каждая вычислительная нить для своей линии поглощения сначала вычисляет интенсивность этой линии по формуле: Sa (T ) = Sa (Tref ) Qa (Tref ) exp ( -C 2 Ean / T)(1 - exp ( - C 2 E a / T)) Qa (T) exp ( - C2 Ean / Tref )(1 - exp ( -C 2 Eой / Tref )) (2) где Tr(g = 296K - нормальная температура, Qa(T ) - произведение вращательной и колебательной статистических сумм молекулы сорта a , EOT - энергия (в см"1) нижнего уровня перехода этой молекулы, E a. - энергия (в см-1) перехода, соответствующая линии поглощения с номером i этой молекулы, C 2 = hc / = 1.438769 см •K - вторая радиационная постоянная, в которой h - постоянная Планка, c - скорость света, kB - постоянная Больцмана. Далее вычисляются полуширина линии контура Лоренца и параметр доплеровской ширины линии по формулам: aL = 7 a self P a \ P r l V ref y + 7 a P_- P P r V ref УУ ref T v T У aD = 2RT Ua (3) где 7se-f - коэффициент самоуширения линии с номером i молекул сорта a , 70ir - коэффициент уширения воздухом этой линии, P ai - коэффициент температурной зависимости этой линии, P , - парциальное давление молекул сорта a , р = 1 атм, P - полное давление атмосферного газа, R - универсальная газовая постоянная, Ua - молярная масса молекул сорта a , Vai - частота центра линии поглощения молекул сорта (X . Рис. 1. Схематичное изображение контуров Фойгта (красные кривые). Черными стрелками изображены центры линий поглощения. После этого вычисляются значения контура Фойгта (рис. 1) в узлах сетки по частоте. Для расчета контура Фойгта используется функция, описанная в работе М. Кунтца [3] и модернизированная авторами данной работы для проведения параллельных расчетов на видеокарте. Контур Фойгта рассчитывается в диапазоне от 186
Made with FlippingBook
RkJQdWJsaXNoZXIy MTUzNzYz