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

Но такое количество слоев в значительной мере увеличивает время вычислений, поэтому в реальности еще используется оптимизация (алгоритм представлен на рис. 2), которая заключается в усреднении параметров по процентному отношению массы вещества в каждом слое от общей массы столба воздуха. Пример значений плотности и температуры воздуха, полученных при помощи NRLMSISE-00 и оптимизированных при помощи вышеописанной обработки, представлены на рис. 3. Высота, км Высота, км Рис. 3. Значения температуры и плотности атмосферы Земли, полученные путем расчетов при помощи модели NRLMSISE-00 для географических координат 67 с.ш. и 33 в.д. с шагом по высоте составляет 0,1 км (а), а также аналогичные параметры, оптимизированные для вычислений. Хорошо видно, что при таком подходе заметно снижается количество используемых слоев, что положительно сказывается на скорости вычислений Fig. 3. The values of the temperature and density of the Earth's atmosphere, obtained by calculation using the NRLMSISE-00 model for geographic coordinates 67 N and 33 E with a step in height is 0.1 km (a), as well as similar parameters optimized for calculations. It is clearly seen that with this approach, the number of layers used is significantly reduced, which positively affects the speed of calculations Для сбора информации о вторичной компоненте KJI используются специальные детектирующие объемы, которые встроены в основные слои геометрии и являются не возмущающими (т.е., не влияют на параметры моделируемого трека частицы). При пролете границ таких чувствительных объемов интересующие параметры энергия Е, количество N, угол относительно направления полета первичной частицы cos(0)) записываются в массивы данных. Условно такой процесс можно представить графически (рис. 4.), либо описать при помощи следующего алгоритма: 1. Вначале задаются параметры и координаты детектирующего объема, происходит его расположение в необходимой области. 2. Во время инициализации определяются переменные HISTMAX, HIST MIN, которые являются границами диапазона энергий, в котором будет собираться вся необходимая информация о вторичных частицах. Ширину столбца при этом можно выразить при помощи формулы: 243