atmosnrlmsise00

Реализуйте математическое представление 2 001 Военно-морского Спектрометра Массы Научно-исследовательской лаборатории Соединенных Штатов и Несвязной Поля точек Радарной Экзосферы

Синтаксис

[T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds) [T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, localApparentSolarTime) [T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, f107Average, f107Daily, magneticIndex) [T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, flags) [T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, otype) [T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, action)

Описание

[T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds) реализует математическое представление 2 001 Военно-морского Спектрометра Массы Научно-исследовательской лаборатории Соединенных Штатов и Несвязной Поля точек Радарной Экзосферы (NRLMSISE-00) модели класса MSIS^®. NRLMSISE-00 вычисляет нейтральную атмосферу эмпирическая модель от поверхности, чтобы понизить экзосферу (от 0 м до 1 000 000 м). Опционально, это выполняет это вычисление включая вклады от аномального кислорода, который может влиять на спутниковое перетаскивание выше 500 000 м.

[T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, localApparentSolarTime) задает массив m локальное очевидное солнечное время (часы).

[T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, f107Average, f107Daily, magneticIndex) задает массивы m 81-дневное среднее значение потока F10.7 (сосредоточенный на doy), m-by-1 ежедневный поток F10.7 в течение предыдущего дня и m-by-7 магнитной информации об индексе.

[T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, flags) задает массив 23, чтобы включить или отключить конкретные изменения для выходных параметров.

[T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, otype) задает вектор символов или строку для общей массовой плотности вывод.

[T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, action) задает из области значений входное действие.

Эта функция имеет ограничения модели NRLMSISE-00. Для получения дополнительной информации см. документацию модели NRLMSISE-00.

Модель NRLMSISE-00 использует UTseconds, localApparentSolarTime и longitude независимо. Эти аргументы не имеют равного значения для каждой ситуации. Для наиболее физически реалистического вычисления выберите эти три переменные, чтобы быть сопоставимыми по умолчанию:

 localApparentSolarTime = UTseconds/3600 + longitude/15

При наличии можно включать отклонения от этого уравнения для localApparentSolarTime, но они имеют второстепенное значение.

Входные параметры

action

Действие для входа из области значений. Задайте тот:

'Error'

'Warning' (значение по умолчанию)

'None'

altitude

Массив m-by-1 высоты, в метрах.

dayOfYear

Массив m-by-1 года.

f107Average

Массив m-by-1 81-дневное среднее значение потока F10.7 (сосредоточенный в день года (dayOfYear)). Если вы specifyf107Average, необходимо также задать f107Daily и magneticIndex. Эффекты f107Average не являются большими или установленными ниже 80 000 м; поэтому, значение по умолчанию равняется 150.

Эти значения f107Average соответствуют радио-потоку на 10,7 см на фактическом расстоянии от Земли от Sun. Значения f107Average не соответствуют радио-потоку на уровне 1 а. е. Следующий сайт обеспечивает оба класса значений: ftp://ftp.ngdc.noaa.gov/STP/space-weather/solar-data/solar-features/solar-radio/noontime-flux/penticton/

Смотрите ограничения в Описании для получения дополнительной информации.

f107Daily

Массив m-by-1 ежедневный F10.7 течет в течение предыдущего дня. Если вы задаете f107Daily, необходимо также задать f107Average и magneticIndex. Эффекты f107Daily не являются большими или установленными ниже 80 000 м; поэтому, значение по умолчанию равняется 150.

Эти значения f107Daily соответствуют радио-потоку на 10,7 см на фактическом расстоянии от Земли от Sun. Значения f107Daily не соответствуют радио-потоку на уровне 1 а. е. Следующий сайт обеспечивает оба класса значений: ftp://ftp.ngdc.noaa.gov/STP/space-weather/solar-data/solar-features/solar-radio/noontime-flux/penticton/

Смотрите ограничения в Описании для получения дополнительной информации.

flags

Массив 23, чтобы включить или отключить конкретные изменения для выходных параметров. Если длина массива flags, m, равняется 23, и вы не задали все доступные входные параметры, эта функция принимает, что flags установлен.

Флаги, сопоставленные с входом flags, включают или отключают конкретные изменения для выходных параметров:

