atmosnrlmsise00

Implement mathematical representation of 2001 United States Naval Research Laboratory Mass Spectrometer and Incoherent Scatter Radar Exosphere

Syntax

[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)

Description

[T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds) implements the mathematical representation of the 2001 United States Naval Research Laboratory Mass Spectrometer and Incoherent Scatter Radar Exosphere (NRLMSISE-00) of the MSIS^® class model. NRLMSISE-00 calculates the neutral atmosphere empirical model from the surface to lower exosphere (0 m to 1,000,000 m). Optionally, it performs this calculation including contributions from anomalous oxygen that can affect satellite drag above 500,000 m.

[T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, localApparentSolarTime) specifies an array of m local apparent solar time (hours).

[T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, f107Average, f107Daily, magneticIndex) specifies arrays of m 81 day average of F10.7 flux (centered on doy), m-by-1 daily F10.7 flux for previous day, and m-by-7 of magnetic index information.

[T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, flags) specifies an array of 23 to enable or disable particular variations for the outputs.

[T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, otype) specifies a character vector or string for total mass density output.

[T rho] = atmosnrlmsise00(altitude, latitude, longitude, year, dayOfYear, UTseconds, action) specifies out-of-range input action.

This function has the limitations of the NRLMSISE-00 model. For more information, see the NRLMSISE-00 model documentation.

The NRLMSISE-00 model uses UTseconds, localApparentSolarTime, and longitude independently. These arguments are not of equal importance for every situation. For the most physically realistic calculation, choose these three variables to be consistent by default:

 localApparentSolarTime = UTseconds/3600 + longitude/15

If available, you can include departures from this equation for localApparentSolarTime, but they are of minor importance.

Input Arguments

action

Action for out-of-range input. Specify one:

'Error'

'Warning' (default)

'None'

altitude

Array of m-by-1 altitudes, in meters.

dayOfYear

Array m-by-1 day of year.

f107Average

Array of m-by-1 81 day average of F10.7 flux (centered on day of year (dayOfYear)). If you specifyf107Average, you must also specify f107Daily and magneticIndex. The effects of f107Average are not large or established below 80,000 m; therefore, the default value is 150.

These f107Average values correspond to the 10.7 cm radio flux at the actual distance of the Earth from the Sun. The f107Average values do not correspond to the radio flux at 1 AU. The following site provides both classes of values: ftp://ftp.ngdc.noaa.gov/STP/space-weather/solar-data/solar-features/solar-radio/noontime-flux/penticton/

See the limitations in Description for more information.

f107Daily

Array of m-by-1 daily F10.7 flux for previous day. If you specify f107Daily, you must also specify f107Average and magneticIndex. The effects of f107Daily are not large or established below 80,000 m; therefore, the default value is 150.

These f107Daily values correspond to the 10.7 cm radio flux at the actual distance of the Earth from the Sun. The f107Daily values do not correspond to the radio flux at 1 AU. The following site provides both classes of values: ftp://ftp.ngdc.noaa.gov/STP/space-weather/solar-data/solar-features/solar-radio/noontime-flux/penticton/

See the limitations in Description for more information.

flags

Array of 23 to enable or disable particular variations for the outputs. If flags array length, m, is 23 and you have not specified all available inputs, this function assumes that flags is set.

The flags, associated with the flags input, enable or disable particular variations for the outputs:

Field	Description
`Flags(1)`	F10.7 effect on mean
`Flags(2)`	Independent of time
`Flags(3)`	Symmetrical annual
`Flags(4)`	Symmetrical semiannual
`Flags(5)`	Asymmetrical annual
`Flags(6)`	Asymmetrical semiannual
`Flags(7)`	Diurnal
`Flags(8)`	Semidiurnal
`Flags(9)`	Daily AP. If you set this field to -1, the function uses the entire matrix of magnetic index information (APH) instead of `APH(:,1)`.
`Flags(10)`	All UT, longitudinal effects
`Flags(11)`	Longitudinal
`Flags(12)`	UT and mixed UT, longitudinal
`Flags(13)`	Mixed AP, UT, longitudinal
`Flags(14)`	Terdiurnal
`Flags(15)`	Departures from diffusive equilibrium
`Flags(16)`	All exospheric temperature variations
`Flags(17)`	All variations from 120,000 meter temperature (TLB)
`Flags(18)`	All lower thermosphere (TN1) temperature variations
`Flags(19)`	All 120,000 meter gradient (S) variations
`Flags(20)`	All upper stratosphere (TN2) temperature variations
`Flags(21)`	All variations from 120,000 meter values (ZLB)
`Flags(22)`	All lower mesosphere temperature (TN3) variations
`Flags(23)`	Turbopause scale height variations