Поле	Описание
`Flags(1)`	Эффект F10.7 на среднее значение
`Flags(2)`	Независимый от времени
`Flags(3)`	Симметричный ежегодник
`Flags(4)`	Симметричный полугодовой
`Flags(5)`	Асимметричный ежегодник
`Flags(6)`	Асимметричный полугодовой
`Flags(7)`	Дневной
`Flags(8)`	Полудневной
`Flags(9)`	Ежедневный AP. Если вы устанавливаете это поле на-1, функция использует целую матрицу магнитной информации об индексе (APH) вместо `APH(:,1)`.
`Flags(10)`	Весь UT, продольные эффекты
`Flags(11)`	Продольный
`Flags(12)`	UT и смешанный UT, продольный
`Flags(13)`	Смешанный AP, UT, продольный
`Flags(14)`	Terdiurnal
`Flags(15)`	Отклонения от распространяющегося равновесия
`Flags(16)`	Все экзосферические температурные изменения
`Flags(17)`	Все изменения от 120 000-метровой температуры (TLB)
`Flags(18)`	Вся более низкая термосфера (TN1) изменения температуры
`Flags(19)`	Весь 120 000-метровый градиент (S) изменения
`Flags(20)`	Вся верхняя стратосфера (TN2) изменения температуры
`Flags(21)`	Все изменения от 120 000 показаний счетчиков (ZLB)
`Flags(22)`	Вся более низкая температура мезосферы (TN3) изменения
`Flags(23)`	Изменения шкалы высот Turbopause

latitude

Массив m-by-1 геодезические широты, в градусах.

longitude

Массив m-by-1 долготы, в градусах.

localApparentSolarTime

Массив m-by-1 локальное очевидное солнечное время (часы). Чтобы получить физически реалистическое значение, функция устанавливает localApparentSolarTime на (sec/3600 + lon/15) по умолчанию. См. Описание для получения дополнительной информации.

magneticIndex

Массив m-by-7 магнитной информации об индексе. Если вы задаете magneticIndex, необходимо также задать f107Average и f107Daily. Эта информация состоит из:

Ежедневный магнитный индекс (AP)

3-часовой AP в течение текущего времени

3-часовой AP в течение 3 часов перед текущим временем

3-часовой AP в течение 6 часов перед текущим временем

3-часовой AP в течение 9 часов перед текущим временем

Среднее значение восьми 3-часовых индексов AP с 12 до 33 часов перед текущим временем

Среднее значение восьми 3-часовых индексов AP с 36 до 57 часов перед текущим временем

Эффекты ежедневного магнитного индекса не являются большими или установленными ниже 80 000 м. В результате функция устанавливает значение по умолчанию на 4. Смотрите ограничения в Описании для получения дополнительной информации.

otype

Общая массовая плотность вывела:

`'Oxygen'`	Общая массовая плотность выходные параметры включает аномальную кислородную плотность номера.
`'NoOxygen'`	Общая массовая плотность выходные параметры не включает аномальную кислородную плотность номера.

UTseconds

Массив m-by-1 секунды в день во всемирное время (UT)

year

Эта функция игнорирует значение year.

Выходные аргументы

T

Массив N-by-2 значения температуры, в кельвине. Первый столбец является экзосферической температурой в кельвине. Второй столбец является температурой на высоте в кельвине.

rho

Массив N-by-9 значения плотности (^кг/м3 или ^1/m3) в выбранных модулях плотности. Порядок следования столбцов:

Плотность его, в ^1/m3

Плотность O, в ^1/m3

Плотность N2, в ^1/m3

Плотность O2, в ^1/m3

Плотность площади, в ^1/m3

Общая массовая плотность, в ^кг/м3

Плотность H, в ^1/m3

Плотность N, в ^1/m3

Аномальная кислородная плотность номера, в ^1/m3

density(6), общая массовая плотность, является суммой массовой плотности Его, O, N2, O2, Площади, H, и N. Опционально, density(6) может включать массовую плотность аномального кислорода, делающего density(6), эффективную общую массовую плотность для перетаскивания.

Примеры

Вычислите температуры, плотность не включая аномальный кислород с помощью модели NRLMSISE-00 на уровне 10 000 м, 45 широт степеней,-50 долгот степеней. Это вычисление использует дату 4 января 2007 в 0 UT. Это использует значения по умолчанию для потока, магнитных индексных данных, и локальное солнечное время с действиями из области значений, генерирующими предупреждения:

[T, rho] = atmosnrlmsise00( 10000, 45, -50, 2007, 4, 0)

T =
   1.0e+03 *

    1.0273    0.2212

rho =

   1.0e+24 *

    0.0000         0    6.6824    1.7927    0.0799    0.0000         0         0         0

Вычислите температуры, плотность не включая аномальный кислород с помощью модели NRLMSISE-00. Используйте модель на уровне 10 000 м, 45 широт степеней, –50 долгот степеней и 25 000 м, 47 широт степеней, –55 долгот степеней.

Это вычисление использует дату 4 января 2007 в 0 UT. Это использует значения по умолчанию для потока, магнитных индексных данных, и локальное солнечное время с действиями из области значений, генерирующими предупреждения:

[T, rho] = atmosnrlmsise00( [10000; 25000], [45; 47], ...
[-50; -55], [2007; 2007], [4; 4], [0; 0])

T =

  1.0e+003 *

    1.0273    0.2212
    1.0273    0.2116


rho =

  1.0e+024 *

    0.0000     0    6.6824    1.7927    0.0799    0.0000     0     0     0
    0.0000     0    0.6347    0.1703    0.0076    0.0000     0     0     0

Вычислите температуры, плотность включая аномальный кислород с помощью модели NRLMSISE-00 на уровне 10 000 м, 45 широт степеней, –50 долгот степеней. Это вычисление использует дату 4 января 2007 в 0 UT. Это использует значения по умолчанию для потока, магнитных индексных данных, и локальное солнечное время с действиями из области значений, генерирующими ошибки:

[T, rho] = atmosnrlmsise00( 10000, 45, -50, 2007, ...
4, 0, 'Oxygen', 'Error')

T =

  1.0e+003 *

    1.0273    0.2212


rho =

  1.0e+024 *

    0.0000     0    6.6824    1.7927    0.0799    0.0000     0     0     0

Вычислите температуры, плотность включая аномальный кислород с помощью модели NRLMSISE-00 на уровне 100 000 м, 45 широт степеней, –50 долгот степеней. Это вычисление использует дату 4 января 2007 в 0 UT. Это использует заданные значения для потока и магнитные индексные данные, и локальное солнечное время по умолчанию. Это указывает, что действие из области значений не должно генерировать сообщение:

aph = [17.375 15 20 15 27 (32+22+15+22+9+18+12+15)/8 (39+27+9+32+39+9+7+12)/8]
f107 = 87.7
nov_6days  = [ 78.6 78.2 82.4 85.5 85.0 84.1]
dec_31daymean = 84.5
jan_31daymean = 83.5
feb_13days = [ 89.9 90.3 87.3 83.7 83.0 81.9 82.0 78.4 76.7 75.9 74.7 73.6 72.7]
f107a = (sum(nov_6days) + sum(feb_13days) + (dec_31daymean + jan_31daymean)*31)/81
flags = ones(1,23)
flags(9) = -1
[T, rho] = atmosnrlmsise00( 100000, 45, -50, 2007, 4, 0, f107a, f107, ...  
aph, flags, 'Oxygen', 'None')
aph =

   17.3750   15.0000   20.0000   15.0000   27.0000   18.1250   21.7500


f107 =

   87.7000


nov_6days =

   78.6000   78.2000   82.4000   85.5000   85.0000   84.1000


dec_31daymean =

   84.5000


jan_31daymean =

   83.5000


feb_13days =

  Columns 1 through 10

   89.9000 90.3000 87.3000 83.7000 83.0000 81.9000 82.0000 78.4000 76.7000 75.9000

  Columns 11 through 13

   74.7000   73.6000   72.7000


f107a =

   83.3568


flags =

  Columns 1 through 17

     1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1

  Columns 18 through 23

     1  1  1  1  1  1


flags =

  Columns 1 through 17

     1  1  1  1  1  1  1  1 -1  1  1  1  1  1  1  1  1

  Columns 18 through 23

     1  1  1  1  1  1


T =

  1.0e+003 *

    1.0273    0.1917


rho =

  1.0e+018 *

    0.0001  0.4241  7.8432  1.9721  0.0808  0.0000  0.0000  0.0000  0.0000

Ссылки

https://ccmc.gsfc.nasa.gov/modelweb/atmos/nrlmsise00.html

Смотрите также

atmoscira

Документация