latitude

Array of m-by-1 geodetic latitudes, in degrees.

longitude

Array of m-by-1 longitudes, in degrees.

localApparentSolarTime

Array of m-by-1 local apparent solar time (hours). To obtain a physically realistic value, the function sets localApparentSolarTime to (sec/3600 + lon/15) by default. See Description for more information.

magneticIndex

An array of m-by-7 of magnetic index information. If you specify magneticIndex, you must also specify f107Average and f107Daily. This information consists of:

Daily magnetic index (AP)

3 hour AP for current time

3 hour AP for 3 hours before current time

3 hour AP for 6 hours before current time

3 hour AP for 9 hours before current time

Average of eight 3 hour AP indices from 12 to 33 hours before current time

Average of eight 3 hour AP indices from 36 to 57 hours before current time

The effects of daily magnetic index are not large or established below 80,000 m. As a result, the function sets the default value to 4. See the limitations in Description for more information.

otype

Total mass density output:

`'Oxygen'`	Total mass density outputs include anomalous oxygen number density.
`'NoOxygen'`	Total mass density outputs do not include anomalous oxygen number density.

UTseconds

Array of m-by-1 seconds in day in universal time (UT)

year

This function ignores the value of year.

Output Arguments

T

Array of N-by-2 values of temperature, in kelvin. The first column is exospheric temperature, in kelvin. The second column is temperature at altitude, in kelvin.

rho

An array of N-by-9 values of densities (kg/m³ or 1/m³) in selected density units. The column order is:

Density of He, in 1/m³

Density of O, in 1/m³

Density of N2, in 1/m³

Density of O2, in 1/m³

Density of Ar, in 1/m³

Total mass density, in kg/m³

Density of H, in 1/m³

Density of N, in 1/m³

Anomalous oxygen number density, in 1/m³

density(6), total mass density, is the sum of the mass densities of He, O, N2, O2, Ar, H, and N. Optionally, density(6) can include the mass density of anomalous oxygen making density(6), the effective total mass density for drag.

Examples

Calculate the temperatures, densities not including anomalous oxygen using the NRLMSISE-00 model at 10,000 m, 45 degrees latitude, -50 degrees longitude. This calculation uses the date January 4, 2007 at 0 UT. It uses default values for flux, magnetic index data, and local solar time with out-of-range actions generating warnings:

[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

Calculate the temperatures, densities not including anomalous oxygen using the NRLMSISE-00 model. Use the model at 10,000 m, 45 degrees latitude, –50 degrees longitude and 25,000 m, 47 degrees latitude, –55 degrees longitude.

This calculation uses the date January 4, 2007 at 0 UT. It uses default values for flux, magnetic index data, and local solar time with out-of-range actions generating warnings:

[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

Calculate the temperatures, densities including anomalous oxygen using the NRLMSISE-00 model at 10,000 m, 45 degrees latitude, –50 degrees longitude. This calculation uses the date January 4, 2007 at 0 UT. It uses default values for flux, magnetic index data, and local solar time with out-of-range actions generating errors:

[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

Calculate the temperatures, densities including anomalous oxygen using the NRLMSISE-00 model at 100,000 m, 45 degrees latitude, –50 degrees longitude. This calculation uses the date January 4, 2007 at 0 UT. It uses defined values for flux, and magnetic index data, and default local solar time. It specifies that the out-of-range action is to generate no message:

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